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TTB  OBJECTIVES 


o ASSESS  PROPULSION  TECHNOLOGY  ADVANCES 
IN  AN  ENGINE  SYSTEMS  ENVIRONMENT 

o ENHANCE  THE  PROCESS  FOR  IMPLEMENTATION  OF  TECHNOLOGY 
INTO  EMERGING  AND  OPERATIONAL  PROGRAMS 

o PROVIDE  SYSTEM  TEST  CAPABILITY  FOR  EVALUATION  OF 
PROTOTYPE  HARDWARE 

o SUPPORT  NASA  PROGRAMS  WITH  ANOMALY  RESOLUTION  ON  AN 
AS-NEEDED  BASIS 

o DEVELOP  AND  MAINTAIN  IN-HOUSE,  HANDS-ON  ROCKET  PROPULSION 
HARDWARE  AND  TEST  EXPERIENCE/CAP  ABILITY  AT  MSFC 
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TTB  MAJOR  ACCOMPLISHMENTS 


o 31  SSME  TESTS  CONDUCTED/  -3000  SECONDS  CUMULATIVE  TEST  TIME 
o EVALUATION  OF  A MODIFIED  SSME 

o ENGINE  3001  (HIGHLY  INSTRUMENTED  SSME)  ENVIRONMENT  CHARACTERIZATION 
o DEMONSTRATION  OF  SSME  ADVANCED  DEVELOPMENT  CONCEPTS 
o SUPPORT  OF  SHUTTLE  FLIGHT  AND  DEVELOPMENT  INVESTIGATIONS 
o ASSESSMENT  OF  ADVANCED  PROPULSION  TECHNOLOGY  CONCEPTS 
o NUMEROUS  IMPROVEMENTS  DERIVED  FROM  IN-HOUSE,  HANDS-ON  INVOLVEMENT 
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ACCOMPLISHMENT:  LARGE  THROAT  SSME  CHARACTERIZATION 

SYSTEM  PERFORMANCE  DEFINED 
MORE  BENIGN  ENVIRONMENT  DEMONSTRATED 
STABLE  COMBUSTION  DEMONSTRATED 
TURBOMACHINERY  REDESIGN  REQUIREMENTS  IDENTIFIED 

ENGINE:  0208 

BENEFIT : EARLY  ASSESSMENT  OF  PROPOSED  SSME  IMPROVEMENTS 

PERFORMANCE  OF  ASSESSMENT  WAS  FREE  FROM  USUAL  DEVELOPMENT 
PRESSURES  AND  WAS  ACCOMPLISHED  WITHOUT  INTERFERENCE  TO  SSC 

IMPACT : TTB  RESULTS  FORMED  THE  BASIS  FOR  ADVANCED  FAB  LTMCC  BASELINE 

STABLE  COMBUSTION  WITHOUT  BAFFLES  OR  ACOUSTIC  CAVITIES  IS 
APPLICABLE  TO  FUTURE  DESIGNS 
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ACCOMPLISHMENT:  PHASE  H SSME  INTERNAL  ENVIRONMENT 

A CHARACTERIZATION 

~ 630  TTB-UNIQUE  ENGINE  MEASUREMENTS 

DEFINES  INTEGRATION  PARAMETERS  FOR 
ALTERNATE  COMPONENTS  (e.g.,  ATDs) 


00 


ENGINE:  3001 


BENEFIT:  INCREASED  UNDERSTANDING  OF  OPERATING  ENVIRONMENT 

BASELINE  DATA  FOR  COMPARISON  WITH  ATDs  AND  PHASE  B+  POWERHEAD 


IMPACT: 


DIGITAL  TRANSIENT  MODEL 

TO  BE  FULLY  REALIZED.  EXAMPLES:  SIGNIFICANT  DEKICENCIES  IN  SSME 
POWER  BALANCE  MODEL  HAVE  BEEN  DISCOVERED.  HAVE  MADE  NEW 
DISCOVERIES  ABOUT  TURBOPUMP  OPERATION. 


George  C.  Marshall  Space  Flight  Center 

Science  and  Engineering  Directorate 

Propulsion  Laboratory 


MAJOR  IMPACTS  OF  3001  TESTING 
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o ACQUISITION  OF  HERETOFORE  UNMEASURED  SSME  HOT-FERE  DATA  - COMPLETE 
MAPPING  OF  ENGINE  OPERATION  FOR  RANGE  OF  MIXTURE  RATIO,  POWER  LEVEL, 
PUMP  INLET  (NPSP),  F7  ORIFICE  AND  REPRESS  CONDITIONS 

MAJOR  FLOWRATES 
INSTRUMENTED  TURBOPUMPS 
OTHER  Ps,  Ts,  STRAINS,  etc. 


o CALIBRATION  AND  IMPROVEMENT  OF  MODELS 

POWER  BALANCE  MODEL 
1-D  TURBOPUMP  MODELS 
CFD  MODELS 
THERMAL  MODELS 
STRUCTURAL  MODELS 
STRESS  MODELS 


o NUMEROUS  LESSONS  LEARNED 

INSTRUMENTATION  DESIGN  AND  IMPLEMENTATION 
TEST  OPERATION  EFFICIENCY  ENHANCEMENTS 


~ 630  MEASUREMENTS 
IN  TOTAL 
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SOME  CHALLENGES  TO  CFD 


o COMPARE  HOT-FIRE  PREDICTIONS  WITH  3001  DATA 

- UTILIZE  DATA  TO  CALIBRATE/IMPROVE  MODELS 

- UTILIZE  MODELS  (+  KNOWLEDGE  OF  FLUID  MECHANICS,  etc.) 
TO  EXPLAIN  DATA 


o TAKE  AN  ACTIVE  ROLE  AND  ENCOURAGE  CONCURRENT  ENGINEERING 

IN  THE  DESIGN  PROCESS 

- PLACEMENT  OF  SENSORS  IN  INSTRUMENTED  TEST  ARTICLES 

- HARDWARE  DESIGN  (EARLY  INVOLVEMENT,  TIMELY  INPUT) 


o 


EXAMPLES ... 
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3001  LESSONS  LEARNED  (CONT.) 
DESIGN  AND  OPERATIONS 


• ENGINE  DESIGN/TESTING 

INSTRUMENTATION: 

• GLASS  BRAIDED  TYPE  THERMOCOUPLE  WIRE  SHOULD  BE 
REPLACED  BY  TEFLON  COATED  T/C  WIRE  TO  AVOID 
UNRAVELING  OF  COATING  DUE  TO  NORMAL  FIELD 
OPERATION,  (i.e.  RUBBING  WIRES  DURING  PUMP 
INSTALLATION) 

• ADDITIONAL  EFFORT  SHOULD  BE  MADE  ON  INSTRUMENTATION 
PLACEMENT  TO  AVOID  LOCATIONAL  EFFECTS  WHICH 
MISREPRESENT  THE  DESIRED  DATA:  (i.e.  ELBOWS, 

STAGNATION  AND  RECIRCULATION  REGIONS) 

• EFFORT  SHOULD  BE  MADE  TO  EQUALIZE  SENSE  LINE  VOLUMES 
ON  BOTH  HIGH  AND  LOW  SIDE  OF  DELTA  PRESSURE 
TRANSDUCERS  WHERE  ACCURATE  TRANSIENT  DATA  IS  REQUIRED 

• PRESSURE  SENSE  LINES  REQUIRED  TO  MEASURE  HOT  GAS 
ENVIRONMENTS  AND  THAT  ARE  ROUTED  ’ NEAR  ’ COOLANT 
CIRCUITS  SHOULD  BE  CONSIDERED  FOR  PURGING  DURING 
TESTING  TO  PREVENT  ICING. 


Georgs  C.  Marshall  Space  Flight  Center 

Science  and  Engineering  Directorate 

Propulsion  Laboratory 


TECHNOLOGY  TEST  BED  REVIEW 


f\JASA 

National  Aeronautics  and 
Space  Administration 


FUTURE  PLANS 


CONTINUATION  OF  PHASE  U ENVIRONMENT  CHARACTERIZATION 
o INSTRUMENTED  RKDN  TURBOPUMPS 
o P&W  ALTERNATE  TURBOPUMPS 

CONTINUATION  OF  TECHNOLOGY  ITEM  INTEGRATION  AND  EVALUATION 
o HEALTH  MONITORING  SYSTEMS 
o TURBOMACHINERY 
o COMBUSTION  DEVICES 
o INSTRUMENTATION 
o CONIROLLER 

CHARACTERIZATION  OF  FUTURE  PROPULSION  SYSTEM  DESIGNS 

o SSME  PROTOTYPE  BLOCK  H (PHASE  11+  POWERHEAD,  ATDs,  LTMCC) 
o SSME'  PRODUCEBILITY  IMPROVEMENTS 
o HLLV  STME  PROTOTYPE 


CONFIGURATION 

■ 

□ 

ADVANCED  ENGINE  (0200) 

DOWNTIME  IN  SUPPORT  OF  17*  DISCONNECT 

E3001  W/  STD  RKDN  PUMPS 

DOWNTIME  IN  SUPPORT  OF  PHASE  II  + HGM 

E3001  W/  INSTR  RKDN  PUMPS 

E3001  W/  HYDROSTATIC  BRQ  HPOTP 

E3001  W / INSTR  P&W  FUEL  PUMP  AND 
INSTR  RKDN  LOX  PUMP 

E3001  W/  SALOX  CAGE/SPEED 
SENSOR  HPOTP 

E3001  W/  INSTR  PAW  PUMPS 

E3001  W / SN  BALLS/BEARING 
DEFLECTOMETER  HPOTP 

E 3003-PHASE  HI  CHARACTERIZATION 

n 

TECHNOLOGY  ITEMS: 

i 

EXIT  PLANE  SPECTROSCOPY  OF  SSME 

PLUME  TEMPERATURE  MEASUREMENTS 

■ 

OPTICAL  PLUME  ANOMALY  DETECTOR 

■ 

SODIUM  RESONANT  LINE  ABSORPTION  MEAS. 

SSME  EXIT  PLANE  HOLOGRAPHY 

LASER  INDUCED  FLUORESCENCE 

SAFD 

TTBIHOSC  EXPERT  SYSTEMS 

NOZZLE  OPTIC  ASSEMBLY 

HPOTP  PREBURNER  PUMP  END  HYDROSTATIC 
BEARING  RETROFIT 

LON  COST  CONTROLLER 

■ 

OPTICAL  PROPELLANT  SENSING 

SOLID  STATE  H2/02  SENSORS 

VORTEX  SHEDDING  FLOWMETER 

i 

IMPROVED  BEARING  CAGE  MATERIAL 

FABRY-PEROT  SPECTROMETER 

NON-INTRUSIVE  SPEED  SENSOR 

ASYM.  AND  MAIN  LEAK  DETEC  PRIOR  TO  TEST 

NON-INTRUSIVE  (I.R.)  GAS  TEMP  SENSOR 

SILICON  NITRIDE  BALLS 

BEARING  DEFLECTOMETER 

■ 

HEAT  FLUX  SENSORS 

■ 

VPS  BLADES 

■ 

ULTRASONIC  FLOWMETER 

BRUSHLESS  TORQUEMETER 

■ 

ELECTROMECHANICAL  ACTUATORS 

■ 

PIEZO-ELECTRIC  SENSOR  AUTO-CALIBRATION 

■ 

ADVANCED  MAIN  COMBUSTION  CHAMBER 

George  C.  Marshall  Space  Flight  Center 

Science  and  Engineering  Directorate 

Propulsion  Laboratory 


TECHNOLOGY  TEST  BED  REVIEW 


(VASA 

National  Aeronautics  and 
Space  Administration 


SUMMARY 


TECHNOLOGY  TEST  BED  HAS  PROVEN  TO  BE  EFFECTIVE  IN  CHARACTERIZING 
PROPULSION  SYSTEM  DESIGNS 

TECHNOLOGY  TEST  BED  HAS  PROVIDED  VALUABLE  HOT-FIRE  DATA  FOR 
INCREASED  UNDERSTANDING  OF  THE  INTERNAL  OPERATING  ENVIRONMENT 
OF  THE  SSME  AND  FOR  CALIBRATION  OF  SSME  MODELS 

TECHNOLOGY  TEST  BED  HAS  ENHANCED  SHUTTLE  DEVELOPMENT  TESTING 
AND  ANOMALY  RESOLUTION 

TECHNOLOGY  TEST  BED  HAS  PROVIDED  A VALUABLE  PLATFORM  FOR 
ASSESSMENT  OF  ADVANCED  PROPULSION  TECHNOLOGIES 

TECHNOLOGY  TEST  BED  FUTURE  PLANS  INCLUDE  CHARACTERIZATION  OF 
FORTHCOMING  PROPULSION  SYSTEM  DESIGNS  AND  EVALUATION  OF 
EMERGING  PROPULSION  TECHNOLOGIES 


ADVANCED  SOLID  ROCKET  MOTOR 

PROJECT  STATUS 

Tenth  Workshop  for  Computational  Fluid  Dynamics  (CFD) 

Applications  in  Rocket  Propulsion 


28  April  1992 

Keith  Coates 
EE  71 
MSFC 


Project  Status 


OUTLINE 

■ Project  Objectives 

■ Team 

■ Locations 

■ Motor  Design 

■ Schedule 

■ Technical  Issues 


15/4/10 


ASRM  3074.00 


Project  Status 


PROJECT  OBJECTIVES 


Improve  System  Safety  and  Reliability 

□ Design  Features 

□ Enhanced  Quality 

□ Reproducibility 

Improve  Shuttle  Payload  Performance:  12,000  lb 
Optimize  Program  Cost 

Promote  Competitive  SRM  Industry 

□ Construct  and  Operate  Government  Owned  Manufacturing 
and  Test  installations 


15/4/10 


ASRM  3074.00 


Project  Status 


ASRM  PROJECT  TEAM 


ASRM  TEAM 


LOCKHEED 


RUST 


LOCKHEED 


PROJECT  MGMT 
SE&I 

SUPPORT  EQUIPMENT 
SRM&QA 

FACILITIES  MGMT 


RUST 


• ARCHITECTURE  & 
ENGINEERING  (A&E) 

• CONSTRUCTION 


Project  Status 


ASRM  PROJECT  TEAM  LOCATIONS 


1 S/4/10 


ASRM  X74.00 


Project  Status 


ASRM  MOTOR  DESIGN  HIGHLIGHTS 


INSULATION 

• Uses  RSRM  J-Seal  Design  at  Field 
Joint 

• Castable  Inhibitor-to-Stress  Relief  - 
Flap  Joint  Location  Minimizes 
Potential  for  Jetting  Into  Sidewall 
Insulation 

• Asbestos  Free  Formulation 


NOZZLE 

• Minimum  Joints  and 
Inlet/throat  rings 

• Improved  process 
ablative  materials 

• Eliminates  flexseal  cowl 
and  boot  assembly 


/CASE  C 

• Two  Field  Joints  • e 

• High  Fracture  Toughness  ( 

• High  Stress  Corrosion  Resistance  • 7 

• Welded  Factory  Joints  [ 

• Integral  Stiffeners  and  ET  Attach  Ring 
In  Aft  Segment  Eliminates  Failure  Points 
Experienced  With  Bolt-On  Stiffeners  and  Ring 

PROPELLANT 

• Industry  Proven  HTPB  Propellant  with  over  60M  lbs 
Successfully  Produced 

• Formulated  and  Proven  for  Continuous  Mix  Process 
• Positive  Margins 


GRAIN/IGNITION 

• High  Safety  Margin  In  Forward 
Fin  Trailing  Edge 

• Minimum  Igniter 
Chamber  Bolt  Leak  Paths 

• Expendable  Carbon  Filament 
Chamber 

• TBI  Initiator  Eliminates  S&A 
Device  Leak  Paths 


ASRM  3074.00 


Project  Status 


ASRM  DESIGN  PARAMETERS 


U> 


Diameter/Length,  in 

Average  Thrust  Vacuum,  Ibf;  Web  time 

Delivered  Isp  Vacuum,  sec 

Area  Ratio,  (Ae/At) 

Motor  Weight,  ib 

Propellant  Weight,  Ib 

Motor  Propellant  Mass  Fraction,  (Wp/Wt) 

Inert  Weight , Ib 

□ Metal  Case  Weight/Number  of  Segments,  Ib 

□ Single  Nozzle  Weight,  Ib 

Solid  Propellant  Type 

Average  Chamber  Pressure,  psia;  Action  Time 

Burn  Rate  at  625  psia,  in/sec 

Action  Time,  sec 

Thrust  Vector  Control 

Recovery/Reuse 


150/1,513.43 

2,654,776 

268.1 

7.48 

1,351,092 

1,209,589 

0.895 

141,503 

98,553/3 

18,800 

HTPB 

612 

0.350 

130.9 

Flexible  Bearing 
Yes 


1 S/4/10 


ASRM  3074.00 


ORIGINAL  page  is 

OF  POOR  QUALITY 


Project  Status 


ORIGINAL  PAGE  IS 

OF  POOR  QUALITY 


35 


Project  Status 


CASE  PREP  BUILDING 


36 


Project  Status 


CASE  PREP  AUTOCLAVE 
ON-SITE  FABRICATION  AT  YELLOW  CREEK 


ORIGINAL  PAGE  lit 
OF  POOR  QUALITY 


Project  Status 


MICHOUD  PILING  AND  PERIMETER  WALL 


Project  Status 


VERTICAL  BORING  MILL  (VBM) 


Project  Status 


HEAT  TREAT/CHILLER  BUILDING 


Project  Status 


Project  Status 

ASRM  TECHNICAL  TOP  5 
March  27, 1992 


1.  Mix/Cast  Construction,  Outfitting  & Process  Verification 

2.  Soluble  Casting  Mandrels 


I 

I 


Project  Status 

SOLUBLE  CORE  DESIGN 


43 


APRIL  28, 1992 


Jan  C.  Monk 
Chief  Engineer 

Space  Transportation  Main  Engine 
Marshall  Space  Flight  Center 


N92-32281 


STME  DEFINITION 


1985  11986  11987  11988  11989  1990  1991  1992 

I STME  Phase  A I I | 

I ” II  STME  Advanced  Development  I 


• Over  7500  STME  Trades  and  Studies 

- Engine  System  Level  Trades 

- Engine  Component  Level  Trades 

- Manufacturing  Studies 

- Assembly  Facilities  Studies 

- Test  Facilities  Studies 

- Overhaul  Facilities  Studies 

• STME  Design  Concept  Formed  in  Total  System  Context 

- Interactive  Engine  /Vehicle/Operability  Trades 

- Broad  range  of  Vehicles  and  Missions 

- QFD  established  Strategic  Quality  Characteristics,  Pugh  Concept  Selection 

- National  Consensus  Cycle  Decision 

• Design  Effort  Supported  by  Extensive  Data  Base 

- ADP  Hardware 

- Lessons  Learned:  SSME,  J-2,  LR87,  F-100,  F-119,  etc. 


DESIGN  PHILOSOPHY 


• ROBUST  DESIGN 

• DESIGN  FOR  OPERABILITY 

• DESIGN  FOR  RELIABILITY 

• DESIGN  FOR  LOW  COST 


ROBUST  DESIGN  g 

• ACHIEVABLE  REQUIREMENTS 

- Mission  Life  - 10 

- Chamber  Pressure  - 2250  psia 

- Weight  - 9100  pounds 

• SERIES  TURBINES 

- Enhances  benign  system  response  to  fuel  turbopump  failures 

• DESIGN  MARGINS 

- Design  based  on  internal  operating  environment  worst  case 
plus  development  margin 

• REDUCED  INTERNAL  ENVIRONMENTS  TO  ASSURE 
ROBUST  PROCESSES 

- No  sheet  metal  liners,  no  overlays,  no  platings 

• NEAR  NET  SHAPE  PROCESSES 

- Minimum  number  of  welds,  capable  processes,  reduced 
process  steps 


eisd  - 


MAXIMUM  DESIGN  CONDITION  || 


• MDC  = NOMINAL  PARAMETER  VALUE 

+ RSS  OF  PREDICTED  COMPONENT  VARIATIONS 
+ THRUST/MR  CALIBRATION  ERRORS 
+ FLIGHT  EFFECTS 
+ DEVELOPMENT  MARGIN 

FUEL  PUMP  DISCHARGE  PRESSURE 
48001 1 


4600 

4400 

10%  Development  Margin 

4200 

Flight  Effects 

F & MR  Tolerance 


CASTINGS  VS  MACHINED  AND 
WELDED  FORGINGS 


Material  costs  $ 15,182 

Machine  and  welding  440,824. 


Total  cost 


$456,006 


SSME  Turbopump  Volute 


1-7 1 8 cast 
Machining 

Total  cost 


$14,527= 

20.700 


IR&D  Cast  Volute 


* Based  on  vendor  quote  for  6 parts 


OPERABILITY  CONSIDERATIONS 


• TANK-HEAD  START 

- No  start  valve  or  ducting 

• ELECTROMECHANICAL  ACTUATORS 

- Simplify  pre-launch  checkout 

- No  pneumatics  or  hydraulics 

• COMMONALITY  FACILITATES  INVENTORY  AND 
ON-PAD  REPLACEMENT 

- Electromechanical  actuators/seals/fastners 

• REDUCTION  OF  POTENTIAL  LEAK  PATHS 

- Proven  joint  and  seal  designs 

- Reduced  number  of  joints 

• SIMPLE  COMPONENT  INSTALLATION 

- No  stretch  bolt  joints 

• AUTOMATED  PRE-LAUNCH  CONTROLS  CHECKOUT 

• MINIMUM  FLIGHT-CRITICAL  INSTRUMENTATION 


HIGH  RELIABILITY  DESIGN  PHILOSOPHY 


• DESIGN  FOR  RELIABILITY  IS  AN  INTEGRAL  PART 
OF  THE  STME  DESIGN  PROCESS 

- Concurrent  engineering 

- Design  to  reliability  goals 

- reliability  lessons  learned 

- Bottoms-up  failure  modes  and  effects  analysis 

- Tops-down  fault  tree  analysis 

- Reliability  tracking 

- New  manufacturing  techniques 

- Reliability  demonstration  program 

- Development  margin 


ENGINE  RELIABILITY  ENHANCEMENT 


• NON-INTRUSIVE  OXIDIZER  HEAT  EXCHANGER 
(FOR  TANK  PRESSURIZATION) 

• SIMPLE  SAFETY  MONITORING  APPROACH.  SAFE 
ENGINE  SHUTDOWN  TRIGGERED  BY  ABNORMAL 
VALUES  OF: 

- Chamber  pressure 

- Interpropellant  seal  purge  pressure 

- Gas  generator  temperature 

• MECHANICALLY  LINKED  GAS  GENERATOR  VALVES 
(TO  PREVENT  GG  MIXTURE  RATIO  EXCURSIONS) 

• PRUDENT  USE  OF  REDUNDANCY  IN  CONTROL  SYSTEM 

- Dual  EMA  motors  and  resolvers 

- Duplex  / triplex  controller  electronics 

- Dual  power  source 


• DESIGN-TO-COST  IS  AN  INTEGRAL  PART  OF 
THE  STME  DESIGN  PROCESS 

- New  manufacturing  techniques 

- Advanced  Development  Programs  used  to  investigate 

low  cost  ideas 

- Suppliers  integral  part  of  design  effort 

- Customer  (Government)  integral  part  of  design  effort 

- Costs  continually  estimated  and  tracked 

- Cost  drivers  identified  and  worked 

- Trade  studies  used  to  select  lowest  cost  concepts 

- Zero  RID'S 

- Zero  MR'S 


Ln 

Ln 


n\m\m 


rSUME 


ENGINE  SYSTEM  REQUIREMENTS 


. . ^ g iaaste-  ^ ^ ^ « 


Dual  Thrust  (Step) 

Normal  Thrust  (100%): 
Minimum  Thrust  (70% ): 
Mixture  Ratio: 

Specific  Impulse  (100%): 
Specific  Impulse  (70%): 
Chamber  Pressure: 

Dry  Weight: 

Area  Ratio: 

Length: 

Nozzle  Diameter: 

Design  Life: 

Verified  Reliability: 

500th  Unit  Cost  Goal: 


650.000  Ibf 

455.000  Ibf 
6.0 

428.5  sec. 
427.3  sec. 
2250  psia 
9100  lbm 
45:1 

160  inches 
96  inches 
10  missions 
0.995 
$5.3M 


P = 45 
T = 37.5 
tn  =217 


N = 24404 1 


P = 3969 
T = 73 


Rated  Conditions 


FimI 

Repress 


GOX 

Repress 


ilJ 

L 

L 

■ 

HEX 

, \ 

_ 

n 

P = 2232 
T = 1600 
m = 74 


Ul 

GGFV  GGOV 


iOTBP 


Internal 

Mixer 


AP  = 1109 


Pc  = 2250 


■P  = 3313 
C> 

■ T=  183 


lllf 


F (vac)  = 650Klb 
Isp  (vac)  = 428.5 
MR  =6.0 
EPS  =45 


P-295 
T = 1190 
m = 74 


§dm<2®  & 


FUEL  TURBOPUMP 


Integrally  Cast 
Volute/Diffuser 
Main  Housing 


Cast  Ti  Impellers 


One  Piece 

Manifold 

Casting 


- Bolt-In 
Nozzle 


' Stator  and  Exit 
Guide  Vanes 


Load  Sharing  - 
Damping  Seals 

High 

Cavitation  N 

Margin 

Inducer 


m 


l Stub  Shaft 
Allows  Removal 
Of  Turbine 
Components 


■ Low  Cost 
Rotor  Blisks 
HEE  Resistant 


Transient  ' j 

Thrust  Control  / 
Bearing  / 

High  Capacity  — ^ 
Hydrostatic  Bearing 
Control  s Rotordynamics 


High  Capacity 
Axial  Thrust 
Balance  Piston 


'Cast 

Inlet  Housing 


'Aluminum 

Crossovers 


- Exhaust 
Elbow 


Lift-Off  Seal 

Eliminates 

Purge 


FUEL  TURBOPUMP 


BASELINE  DESIGN  FEATURES 

• CAST  HOUSINGS 

• INTEGRALLY  BLADED  TURBINE 

• CAST  PUMP  IMPELLERS 

• HYDROSTATIC  BEARINGS 

OPTIONAL  FEATURES 

• LOW  COST  MACHINED  IMPELLERS 

• BACKUP  DESIGNS  FOR  TURBINE  AND  BEARINGS 

DEVELOPMENT  STATUS 

• FABRICATION  PROCESSES  IN  TRIAL 

• HYDROSTATIC  BEARING  RIG  TESTING  UNDERWAY 

• SUPPLIERS  ON  CONCURRENT  ENGINEERING  TEAM 

92  PLANS 

• CONTINUED  FAB  DEVELOPMENT 

• HYDROSTATIC  BEARING  TESTING 


LEAD:  STPT-ROCKETDYNE/CANOGA  PARK  CA, 
AEROJET/SACRAMENTO  CA 


Discriminators 

SSME 

STME 

Number  of  parts 

198 

60 

Number  of  weld  counts 

801 

1 

Weld  overlays 

Yes 

No 

Sheet  Metal  Liners 

Yes 

No 

Cost  (TFU) 

$4.1M 

$ 1.02M 

Discharge  Pressure 

6320  psia 

3969  psia 

Turbine  Inlet  Temperature 

1860  deg  R 

1600  deg  R 

Speed 

35180  rpm 

23000  rpm 

Weight 

771  lbs 

1877  lbs 

§d®nn<c©  <&  EnngkKMirtiDg 


LOX  TURBOPUMP 


Vaneless  diffuser  mixes 
and  equalizes  inlet  flow 


Volutes  reduce  turbine  - 
side  loads  and  minimize 
flow  disturbances 


Materials 
• Low  Cost 


Radial  Vanes 
Aid  in  thrust 
balance  by 
equalizing 
pressure 


• Hollow 

• Damped 

• Unshrouded 

Rotordvnamics 

• Subcritical 

• Damper  Seals 


Axial 

Pump 

Inlet 


Cast  Impeller  — 
Low  Tip  Speed 


Splitter 
to  reduce 
side  loads 


Cast  Housings  ' 

• Turbine  & Pump 


Bearings 
• Same  at  Both 
Locations 


LOX  TURBOPUMP 
1 — 


BASELINE  DESIGN  FEATURES 

• CAST  HOUSINGS 

• SUB-CRITICAL  SPEED  DESIGN 

• CAST  PUMP  IMPELLORS 

• SINGLE  PIECE  FORGED  DISK/SHAFT 

OPTIONAL  FEATURES 

• LOW  COST  MACHINED  IMPELLORS 

• BACKUP  DESIGNS  FOR  SHAFT  AND  BEARINGS 

DEVELOPMENT  STATUS 

• FABRICATION  PROCESSES  IN  TRIAL 

• SPIN  RIG  INSTRUMENTATION  UNDERWAY 

• SUPPLIERS  ON  CONCURRENT  ENGINEERING  TEAM 

92  PLANS 

• CONTINUED  FAB  DEVELOPMENT 

• SPIN  TESTING  OF  IMPELLER 

LEAD:  STPT-PRATT  & WHITNEY/WEST  PALMBEACH  FL 


Discriminators 

SSME 

STME 

Number  of  parts 

153 

87 

Number  of  weld  counts 

128 

0 

Weld  overlays 

Yes 

No 

Disk  gold  plating 

Yes 

No 

Cost  (TFU) 

5.16M 

1.38M 

Discharge  Pressure 

4360/7340  psia 

3313  psia 

Turbine  Inlet  Temperature 

1510  deg  R 

1190  deg  R 

Speed 

28200  rpm 

7913  rpm 

Weight 

570  lbs 

1712  lbs 

§dl(Mi<£©  <&  IEimgki©sirfiing 


MAIN  INJECTOR 


BASELINE  DESIGN  FEATURES 

• TWO  MAJOR  COMPONENTS  (DOME  AND  BODY) 

• BRAZED  IN  ELEMENTS 

• 100  % INSPECTIBLE 


OPTIONAL  FEATURES 

• ELEMENTS  INTEGRAL  TO  BODY  (SINGLE  PIECE) 
DEVELOPMENT  STATUS 

• FABRICATION  PROCESSES  IN  TRIAL  DEVELOPMENT 

• 40K  SUBSCALE  AND  UNIELEMENT  TESTING  INITIATED 

• SUPPLIERS  ON  CONCURRENT  ENGINEERING  TEAM 

• FIRST  FULL  SCALE  TEST  INJECTOR  DELIVERED 

92  PLANS 

• PATTERN  OPTIMIZATION 

• 40K  TEST  FOR  PERFORMANCE  VERIFICATION 

• 1ST  FULL  SCALE  TEST  (STABILITY  DATA) 

LEAD:  STPT-AEROJET/SACRAMENTO  CA 


. hfjh 


Sdtenn©©  & Enng8nn®®irfiing 


Discriminators 


Number  of  parts 
Number  of  processes 
Number  of  welds 
Number  of  inspections 
Cost  (TFU) 

Chamber  Pressure 
Weight 


SSME 


90 
2.71M 
3126  psia 
394  lbs 


STME 


2213 
151 
13 
58 
0.88M 
2250  psia 
1339  lbs 


COMBUSTION  CHAMBER 
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COMBUSTION  CHAMBER 


BASELINE  DESIGN  FEATURES 

• CAST  STRUCTURAL  JACKET  AND  MANIFOLD 

• VACUUM  PLASMA  SPRAY  NARLOY  Z LINER 

OPTIONAL  PROCESSES 

• LIQUID  INTERFACE  DIFFUSION  BOND  ASSY 

• PLATELET  LINER 

DEVELOPMENT  STATUS 

• FABRICATION  PROCESSES  IN  TRIAL 

• 1ST  FAB  DEVELOPMENT  JACKET  CAST  @ PCC  -PORTLAND,  OR 

• 1ST  LIDB  CHAMBER  FAB  COMPLETE 

• VPS  SPRAY  TESTS  UNDERWAY 

92  PLANS 

• LIDB  CHAMBER  HOT  FIRE  TEST  WITH  LSI 

• CONTINUE  CASTING  AND  VPS  DEVELOPMENT 

LEAD:  STPT-ROCKETDYNE/CANOGA 
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NOZZLE 


CAST 

NOZZLE  COOLANT 
MANIFOLD 


STIFFENERS 


i ATTACHMENT 

L CONE  — * 

VJ 


CONVECTIVE  COOLANT 
GAS  OUT  TO  AMBIENT 


TUBULAR  SKIRT  / JACKET 

INFLATION  FORMED  / 
l~~~  heat  SHIELD  DIFFUSION  BONDED  - 
/ TUBES,  JACKET 


SECTION  A-A 

• ALL  MATERIALS  ARE  INCO  625 


SECONDARY 


PRIMARY  COOLANT  FILM 
INJECTION 


NOZZLE | 


DESIGN  FEATURES 

• ALL  COMPONENTS  MADE  OF  SAME  ALLOY 

• TUBULAR  SKIRT/JACKET  - NO  WELDS 

• BOLT  TO  CHAMBER 

• CAST  MANIFOLDS 

DEVELOPMENT  STATUS 

• FABRICATION  PROCESSES  IN  TRIAL 

• BASELINE  TO  BE  SELECTED  EARLY  92 

• SUBSCALE  NOZZLE  IN  FABRICATION 

92  PLANS 

• SUBSCALE  NOZZLE  HOT  FIRE  TESTING 

• DESIGN  UPDATE  AND  CONTINUED  FAB  DEVELOPMENT 


LEAD:  STPT-PRATT  & WHITNEYAV.  PALM  BEACH  FL 


DESIGN  FEATURES 

• UNCOOLED  CHAMBER 

• OPERATING  PRESSURE  SAME  AS  MAIN  CHAMBER 

DEVELOPMENT  STATUS 

• WORKHORSE  GAS  GENERATOR  TESTED 

• DESIGN  OPTIONS  IDENTIFIED 

92  PLANS 

• IGNITER  DESIGN  • CONCEPTUAL  GAS  GENERATOR  DESIGN 


LEAD:  STPT-AEROJET/SACRAMENTO  CA 


GAS  GENERATOR/PREBURNER  COMPARISON 


STME 


SSME 


Discriminators 

SSME* 

STME  I 

Number  of  parts 

591 

70 

Number  of  weld  counts 

20 

5 

Weld  overlays 

2 

0 

Sheet  Metal  Liners 

2 

0 

Cost  (TFU) 

1.17 

243  K 

Discharge  Pressure 

5460  PSIA 

2205  PSIA 

Turbine  Inlet  Temperature 

1900  R 

1600  R 

Weight 

143  LBS 

111  LBS 

Fuel  Prebumer  Only 


SPACE  TRANSPORTATION  PROPULSION  TEAM  (STPT) 


Senior  Executive 
Council 


• “Board  of  Directors* 

• Delegates  Authority  to  team 

• Establishment  of  Policy 

• Provision  of  Overall  Guidance 

• Final  Resolution  for  Team  Conflicts 

• One  Senior  Executive  from  each  Company 


Team 

Program  Office 
(TPO) 


Aerojet 

Propulsion  Division 
ALS/STME  Program 


• Thrust  Chamber  Assembly  (L) 

• Control  System  (L) 

• Injector  (L) 

• Gas  Generator 

• Igniter 

• Engine  Systems  D&D 

• Fuel  Turbopump  (Turbine) 

• System  Components 


• Project  Director  is  Final  Authority 

• Overall  Program  Management 
Responsibility 

• Principal  Director  from  each 
Company 


Pratt  & Whitney 
Space  Propulsion 
& Systems 
ALS/STME  Program 


» Engine  Systems  D&D  (L) 

► Oxidizer  Turbopump  (L) 

> Nozzle  (L) 

» Engine  Systems  Assy  & Test 
* Control  Systems 
' System  Components 


Rocketdyne 

Division 

ALS/STME  Program 


Fuel  Turbopump  Assy  (L) 

Main  Combustion  Chamber  (L) 
Engine  System  D&D 
Engine  Systems  Test 
Control  Systems 
System  Components 


§<d®m©5  <&  Enng8nn®®]r2iing 


SPACE  TRANSPORTATION  MAIN  ENGINE 

PARTICIPANTS 


•AEROJET 

\ / T* 

PHILLIPS 
. LABORATORY- 
'S ROCKETDYNE 


"Ti»  LEWIS  RESEARCH  CENTER 

1 / 

SPACE  TRANSPORTATION 
'PROPULSION  TEAM 
HEADQUARTERS 
]?/—  s' 

^"MARSHALL  SPACE 
W FLIGHT  CENTER 


STENNIS 

SPACE 


‘PRATT  & WHITNEY 


PROGRAM  STATUS 


THRUST  LEVEL  INCREASED  TO  650K 

• REQUIRED  TO  SUPPORT  20K  VEHICLE 

• DESIGNS  BEING  REVISED  TO  REFLECT  INCREASE 

FULL  SCALE  DEVELOPMENT 

• REQUEST  FOR  PROPOSALS  IN  PREPARATION 

• SCHEDULE  FOR  RELEASE  - JUNE  1992 

• FIRST  ENGINE  SYSTEM  TEST  ■ MAY  1996 


OO 

N5 


arshall  Space  C enlcr 


i SdtenD®®  <&  IE 


SUMMARY  | 

• STME  IS  PROCESS  FOCUSED 

• PRODUCT  OF  INTEGRATED  SYSTEMS  PROCESS 

• INNOVATIVE  APPROACH  TO  PROJECT 

• DEFINITION  MATURE 

• BENEFITS 

- Robust  Propulsion 

- Legacy  of  New  Process 
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THE  IMPACT  OP  TIME  STEP  DEFINITION  ON 
CODE  CONVERGENCE  AND  ROBUSTNESS 

S . Venkateswaran,  J.  M.  Weiss  and  C.  L.  Merkle 
Propulsion  Engineering  Research  Center 
The  Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

A large  fraction  of  the  Navier-Stokes  codes  in  use  today  are 
based  on  the  so-called  'time-marching'  procedures,  wherein  the 
unsteady  fora  of  the  governing  equations  are  solved  in  time. 

For  compressible  flows,  these  methods  perform  very  well  in  the 
transonic  and  supersonic  flow  regimes  and  have  been  applied  to 
solving  a wide  variety  of  problems.  There  are,  however,  several 
disadvantages  associated  with  these  methods.  It  is  well  known  that 
at  low  Mach  numbers,  the  convergence  of  these  schemes  deteriorate 
dramatically.  The  reason  for  the  behaviour  is  the  wide  disparity 
in  the  eigenvalues  of  the  system  at  low  speeds.  Furthermore, 
highly  viscous  regions  of  the  flow  and  the  presence  of  strong 
source  terms  also  introduce  convergence  difficulties.  These  again 
are  due  to  the  very  disparate  time  scales  involved  in  these 
processes. 

Preconditioning  offers  a means  of  controlling  the  time-step 
size  for  a wide  variety  of  flow  situations.  Originally, 
preconditioning  methods  were  developed  as  a means  of 
circumventing  the  disparity  in  the  eigenvalues  at  low  Mach  numbers. 
Essentially,  this  involves  altering  the  time  derivative  of  the 
equations  of  motion  such  that  the  acoustic  speed  is  scaled  down  to 
the  level  of  the  fluid  velocity.  This  'inviscid'  preconditioning 
enables  Mach  number- independent  convergence  to  be  obtained. 

We  have  also  extended  the  preconditioning  approach  to  handling  very 
viscous  flows.  Here,  the  acoustic  speed  is  altered  such  that  the 
local  CFL  number  is  approximately  the  same  order  as  the  viscous 
time  step.  This  enables  excellent  convergence  rates  to  be 
maintained  over  a wide  range  of  Reynolds  numbers.  We  are  currently 
investigating  extending  this  approach  to  various  source  terms  of 
interest particularly  related  to  reacting  flowfields. 

We  have  implemented  preconditioning  for  multi-species  reacting 

flows  in  two  independent  codes an  implicit  (ADI)  code  developed 

in-house  and  the  RPLUS  code  (developed  at  NASA- Lewis  Research 
Center) . The  RPLUS  code  has  been  modified  to  work  on  a 4-stage 
Runge-Kutta  scheme.  The  performance  of  both  the  codes  have 
been  tested  and  show  that  preconditioning  can  improve  convergence 
by  a factor  of  two  to  a hundred  depending  on  the  problem. 

Our  efforts  are  currently  focussed  on  evaluating  the  effect  of 
chemical  sources  and  on  assessing  how  preconditioning  may  be 
applied  to  improve  convergence  and  robustness  in  the  calculation  of 
reacting  flows. 
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The  Impact  of  Time-Step  Definition  on 
Code  Convergence  and  Robustness 


S.  Venkateswaran 
Jonathan  M.  Weiss 
Charles  L.  Merkle 

Propulsion  Engineering  Research  Center 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Problems  Associated  with  Conventional  Codes 


• Convergence  is  poor  at  low  Mach  Numbers. 

e.g.,  the  combustion  zone  in  rocket  engines. 

• Convergence  is  poor  in  viscous  regions. 

e.g.,  boundary  layers,  recirculation  zones,  etc. 

• Large  source  terms  induce  instabilities. 

e.g.,  combustion,  turbulence,  axisymmetry,  etc. 

• High  aspect  ratio  grids  cause  poor  convergence. 

in  regions  where  strong  local  grid  stretching  is  used. 


Time- Step  Definition  in  Conventional  Codes 


Eigenvalues  of  Jacobian  A define  inviscid  time-step. 
— Eigenvalues  are  u-fc,  u-c,  u,  u,  u,  etc. 


CFLX  = 


XAt 

Ax 


• In  viscous  regions,  a viscous  time-step  is  defined. 


VNN  = 


vAt 

A? 


• Time  step  is  fixed  based  on  CFL  and  VNN  conditions. 


• Variable  time  stepping  is  normally  used. 


Problems  with  Conventional  Definition 


• At  low  Mach  numbers,  CFLu+c  > CFLU. 

— Acoustic  wave  speeds  dominate  over  particle  convection  speed. 

• In  viscous  regions,  VNN  > CFL. 

— Diffusion  time  scales  are  much  larger  than  wave  speeds  for 
acoustic  and  particle  convection. 

• For  high  aspect  ratio  cells,  CFLv+c  > CFLu+c. 

— Wave  speeds  are  much  higher  in  the  direction  normal  to  the 
flow  direction. 


• Source  terms  introduce  additional  time  scales. 


Preconditioning  the  Equations  of  Motion 


• Alter  the  time  derivative  by  multiplying  with  a preconditioning 
matrix  T. 

dQv  dE  dF  dEv  dFv 

dt  dx  + dy  + dx  + dy 

— Qv  = ( P ; «,  v,  T,  Yi,  F2,  • • • ) 


• Steady  state  solution  remains  unaltered. 


• Eigenvalues  of  T-1A  now  define  CFL  number. 


Choosing  the  Preconditioning  Matrix 


• Define  T so  that  the  acoustic  speeds  are  altered. 

• At  low  Mach  numbers,  keep  acoustic  speeds  of  the  same  order 
as  fluid  velocities. 

• In  viscous  regions,  alter  acoustic  speed  so  that  the  inviscid  time- 
step  is  of  the  same  order  as  the  viscous  time-step. 

• Scale  acoustic  speed  in  a similar  manner  to  account  for  large 
source  terms. 

• Extend  philosophy  to  high  aspect  ratio  grid  cells. 


Example  of  Preconditioning 


i 


Comparison  of  Convegence 


r ra ration  step 


Implementation  of  Preconditioning 


• Incorporate  preconditioning  in  implicit,  reacting  flow  code. 

— Euler  Implicit /ADI  Algorithm 

— Two  D/ Axisymmetric  Code 

— Multi- Component  Species  Transport 

— Multi-Step  Finite  Rate  Chemistry 

• Incorporate  preconditioning  into  RPLUS  code. 

— Developed  at  NASA-Lewis  Research  Center. 

— Multi-Stage  Explicit  Runge-Kutta  (RPLUS/RK) 


Preconditioning  in  the  RPLUS  Code 


• Modified  Runge-Kutta  Stage: 


A Qv  = -afcA t r 1 


,dE  dF  dEv 

dx  dy  dx 


dFv) 
dy  > 


• Time-Step  definition  altered  using  the  new  eigenvalues  of  the 
system. 


Density  Residual 


Convergence  of  RPLUS/RK 


Density  Residual 


Density  Residual 


Convergence  of  RPLUS/RK 


Reacting  Shear  Layer 


i 1 1 1 1 1 1 1 i I 

1000.  / 2000.  3000.  4000.  5000. 

Iteration 


26  Apr  92  20:50: 14 


Convergence  of  Implicit  ADI 


27  Apr  92  10:25:26 
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Definition  of  Time- Step 


A t = 


CFLAx 

A 


• Based  on  the  maximum  eigenvalue: 


A = Max  (w  + c,  v + c) 

• Based  on  an  average  eigenvalue: 

A = +(u  + 

A = l/2(  u + c + v + c ) 

• Based  on  the  eigenvalue  in  the  direction  of  flow: 

A = u + c 


i 


Axisymmetric  Convergence 


Maximum  Eigenvalue 


X- Eigenvalue 


? -50  V'-.- 


-9.0 

*V 


200.  400.  600.  800.  1000.  12(H). 


j \ 71X31  CFL=3 

*15.0  -| 1 1 1 1 j 1 1 1 1 j 1 1 1 

/,(X’  0.  200.  400.  600.  800.  1000.  1200.  1< 


~TtEA-/'TI0'<'  SrE-?> 


Convergence  Based  on  X- Eigenvalue 


Axisymmetric 


-J  r£ KATio/Y  SreF 
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Time  Step  Definition 


• Important  to  use  proper  eigenvalue  in  CFL. 

— u-|-c 

— Max(u+c,v4-c) 

— + c)2  + (v  + c)2 

• Preliminary  results  for  H-grids  show: 

— best  choice  is  u+c 

— control  convergence  with  grid  refinement 

— control  convergence  in  near  wall  region 

• Additional  work  needed  to  generalize 


Conclusions 


• The  definition  of  time-step  has  a profound  impact  on  the  per- 
formance of  time-marching  codes. 

• Preconditioning  is  a powerful  method  of  controlling  the  time- 
step. 

— Low  Mach  number  preconditioning  or  characteristic  time- 
stepping has  been  used  widely. 

— Preconditioning  lias  been  successfully  extended  to  viscous 
dominated  flows. 

— Similar  extensions  are  currently  being  investigated  for  com- 
bustion and  other  sources  of  interest. 

• Time-step  should  be  defined  based  on  the  eigenvalue  in  the  di- 
rection of  flow. 

— Important  when  the  grid  aspect  ratio  is  very  high. 
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Development  of  CFD  Code  Evaluation  Criteria  and  a Procedure  for 
Assessing  Predictive  Capability  and  Performance 

S.J.  Lin,  D.C.  Chan,  M.M.  Sindir,  and  S.L.  Barson 
Rockwell  International,  Rocketdyne  Division 
Canoga  Park,  California 


Careful  validation  of  Computational  Fluid  Dynamic  codes  is  essential  if  they  are  to 
be  used  as  engineering  design  tools.  Validation  must  be  carried  out  in  a 
systematic  manner  to  ensure  that  all  code  aspects  as  they  apply  to  the  application 
of  interest  are  understood  and,  to  the  greatest  extent  possible,  quantified. 

A study  is  being  conducted  in  which  a general  code  validation  procedure  is  defined 
and  demonstrated.  A four  phase  validation  procedure  is  defined  in  which  a series 
of  validation  test  cases  are  computed  and  compared  with  available  analytical 
solutions  and  test  data.  The  procedure  is  demonstrated  using  the  REACT  CFD 
code  to  compute  validation  cases  for  each  of  the  four  phases.  For  phase  4,  the 
application  of  interest,  the  SSME  high  pressure  fuel  turbopump  impeller  flowfield  is 
computed. 


preceding 


"•’AGE 


blawk  not  filmed 
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DEVELOPMENT  OF  CFD  CODE  EVALUATION  CRITERIA 
AND  PROCEDURE  FOR  ASSESSING  PREDICTIVE 
CAPABILITY  AND  PERFORMANCE 


S.J.  Lin,,  D.C.  Chan,  M.M.  Sindir,  and  S.L.  Barson 
Rockwell  International,  Rocketdyne  Division 


Workshop  for  Computational  Fluid  Dynamic 
Applications  in  Rocket  Propulsion 


April  28-30, 1992 

NASA  Marshall  Space  Flight  Center 


Rockwell  International 

Rocketdyne  Division 


CFO  92-030-00 1/02/SLB  J 


DEVELOPMENT  OF  CODE  EVALUATION 
CRITERIA  AND  A PROCEDURE 


• TASK  OBJECTIVES 

• PROVIDE  CODE  EVALUATION  CRITERIA,  CLASSIFICATION  SCHEME, 
NUMERICAL  ERROR  ASSESSMENT  TECHNIQUES,  AND  A PROCEDURE 
FOR  COMPREHENSIVE  CODE  EVALUATION  AND  CERTIFICATION 

• ENSURE  INTEGRITY,  ACCURACY,  AND  APPLICABILITY  OF  CFD  CODES 

• PROVIDE  PROCEDURES  AND  GUIDELINES  FOR  CFD  SOFTWARE 
QUALITY  CONTROL 

• DEMONSTRATE  CODE  EVALUATION  PROCEDURE  USING  2-D  AND  3-D 
BENCHMARK  EXPERIMENTS. 


• PRESENTATION  FOCUS 

• CODE  VALIDATION  PROCEDURE 

• DEMONSTRATION  OF  PROCEDURE 


Rockwell  International 

Rocketdyne  Division 


CFD  92 -030-002/D  2/SLB 
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THOUGHTS  ON  VALIDATION 

• GENERAL  VALIDATION  PROCEDURE  FOR  ALL  APPLICATIONS 
IS  POSSIBLE 

• NO  GENERAL  AND  ABSOLUTE  VALIDATION  POSSIBLE  FOR 
ALL  CASES 

• QUANTITATIVE  VALIDATION  ONLY  MEANINGFUL  WITHIN  A 
LIMITED  CLASS  OF  APPLICATIONS 

• LEVEL  OF  VALIDATION  DEPENDS  ON  FINAL  APPLICATION 

• VALIDATION  PROCESS  MUST  BE  REALISTICALLY  ACHIEVABLE 


Rockwell  International 

Rocketdyne  Division 


CFD  02-030-003/D  2/S  LB 
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FOUR  PHASE  CODE  VALIDATION  PROCEDURE  DEFINED 


PHASE  1 
UNIT  PROBLEM 


PHASE  2 

BENCHMARK  CASES 


PHASE  3 

SIMPLIFIED  PARTIAL 
FLOWPATH 


PHASE  4 

ACTUAL  HARDWARE 


• SINGLE  FLOW 

• MORE  THAN  ONE 

• MULTIPLE  RELEVANT 

• COMPLETE  FLOW 

FEATURE 

FLOW  FEATURE 

FLOW  FEATURES 

PHYSICS 

• ANALYTIC  SOLUTION 

• SIMPLE  FLOW 

• ACTUAL  FLOW 

• HARDWARE  TEST 

OR  HIGH  FIDELITY 
COMPUTATIONAL 
SOLUTION  (DNS) 
AVAILABLE 

PHYSICS 

• BENCHMARK 
EXPERIMENT 
DATA 

PHYSICS 

• HIGH  QUALITY 
TEST  DATA 

DATA 

INCREASING  GEOMETFlIC  AND  FLOW  COMPLEXITY 


RUN  UNIT  PROBLEMS 
VERIFY  INTEGRITY 
ASSESS  ACCURACY, 
CONVERGENCE, 

AND  FUNCTIONALITY 


• RUN  BENCHMARK 

• RUN  SIMPLIFIED 

• RUN  ACTUAL 

CASES 

PARTIAL  FLOWPATH 

CONFIGURATION 

• ASSESS  PHYSICAL 

• ASSESS  AGREEMENT 

• COMPARE  WITH 

MODELS 

WITH  DATA 

TEST  DATA 

• ESTABLISH  GRID 

• ESTABLISH  GRID 

DISTRIBUTION 

DISTRIBUTION 

REQUIREMENTS 

REQUIREMENTS 

DECREASING  DATA  AVAILABILITY  AND  ACCURACY 


Rockwell  International 

Rocketdyne  Division 


CFD  82-030-004/D2/SL 


00  0(0)000X000 


DECREASING  NUMBER  OF  CASES  REQUIRED 
FOR  LATTER  VALIDATION  PHASES 


UNIT 

PROBLEMS 


Rockwell  International 

Rocketdyne  Division 


CFO  B2-0M-005/D2/SLB 


PROBLEM  OF  INTEREST  SUCCESSIVELY 
DECOMPOSED  INTO  LESS  COMPLEX  CASES 

EXAMPLE:  SSME  HPFTP  IMPELLER 


PHASE  1 

UNIT 

PROBLEMS 

• flat  plate 

• straight  duct 

• diffuser 

• sudden 
contraction  (lam.) 

• backward  facing 
step  (lam.) 

• driven  cavity 

• rotating 
concentric 
cylinders  (Taylor- 
Couette  flow) 


PHASE  2 

BENCHMARK 

CASES 

• square  duct  with 
90°  bend 

• S-shaped  duct 

• backward  facing 
step  (turb.) 

• orifice  flow  (turb.) 

• flow  around 
confined  bluff 
bodies 

• 2-D  turbine  cascade 

• rotating  disk 


PHASE  3 

SIMPLIFIED 

FLOWPATHS 


• 3-D  turbine  blade 
cascade 

• rotating  curved 
duct 


PHASE  4 

ACTUAL 

HARDWARE 


• SSME  HPFTP 
impeller  (2  sets 
partial  blades) 


Rockwell  International 

Rocketdyne  Division 


CFD  92-030-006/D2/SLB 


REACT*  CODE  DESCRIPTION 


• CO-DEVELOPED  BY  ROCKETDYNE/UNIVERSITY  OF  LONDON 

• 2-D/3-D,  STEADY  STATE  OR  TRANSIENT,  FULL  NAVIER-STOKES 

• MULTI-ZONE  FINITE-VOLUME  IN  GENERALIZED  COORDINATES 

• PRESSURE-VELOCITY  COUPLING  THROUGH  "SIMPLE"  AND  "PISO" 

• STONE'S  STRONGLY  IMPLICIT  AND  CONJUGATE  GRADIENT  SOLVERS 

• VARIOUS  2-EQUATION  TURBULENCE  MODELS 

• CONJUGATE  FLUID-SOLID  HEAT  TRANSFER  CAPABILITY 

• MULTI-SPECIES  CAPABILITY 

• PRIMARY  USE  FOR  TURBOMACHINERY  APPLICATIONS 


Residual 


UNIT  PROBLEMS  COMPUTED  FOR  PHASE  1 

DUCT  FLOW 


Verify  Single  Zone  and 
Multizone  Convergence 


Establish  Multizone  Grid  Matching  Requirements 


Rockwell  International 

Rocketdyne  Division 


CFD  92-030-008/D  2/S  LB 
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BENCHMARK  CASES  COMPUTED  FOR  PHASE  2 

CURVED  DUCT  FLOW 


90°  Bend  With  Square  Cross  Section 


Secondary  Flow  Features  Captured 


Refined  Grid  Yields  Improved 

Agreement  With  Experimental  Data 


Experimental  Data  * 
88x22x12  Grid 


88x42x22  Grid 

88x82x42  Grid 


04  04  00  Ut  04  00  0* 

radial  veiocity 


Pressure  Contours 


Rockwell  International 

Rocketdyne  Division 


Velocity  Profiles  At  77.5°  Location 

A.M.K.  Taylor,  J.H.  Whitelaw,  And  M.  Yianneskls,  “Measurements  of 
Laminar  and  Turbulent  Flow  In  a Curved  Duct  with  Thin  Inlet  Boundary 
Layers,"  NASA  Contractor  Report  3367 

CFD  02-O3O-OO5VD2/SLB  J 
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PARTIAL  FLOWPATH  COMPUTED  FOR  PHASE  3 

HPFTP  FIRST  STAGE  TURBINE  CASCADE 


SSME  HPFTP  First  Stage  Stator 
Multlzone  O-H  Grid 


• STUDY  EFFECTS  OF  SINGLE  ZONE  H GRID  AND  MULTIZONE  O-H  GRID 

• PERFORM  COARSE  AND  FINE  GRID  COMPUTATIONS  ON  EACH  AND 
COMPARE  WITH  TEST  DATA 


Static  Pressure  At  Stator  Mid-section 
Fine  Grid  Solution 


in 


Rockwell  International 

Rocketdyne  Division 


CFD  «-OXM)1(yp2/SLBy 
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COMPUTATIONS  ON  ACTUAL  HARDWARE 
CONFIGURATION  IN  PROGRESS  FOR  PHASE  4 


PROBLEM  FEATURES 

• HIGHLY  THREE-DIMENSIONAL 

• TWO  SETS  OF  PARTIAL  BLADES 

• STRONG  CURVATURE 

• HIGH  ROTATIONAL  SPEEDS 

• TURBULENT  FLOW 

MODELING  APPROACH 

• 3-D  MULTIZONE  MODEL 

• k-e  TURBULENCE  MODEL 

STATUS  AND  PLANS 

• COARSE  GRID  SOLUTION  COMPLETED 
WITH  ASSUMED  INLET  CONDITION 

• TWO  FINE  GRID  COMPUTATIONS 
PLANNED 

• ASSUMED  INLET  CONDITION 

• INLET  CONDITION  FROM  TEST  DATA 

• DATA  ACQUISITION  IN  PROGRESS 

• INLET  FLOW  DATA  NOW  AVAILABLE 

• OUTLET  DATA  AVAILABLE  SOON 


Rockwell  International 

Rocketdyne  Division 


CFO  92-030-01 1/D2/SLB  \ 
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SUMMARY 


• GENERAL  VALIDATION  PROCEDURE  DEFINED  FOR 
ALL  APPLICATIONS 

• FOUR  PHASE  PROCEDURE  OUTLINED 

• QUANTITATIVE  VALIDATION  ONLY  MEANINGFUL 
WITHIN  LIMITED  CLASS  OF  APPLICATIONS 

• CRITERIA  BEING  DEFINED 

I 

i 

• REALISTIC  VALIDATION  PROCEDURE  DEMONSTRATED 
ON  ACTUAL  HARDWARE! 

• SSM$HPFTP  IMPELLER 

• DATA  ACQUISITION  IN  PROGRESS 

• FINAlL  COMPUTATIONS  AND  DATA  COMPARISONS 
TO  FOLLOW 


Rockwell  International 

Rocketdyne  Division 


CFO  92-03Q-01 2/0  2/S  LB 
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COMPARISON  BETWEEN  THE  PISO  ALGORITHM  AND 
PRECONDITIONING  METHODS  FOR  COMPRESSIBLE  FLOW 


Charles  L.  Merkle,  Philip  E.  O.  Buelow  and  S.  Venkateswaran 
Propulsion  Engineering  Research  Center 
The  Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802. 


Two  widely  used  family  of  algorithms,  pressure-based  and  density-based 
methods,  have  been  developed  for  CFD  problems  over  the  years.  Pressure-based 
methods  (such  as  SIMPLE  and  PISO)  use  a Poisson-like  equation  for  updating 
pressure  instead  of  the  continuity  equation,  while  density-based  methods  use  the 
continuity  equation  to  update  density  (an  equation  of  state  is  used  to  provide  density 
in  pressure  based  schemes  and  pressure  in  density  based  schemes).  Pressure-based 
methods  were  developed  originally  for  incompressible  flows  at  low  Reynolds  numbers 
and  were  then  extended  to  high  Reynolds  numbers  and  compressible  applications. 

On  the  other  hand,  density  based  methods  were  originally  developed  for  transonic 
flows  and  have  been  extended  down  to  low  Mach  numbers  through  the  use  of 
preconditioning  techniques.  Both  methods  have  enjoyed  considerable  success  in 
solving  complex  flowfields,  though  the  relative  effectiveness  of  the  schemes  has  long 
been  argued.  Generally,  pressure-based  methods  are  more  robust  while  density-based 
schemes  are  more  temperamental  but  provide  more  accurate  solutions. 

In  the  present  paper,  we  compare  these  two  very  different  approaches  to  solving 
the  Navier-Stokes  equations  in  order  to  gain  an  understanding  of  their  similarities 
and  differences.  Specifically,  we  consider  the  PISO  scheme  as  a representative 
pressure-based  method  and  contrast  it  with  a recently  developed  preconditioning 
scheme  To  facilitate  the  comparison,  we  write  both  schemes  in  a vector  formulation. 
Our  findings  indicate  that  the  PISO  scheme  is  very  closely  related  to  the  philosophy 
of  the  preconditioning  scheme.  In  particular,  preconditioning  causes  the  density- 
based  scheme  to  appear  pressure-based  at  low  speeds  but  to  remain  density-based  at 
high  speeds.  Furthermore,  both  schemes  alter  the  sonic  speed  so  that  the  equations 
stay  well  conditioned  in  the  limit  of  low  Mach  numbers. 

We  also  compare  the  relative  performance  of  the  PISO  algorithm  with  an  Euler 
implicit  algorithm  that  is  employed  to  solve  the  preconditioned  equations  by  means 
of  a vector  stability  analysis.  The  results  of  the  stability  analysis  indicate  that 
the  PISO  algorithm,  which  is  a multi-step  (one  predictor  step  followed  by  severed 
corrector  steps),  uncoupled  (i.e.,  sequential)  solution  procedure,  is  conditionally 
stable.  Good  convergence  is  promised  at  low  CFL  numbers,  while  at  high  CFL 
numbers,  both  low  wave  number  and  high  wave  number  instabilities  are  present. 

The  high  wave  number  instability  appears  to  be  ‘compressible’  in  origin,  arising 
from  the  treatment  of  the  equation  of  state.  The  low  wave  number  instability  is 
‘incompressible’  in  origin  since  it  is  present  when  the  incompressible  limit  of  the 
equations  are  examined.  An  important  finding,  in  this  regard,  is  that  the  overall 
scheme  may  be  unstable  even  when  the  individual  predictor  and  corrector  stages 
are  themselves  stable.  In  contrast,  the  Euler  implicit  algorithm  shows  unconditional 
stability.  It  should  be  noted,  however,  that  multi- dimenisional  solution  of  the 
equations  demands  the  use  of  approximate  factorization  which  limits  CFL  numbers 
to  about  10.  Thus,  the  two  algorithms  still  remain  quite  competetive  in  solving 
practiced  flow  problems. 


'IK  MOT  FILMED 


PRECEDING  PAGE  BLAr 


123 


124 


The  Relationship  Between 
Pressure-  and  Density-Based  Algorithms 


Charles  L.  Merkle,  Sankaran  Venkateswaran 
and  Philip  E.  O.  Buelow 
The  Pennsylvania  State  University 
Propulsion  Engineering  Research  Center 
Department  of  Mechanical  Engineering 


Presented  at 


Computational  Fluid  Dynamics  Workshop 
Marshall  Space  Flight  Center 
April  28-30,1992 


125 


Introduction 


Compare  Pressure-Based  and  Density-Based  Methods 
Pressure-Based 

- SIMPLE,  MAC.PISO  etc. 

- Replace  Continuity  by  Poisson  Equation 

- Solve  by  Sequential  Procedure 
Density-Based 

- ADI,  LU,  Lax-Wendroff 

- Solve  Continuity  Directly 

- Solve  by  Simultaneous,  Coupled  Procedure 
Express  in  Common  Vector  Form  for  Comparison 
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Development  of  Pressure  Poisson  Relation 


Use  Continuity  with  Source: 

v • y = d 


Discretize  Momentum: 

un+I  - iT  (du2  c)uv  3p 
+ — + — +— 

At  ^ dx  dy  dx 

vn+1_vn  3UV  0V2  0p 

+ — -f — +— 

At  dx  dy  dy 


Solution  of  Poisson  Equation 

Solve  by  Point  Jacobi  with  OverRelaxation: 


- Express  as  Equivalent  Time  Marching 


4At  9pn+' 
toAx2  at 


= V2p*  +o"  fv*  v) 

At'  1 


where: 


n-1 


By  Comparison  of  Equations: 

4At  3pn+'  1 / %n+l  I 

7—7 — = — (V  • V)  =- 

WAX2  at  At'  ' At 

Hence  the  Equivatent  Equation  Can  Be  Written 


D 


4At2  I 

COAX2  At 


ap°"  apn 
at  at 


= V2p*  +cr* 


Poisson  Method  is  Hyperbolic  if  Is  Retained 
Poisson  Method  is  Parabolic  i f —;Dn  |s  Dropped 
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Poisson  Equation  in  Compressible  PISO  Method 


Continuity  Equation: 
p -p 


At 


+ 


3pu  3pv 
^ dx  dy  J 


= 0 


Or,  Using  Perfect  Gas  Relation 

P -P 


\n  f , "N 


RTnAt 


+ 


dpu  ^pv 


V 


dx  dy  J 


= 0 


Momentum  Equations: 


n,  ,n 


p U -p'V 
At 


■+ 


rd  u2 


+ 


9uv 


\ 


\ 


dx  dy 


J 


uv  ww 


p v — pnvn 
At 


+ 


fauv  av2) 

H 

ax  dy  J 


+ * .0 

dX 

dp 


*-* 


+ 


ay 


= 0 


*-* 

P 


=p**/rt* 


Combine  from  Divergence  of  Momentum 


Characteristic  Speed  in  PISO  Poisson  Equation 

• Combined  Continuity  and  Momentum  Equations 

P**~Pn  +_Lfej+jPXl  - a*  - V2p*“  =0 

RTnAt2  Atl^ax  ay  J 

• Replace  Divergence  with  Density  Derivative 

p~_2p"+pn-'=  o„ 

RTnAt2 

• PISO  Poisson  Equation  is  Hyperbolic 

- Characteristic  Speeds  Are  the  Acoustic  Speeds 

• LowMach  Number  Convergence  Requires: 

- Multiple  Sweeps  of  Continuity  Equation 

- Re-scaling  of  Time  Derivative-Preconditioning 


Equations  of  Motion 


3Q  3E  3F 

at+ax  ay 


where 

Q=(p,pu,pv;e)T 

E=(pu,pu2+p,puv;eu+pu)T  (2) 


h =(pV)pUV,pV2+P,6V+pV) 


T 


(Qv)= 


— RXv^-Qv 

ax  xxax  v 


aD  a 

+axRxyay 


Q 


V 


R 


ay  yx3x 


Q 


131 


Formulation  of  PISO  Algorithm 


Split  Flux  Vectors: 

EL=(pu,p,0,o£ 


E=El+En 

EN=(0,pu2,puv,(e+p)u) 

C.U. 


Use  Predictor-Corrector  Procedure: 


d A n 9n  n 
8xAn  +3vBn 


c+Vn  r+At  — An  +— Bn  -Lv  fAQ  ——At lm 


3E  3F  . ,n  , 

z -+ r LV(QV) 

ax  ay 


a . * d 


r+  +^bl  acT  +U^Am  +^Bn  -L 


a . * a „ * 


a * * P)  * * | ( P)  * * p±  * * - 

— A|_  Qv  +— Bl  Qv  j +le  Qv  Qy  -LV(QV 


Stability  Analysis 

Represent  Disturbance  Growth  by  Amplification  Matrix 

Qn+'=GQn 

Result  Provides  Four  Amplification  Factors 
Plot  Maximum  of  These  Four 


For  PISO  Scheme, 


y-wavenumber 


CFL=1,  M=0.3 


CFL=1,  M=0.03 


y-wavenumber 


MACH  NUMBER  EFFECT  FOR  PISO  SCHEME 


CFL=1,  M=0.3  CFL=1,  M=0.03 


x-wavcnumber 


x-wavenumber 


y-wavenumber 


TIME  DER1VA' 


COUPLED  SYSTEM 


0.0  TT 


x-wavenumber 


y-wavenumber 


EFFECT  ON  PISO  SCHEME 


0.0 


x-wavenumber 


log(residual) 


log(residual) 


Comparis<  n of  Convergence 


inflo  v/ outflow  B.C.’s 


Euler  Implicit 


PISO 


LO 

^1 


steps 


I I 

5600. 


Y wavenumber 


PISO:  Gauss  Seidel  (10  sweeps) 


CFL  = 1 


X wavenumber 
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Summary 


♦ Differences  in  Pressure-  and  Density-Based  Methods. 

Upwind  Direction 

Choice  of  Solution  Variables 

Coupled  vs.  Uncoupled  Equations 

. Pressure-Based  Methods  for  Incompressible  Flows  Are  Hyper- 
bolic 

Not  Parabolic 

* Pressure-Based  Methods;  for  Compressible  Flow  are  Hyperbolic 
Characteristics  of  Poisson  Equation  are  Stiff 

— Time-step  Control  is  Needed  for  Convergence 

— Can  be  Offset  by  Multiple  Sweeps  of  the  Poisson  Equation 


Summary  (Contd.) 


• Vector  Form  of  Pressure-Based  Method  Facilitates  Comparison 

• PISO  Vector  Stability  Analysis  Indicates: 

— Conditionally  Stable 

— Low  Wave  Number  Instability  (Incompressible) 

— High  Wave  Number  Instability  (Compressible) 

• Code  Convergence  Verifies  Stability  Predictions 

• Approximate  Factorization  of  Poisson  Equation 

— Low  Mach  Number  Stiffness 

— Mitigate  by  Scaling  Time  Step 

— Circumvented  by  Gauss-Seidel 
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A Comparison  of  Artificial  Compressibility  and  Fractional  Step 
Methods  for  Incompressible  Flow  Computations 

Daniel  C.  Chan 

Department  of  Aerospace  Engineering,  University  of  Southern 

California, 

Los  Angeles,  California 
and 

Rocketdyne  Division,  Rockwell  International  Corporation 
Armen  Darian  and  Munir  Sindir 

Rocketdyne  Division,  Rockwell  International  Corporation 

Canoga  Park,  California 

We  have  applied  and  compared  the  efficiency  and  accuracy  of 
two  commonly  used  numerical  methods  for  the  solution  of  Navier- 
Stokes  equations.  The  artificial  compressibility  method,  postulated  by 
Chorin,  augments  the  continuity  equation  with  a transient  pressure 
term  and  allows  one  to  solve  the  modified  equations  as  a coupled 
system.  Due  to  its  implicit  nature,  one  can  have  the  luxury  of  taking  a 
large  temporal  integration  step  in  the  expenses  of  higher  memory 
requirement  and  larger  operation  counts  per  step.  Meanwhile,  the 
fractional  step  method,  developed  independently  by  Chorin  and 
Temam,  splits  the  Navier-Stokes  equations  into  a sequence  of 
differential  operators  and  integrates  them  in  multiple  steps.  The 
memory  requirement  and  operation  count  per  time  step  are  low, 
however,  the  restriction  on  the  size  of  time  marching  step  is  more 
severe. 

To  explore  the  strength  and  weakness  of  these  two  methods,  we 
used  them  for  the  computation  of  a two-dimensional  driven  cavity  flow 
with  Reynolds  number  of  100  and  1000,  respectively.  Three  grid 
sizes,  41  x41 , 81  x81  and  161x161  were  used.  The  computations  were 
considered  converged  after  the  L2-norm  of  the  change  of  the 
dependent  variables  in  two  consecutive  time  steps  has  fallen  below 
10‘5.  Same  programming  style  is  applied  to  the  development  of  these 
codes.  All  computations  were  performed  on  the  NASA-Marshall 
Convex  C240  computer  with  double  precision  arithmetic. 

In  summary,  we  find  that  the  artificial  compressibility  method 

requires  twice  as  much  memory  per  grid  points  and  is  less  efficient  for 
grid  resolution  below  81x81.  Fractional  step  method,  on  the  other 
hand,  is  more  efficient  in  both  memory  requirement  and 
computational  speed  for  coarse  grid  computations,  however,  due  to 
its  explicit  nature,  its  convergence  rate  deteriorates  dramatically  for 
fine  grid  computations 
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AGENDA 


• MOTIVATION 

• APPROACH 

• TEST  CASE  DESCRIPTION 

• RESULTS 
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MOTIVATION 


• EXPLORE  AN  ALTERNATIVE  NAVIER-STOKES  SOLVER 
FOR  ENGINEERING  APPLICATIONS 

• INTERNAL,  INCOMPRESSIBLE  FLOWS 

• IDENTIFY  THE  STRENGTH  AND  WEAKNESS 

• LEVEL  OF  EXPERTISE  REQUIRED  ' 

• ACCURACY 

• SPEED 


Rockwell  International 
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APPROACH 


• COUPLED  METHOD 

• EXTENSION  OF  COMPRESSIBLE  FLOW  FORMULATION 

• HIGHLY  IMPLICIT 

• ELEGANT  MATHEMATICAL  FORMULATION 

• ITERATIVE  METHOD 

• LESS  IMPLICIT 

• LONGER  TRACK  RECORD 

• LESS  COMPLICATED  TO  FORMULATE 
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GOVERNING  EQUATIONS  FOR 
ARTIFICAL  COMPRESSIBILITY  METHOD 


1 I [D]V24  = 0 
dt  dx  dy  Re1  J 


p 1 

‘ «2  ' 

9 

M 

F = 

+ p 

V 

UV 
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_ 

V 

* 

0 0 O' 

G = 

UV 

v^  + p 

[°1= 

0 1 0 
0 0 1 
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NUMERICS 


• SECOND  ORDER  CENTRAL  DIFFERENCING  FOR  ALL  SPATIAL 
DERIVITIVES 


• IMPLICIT  TEMPORAL  INTEGRATION  IS  NEEDED  TO  OVERCOME 
STIFFNESS 


• APPROXIMATE  FACTORIZATION  IN  DELTA  FORM  PLUS  SECOND 
ORDER  IMPLICIT  AND  FOURTH  ORDER  EXPLICIT  DAMPING 
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BOUNDARY  CONDITIONS 


VELOCITY 

• NO-SLIP  ALONG  SOLID  WALLS 


• PRESSURE 

• NEUMANN  CONDITION  DERIVED  FROM  MOMENTUM 
EQUATIONS 

• EXPLICIT  IMPLEMENTATION 
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FRACTIONAL  STEP  METHOD 


• INTEGRATE  DIFFERENTIAL  OPERATORS  IN  A SEQUENCE  OF 
STEPS 

u.  -un 

_! — J-=L(u.)-N(u.) 

At  1 1 


un+l-un 

— = — — --  where 

At  8x. 

i 


88 


88  _ 

8x. 


8P  I ( ^ 
— -L(u.) 

8x.  1 


CONTINUITY  REQUIRES 


8x.8x. 

i i 


At  8x. 
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FRACTIONAL  STEP  METHOD 
(CONT'D) 


• VARIOUS  MULTI-STEP  METHODS  CAN  BE  INCORPORATED  TO 
ACHIEVE  DESIRED  TEMPORAL  ACCURACY  AND  STABILITY 

• FULLY  IMPLICIT,  LEAP-FROG,  CRANK-NICOLSON, 
RUNGE-KUTTA 

• BOUNDARY  CONDITION 

• NO-SLIP  ALONG  WALL 

. GREEN'S  THEOREM  REQUIRES: 

j!^ds=funds=$^  = 0 
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A COMPARISON  OF  THE  TWO  METHODS 


TIME  INTEGRATION 


AFITIFICIAL 

COMPRESSIBILITY 

FOII Y TMPUC1T 

(NEWTON'S;  LINEARIZATION) 


FRACTIONAL  STEP 


DIFFUSION:  IMPLICIT 
PRESSURE:  IMPLICIT 


SPATIAL  DISCRETIZATION  CENTRAL  DIFFERENCING 


ARTIFICIAL  DAMPING 


4Th  Order  explicit; 

2ND  ORDEiR  IMPLICIT  FOR 
ALL  EQUATIONS 

explicit 

2IT70RDER 

30WORDS/POINT 

a1,  At,  e.,  e 
i e 


UPWIND 

DIFFUSION:  CENTRAL 
PRESSURE:  CENTRAL 

FINITE  VOLUME 

NON-STAGGERED 
NON-LINEAR,  EXPLICIT  4TH 
AND  2ND  ORDER  FOR 
PRESSURE  STEP 

EXPLICIT  ON  V 

IMPLICIT  ON  4> 

2ND  ORDER 

~20  WORDS/POINT 

At 
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COMPUTATIONAL  MODEL  FOR 
THE  FRACTIONAL  STEP  METHOD 
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COMPUTATIONAL.  MODEL  FOR  THE 
ARTIFICIAL  COMPRESSIBILTY  METHOD 


On 

VD 


11=1  dP  1 


v=0  ay  Re  dy2 


(0,0) 
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METHOD  OF  COMPARISON 


• CONSISTENT  PROGRAMMING  STYLE 

• ONE  DEVELOPER 

• NO  EXPLICIT  VECTORIZATION 

• USE  ONE  COMPUTER  WITH  SAME  COMPILATION  OPTION 

• CONVEX  C-240 

• fc  -02  -pd8 

• MAINTAIN  CONSTANT  AT  AT  EACH  GRID  POINT 

• 'OPTIMIZE'  INPUT  PARAMETERS  WITH  41x41  GRID  RESOLUTION 


v 

• ei  AND  Ce  FIXED  AT  0.2  AND  0.1 , RESPECTIVELY 


PARAMETERS  TO  COMPARE 


• RANGE  OF  OPERATION 

• CPU  TIME  REQUIREMENT  AS  A FUNCTION  OF 
GRID  POINTS  USED 

• REYNOLDS  NUMBER  DEPENDENCY 
. ACCURACY 
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L2  NORM 


CONVERGENCE  HISTORY  FOR 
41  x 41  GRID,  Re=1 00 


0=4.0,  AT=0.01  0= 2.0,  AT=0.01 


NUMBER  OF  TIME  STEPS  NUMBER  OF  TIME  STEPS 
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CONTINUITY  DEFECT 


L2  NORM  OF  AQ 


CONVERGENCE  HISTORY  FOR 
41  x 41  GRID,  Re=100 


0=2.0,  AT=0.025 


0=2.0,  AT=0.05 


NUMBER  OF  TIME  STEPS 


NUMBER  OF  TIME  STEPS 
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CONVERGENCE  HISTORY  FOR 
41  x 41  GRID,  Re=100 


jS=1.0,  AT=0.01 


fi= 1.0,  AT=0.025 


J io~5-- 


-0.01  8 


# * 

• » 

• • 

« * 

i • 

• • 

i # 

i • 

• « 

> 

i » 

i < 

1 < 

a • 

i ft 

' • ! 
» • 

9 9 i 

i * 

ji 

# ft 

i • 

• • 

/ 

• 

— 

a ft 

ft  ft 

ft  ft 

» > 

i 

200  400  600  800  1000  1200 

NUMBER  OF  TIME  STEPS 


200  400  600  800  1000  1200 

NUMBER  OF  TIME  STEPS 


Rockwell  International 

Rocketdyne  Division 


CFD-92-030-01 5/D2/AD 


CONTINUITY  DEFECT 


L2  NORM  OF  AO 


CONVERGENCE  HISTORY  FOR 
41  x 41  GRID,  Re=100 


0=1.0,  AT=0.05 


NUMBER  OF  TIME  STEPS 


0=1.0,  AT=0.1 
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CONTINUITY  DEFECT 


L2  NORM  OF  AQ 


CONVERGENCE  HISTORY  FOR 
81  x81  GRID,  Re=1000 


A 


AC  METHOD 
AT =0.05 


ON 

ON 


NUMBER  OF  TIME  STEPS 
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FS  METHOD 
AT=0.004 


NUMBER  OF  TIME  STEPS 
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CONTINUITY  DEFECT 


CONVEX  C240  CPU  TIME(SECOND) 


EFFECT  OF  TIME  STEP  ON  CONVERGENCE  RATE 

Re-100, 41X41  GRID 
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MAXIMUM  CFL  NUMBER 


ON 

VO 


DISTANCE  ALONG  VERTICAL  CENTERLINE 


VELOCITY  PROFILES  ALONG  THE  GEOMETRIC 
CENTERLINES  OF  A DRIVEN  CAVITY  FLOW 

Re=1 000 
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SUMMARY  OF  RESULTS 


. APPROXIMATE  FACTORIZATION/ARTIFICIAL 
COMPRESSIBILITY  METHOD 

• MEMORY  INTENSIVE 

• SENSITIVE  TO  REYNOLDS  NUMBER 

• REQUIRES  MORE  USER  INTERACTIONS 

• FRACTIONAL  STEP  METHOD 

• MEMORY  EFFICIENCY 

• SENSITIVE  TO  TIME  STEP  USED 

• COMPUTING  INTENSIVE  FOR  FINE  GRID 
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A Status  of  the  Activities  of  the  NASA/MSFC 
Pump  Stage  Technology  Team 


R.  Garcia,  R.  Williams,  and  Y.  Dakhoul 


The  Consortium  for  Computational  Fluid  Dynamics  (CFD) 
Application  in  Propulsion  Technology  was  established  to  aid  the 
transfer  of  CFD  related  advancements  among  academia,  government 
agencies,  and  industry.  The  specific  goals  of  the  Consortium 
are  to  develop  CFD  methodologies  necessary  to  solve  propulsion 
problems,  to  validate  these  methodologies,  and  to  apply  these 
methodologies  in  the  design  process.  To  accomplish  these 
goals,  a team  of  experts  in  various  related  fields  has  been 
formed,  a schedule  of  activities  necessary  to  meet  the  goals 
has  been  generated,  and  funding  for  the  activities  has  been 
obtained  from  NASA.  During  the  past  year  (3/91-3/92)  the 
team's  activities  have  focused  on  preliminary  code  validation 
and  on  the  design  of  an  advanced  impeller.  Six  codes  were  used 
to  calculate  the  flow  in  a Rocketdyne  0.3  flow  coefficient 
inducer  and  the  results  were  compared  to  L2F  data  available  for 
fho  inducer  _ This  ■iHonti  -F  i oH  Qhnrfrnim’  nrrc:  -In  -hhr* 

experimental  data  sets  and  in  the  analytical  solutions  which 
must  be  surmounted  in  any  future  team  activity.  The  design  of 
the  advanced  impeller  relied  heavily  on  CFD  results  to  obtain 
an  optimized  geometry.  The  optimized  geometry  has  been  ^ 

analyzed  using  four  different  codes  and  at  design  and 
off -design  conditions.  Activities  for  the  next  year  include 
the  optimization  of  a tandem  blade  impeller  design,  benchmark 
of  CFD  codes  for  diffuser  and  volute  flows,  the  collection  of 
L2F  data  for  "state-of-the-art"  impeller  and  inducer,  and  the 
verification  of  the  advanced  pump  team  impeller  design  in  a 
water  rig. 
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A Summary  of  the  Activities  of  the  NASA/MSFC 


Pump  Stage  Technology  Team 


R.  Garcia 

NASA,  Marshall  Space  Flight  Center 
R.  Williams 

NASA,  Marshall  Space  Flight  Center 
Y.  Dakhoul 

Sverdrup  Technologies 


Presented: 

CFD  Applications  Workshop 
MSFC,  April  28-30,  1992 


NASA 

National  Aeronautics  and 
Space  Administration 


A Summary  of  title  Activities  of  the  NASA/MSFC 
Pump  Stage  Technology  Team 


Overview 

• Structure/objectives 

• Approach 

• Validation  data 

• CFD  Analysis: 

- Benchmark  activity 

- Advanced  hardware 


Summary/conclusions 


176 


NASA  A Summary  of  the  Activities  of  the  NASA/MSFC 
National  Aeronautics  and  Pump  Stage  Technology  Team 

Space  Administration 


Structure/Objectives 


NASA 

National  Aeronautics  and 
Space  Administration 


A Summary  of  the  Activities  of  the  NASA/MSFC 
Pump  Silage  Technology  Team 


Approach 

* Assemble  a team  of  experts 

- Team  members  from  academia,  industry,  and  government  agencies 

* Implement  a plan  to  coordinate  pump  team  activities 

- Set  milestone  dates  consistent  with  rocket  engine  development  requirements 

* Hold  quarterly  meetings  to: 

- Critique  activities 

- Raise  unexpected/new  issues  and  requirements 

- Maintain  focus  on  the  deliverable  product  and  on  the  schedule 


178 


NASA  A Summary  of  the  Activities  of  the  NASA/MSFC 

National  Aeronautics  and  Pump  Stage  Technology  Team 

Space  Administration 


Pump  Team  Members 

Consortium  for  CFD  Application  in  Propulsion  Technology  Pump  Stage  Technology  Team 

• NASA  Marshall  Space  Right  Center  (MSFC) 

• NASA  Ames  Research  Center  (ARC) 

• NASA  Lewis  Research  Center  (LeRC) 

• David  Taylor  Research  Center 

• Rocketdyne  (RDYN) 

• Pratt  & Whitney  (P&W) 

• Aerojet 

• Ingersoll-Rand 

• Computational  Fluid  Dynamics  (CFD)  Research  Corporation 

• SECA 

• Scientific  Research  Associates  (SRA) 

• The  University  of  Alabama  in  Huntsville  (UAH) 

• Pennsylvania  State  University  (PSU) 

• University  of  Cincinnati 

• Virginia  Polytechnic  Institute 

• California  Institute  of  Technology 
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PUMP  DESIGN  TECHNOLOGY 
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PUMP  DESIGN  TECHNOLOGY 
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NASA  A Summary  of  the  Activities  of  the  NASA/MSFC 

National  Aeronautics  and  Pump  Stage  Technology  Team 

Space  Administration 


Pump  CFD  Code  Validation  Tests  SMSME  HPFTP  Impeller 
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NASA  A Summary  of  the  Activities  of  the  NASA/MSFC 

National  Aeronautics  and  Pump  Stage  Technology  Team 

Space  Administration 


CFD  Analysis 


• Advanced  hardware  development 

- Conventionally  designed  advanced  impeller  optimized  using  CFD 

• Impeller  design  to  satisfy  STME  fuel  pump  requirements  with  two  stages 

• CFD  study  of  1 5 parameters:  b2,  P2,  axial  length,  total  wrap  angle,  and  discharge  wrap 
angle  difference  (hub-to-tip) 

• Viability  of  CFD  parametrics  demonstrated 

• Baseline  and  optimized  geometry  analyzed  by  five  team  members 

• All  solutions  show  higher  efficiency  and  reduced  impeller  discharge  flow  distortion 

• Off-design  analysis  under  way 

• Impeller  being  manufactured;  performance  to  be  verified  in  water  rig  in  the  fall  of  1 992 

- Tandem  blade  impeller  concept 

• Concept  has  potential  for  increased  head  coefficient  and  efficiency 

• CFD  parametric  study  to  begin  in  May  1992 

• Study  to  include  position  of  blade  split,  blade  clocking,  and  chordwise  spacing 

• Final  configuration  to  rely  entirely  on  results  of  parametric  study 

• Impeller  will  be  sized  to  satisfy  STME  fuel  pump  requirements 
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Organization 
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OPTIMIZED  IMPELLER:  CM  VS.  X 
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AMES  RESULTS:  EXIT  CM  VS.  X 
FOR  VARYING  FLOW  RATE 
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NASA  A Summary  of  the  Activities  of  the  NASA/MSFC 
National  Aeronautics  and  Pump  Stage  Technology  Team 

Snare  Administration  1 ^ & * 


Summary/Conduslons 

• Technology  team  in  place  and  functioning  efficiently 

- Participation  by  industry,  universities,  and  government 

• Detailed  experimental  data  sets  suitable  for  benchmarking  have  been  or  are  being  generated 

• Preliminary  evaluation  of  six  different  codes  complete 

• CFD  codes  being  used  to  reduce  the  design  development  time  and  improve  performance  of 
advanced  impellers 

• Verification  of  advanced  impeller  predictions  planned  for  the  fall  of  1 992 

• Future  work  to  include  impeller-diffuser  interaction  and  inducer  non-cavitating  analysis 
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CFD  ANALYSIS  OF  PUMP  CONSORTIUM  IMPELLER 
Gary  C.  Cheng*,  Y.S.  Chen1,  and  R.W.  Williams1 
Abstract 

Current  design  of  high  performance  turbopumps  for  rocket  engines  requires  effective  and 
robust  analytical  tools  to  provide  design  impact  in  a productive  manner.  The  main  goal  of  this 
study  is  to  develop  a robust  and  effective  computational  fluid  dynamics  (CFD)  pump  model  for 
general  turbopump  design  and  analysis  applications.  A Navier-Stokes  flow  solver,  FDNS, 
embedded  with  the  extended  k-e  turbulence  model  and  with  appropriate  moving  interface 
boundary  conditions,  is  developed  to  analyze  turbulent  flows  in  the  turbomachinery  devices.  The 
FDNS  code  has  been  benchmarked  with  its  numerical  predictions  of  the  pump  consortium 
inducer,  and  provides  satisfactory  results.  In  the  present  study,  a CFD  analysis  of  the  pump 
consortium  impeller  will  be  conducted  with  the  application  of  the  FDNS  code.  The  pump 
consortium  impeller,  with  partial  blades,  is  the  new  design  concept  of  the  advanced  rocket  engine. 
A 3-D  flow  calculation  with  81  x 41  x 41  grid  system  was  conducted  for  the  team  base-line 
impeller.  The  result  shows  a massive  flow  separation  occurs  between  the  full-blade  pressure 
surface  and  the  partial-blade  suction  surface.  Similar  result  was  predicted  by  the  other 
consortium  members.  A pump  consortium  optimized  impeller,  a revision  based  on  the  base-line 
impeller,  was  then  designed  by  Rocketdyne  to  remove  the  flow  separation.  A 3-D  flow  analysis, 
with  103  x 23  x 30  mesh  system  and  with  the  inlet  flow  conditions  provided  by  Rocketdyne,  was 
performed  for  the  optimized  impeller.  The  numerical  result  indicates  no  flow  separation  occurs 
inside  the  flow  passage,  which  is  also  consistent  with  the  other  consortium  members’  predictions. 
However,  the  flow  field  inside  the  optimized  impeller  as  calculated  by  the  team  members  showed 
great  variations,  especially  near  the  exit  shroud  region.  The  discrepancy  is  suspected  to  h<*  due 
io  different  exit  boundary  conditions  used  by  the  consortium  members.  Therefore,  three  different 
exit  wall  boundary  conditions  will  be  further  examined  by  the  FDNS  code,  those  are  fixed-wall, 
wall-slip  (symmetry),  and  rotating  wall  boundary  conditions.  The  computed  results  will  be 
compared  in  order  to  address  the  effect  of  exit  boundary  conditions  on  the  impeller  flow  field. 
Meanwhile,  two  off-design  cases  of  the  optimized  impeller,  80%  and  120%  of  the  design  flow, 
will  also  be  analyzed  with  a particular  exit  boundary  condition.  All  CFD  analysis  of  the  pump 
consortium  base-line  impeller,  and  the  optimized  impeller  with  various  exit  boundary  conditions 
will  be  presented  in  the  coming  CFD  workshop  meeting. 


SECA,  Inc.,  3313  Bob  Wallace  Ave.,  Suite  202,  Huntsville,  AL 
Engineering  Sciences,  Inc.,  4920  Corporate  Dr.,  Suite  K,  Huntsville,  AL 
ED  32,  NASA/Marshall  Space  Right  Center,  Huntsville,  AL 
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CFD  ANALYSIS  OF  PUMP  CONSORTIUM  IMPELLER 
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INLET/EXIT  WALL  B.C.  TESTED 
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DEFINITION  OF  PERFORMANCE  PARAMETERS 


• Cu  = cu  / Utip;  CM  = cM  / Utip  where  cu  = Absolute  Tangential  Velocity, 
cM  = Meridional  Velocity,  Utlp  = Wheel  Tip  Velocity 

• p = Relative  Flow  Angle  Relative  to  Tangential  Direction 

• Relative  Radius  = (R,  - Rhub)  / (Rshroud  - RhJ 

• Relative  X = (X,  - Xshroud)  / 

(^hub  “ ^shroud) 

® Relative  Angle  = (Anglej  ~ Anglesuctj0n)  / (Anglepressure  - Anglesuctj0n) 

• ¥ (Head  Coefficient)  = AHtg  / Utip2 

• t|  (Efficiency)  = Head  Rise  / Euler  Head  Rise 
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VELOCITY  VECTORS  NEAR  PRESSURE  SIDE  OF  SPLITTER 


Case  2 Case  5 
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CM/U.  CM/U 


OPTIMIZED  IMPELLER:  BLADE-TO-BLADE  CM 
AT  THE  IMPELLER  EXIT,  NEAR  THE  SHROUD 
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OPTIMIZED  IMPELLER:  BLADE-TO-BLADE  CM 
AT  IMPELLER  EXIT,  NEAR  THE  HUB 
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HEAD  COEFFICIENT  EFFICIENCY 


OPTIMIZED  IMPELLER: 
EFFICIENCY  VS.  X 
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CONCLUSIONS 


• THE  PRESENT  CFD  RESULTS  HAVE  SHOWN  SENSITIVITY 
OF  INLET  AND  EXIT  WALL  BOUNDARY  CONDITIONS  ON 
THE  FLOW  STRUCTURE  INSIDE  THE  OPTIMIZED 
CONSORTIUM  IMPELLER  DESIGN 

• INLET  SHROUD  WALL  BOUNDARY  TREATMENTS  HAVE 
SIGNIFICANT  EFFECT  ON  THE  FLOW  SPLIT  AROUND 
THE  PARTIAL  BLADE  (MORE  FLOW  THROUGH  THE 
PARTIAL/FULL-PRESSURE  PASSAGE  WHEN  THE  INLET 
SHROUD  WALL  IS  ASSUMED  ROTATING) 

• ONLY  MINOR  IMPACT  ON  THE  OVERALL  IMPELLER 
PERFORMANCE  DATA  WAS  REVEALED  FOR 
DIFFERENT  BOUNDARY  CONDITIONS  IMPOSED 


I 

I 

N92-32288  ' 


Abstract  of  a proposed  paper  for  the  presentation  at  workshop 
for  CFD  Applications  in  Rocket  Propulsion  to  be  held  at 
NASA  Marshall  Space  Flight  Center,  AL,  April  28-30,  1992 

i 

i CFD  APPLICATIONS  IN  PUMP  FLOWS 

! Cetin  Kiris,  Liang  Chang 

MCAT  Institute,  Moffett  Field,  CA 

and 

Dochan  Kwak 

NASA- Ames  Research  Center,  Moffett  Field,  CA 

I 

I The  objective  of  the  proposed  paper  is  to  develop  a computational  procedure 

that  solves  incompressible  Navier-Stokes  equations  for  pump  flows.  The  solution 
method  is  based  on  the  pseudocompressibility  approach  and  uses  an  implicit-upwind 
differencing  scheme  together  with  the  Gauss-Seidel  line  relaxation  method.  The  equa- 
tions are  solved  in  steadily  rotating  reference  frames  and  the  centrifugal  force  and  the 
Coriolis  force  are  added  to  the  equation  of  motion.  As  a bench  mark  problem,  the 
flow  through  the  Rocketdyne  inducer  is  numerically  simulated.  A coarse  grid  solution 
is  obtained  with  a single  zone  by  using  an  algebraic  turbulence  model.  In  multi-zone 
fine  grid  computation,  one-equation  Baldwin-Barth  turbulence  model  is  utilized.  Nu- 
merical results  are  compared  with  experimental  measurements  and  a good  agreement 
is  found  between  the  two.  The  resulting  computer  code  is  then  applied  to  the  flow 
analysis  inside  two-stage  fuel  pump  impeler  operating  at  80  %,  100  %,  and  120  % of 
design  flow. 
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CFD  APPLICATIONS  IN  PUMP  FLOWS 


Cetin  Kiris,  Dochan  Kwak,  and  Leon  Chang 
NASA- Ames  Research  Center 


Workshop  for  CFD  Applications  in  Rocket  Propulsion 
NASA-MSFC,  April  28-30,  1992 


Introduction 


• Motivation 

Increasing  efficiency  and  reliability  of  the  liquid  rocket  engine 
components  is  an  important  task. 

> Understand  fluid  dynamics  of  fuel  and  oxidizer  flows  from  fuel 
tank  to  plume. 

> Role  of  CFD  toward  a better  design. 

• Goal 

To  implement  CFD  technology  to  simulate  the  flow  through  the 
pump  components. 

First  Step  : Bench  mark  problems  and  component  analysis. 
Second  Step  : Unsteady  flows  through  the  entire  pump  (future 
work) . 
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Method  of  Solution 


• Available  algorithms  for  pump  applications  are  : INS3D-UP  and 
IND3D-LU. 

• Currently  INS3D-UP  is  used. 

• Based  on  method  of  artificial  compressibility. 

• Both  steady-state  and  time-accurate  formulation. 

• Steady-state  formulation  in  steadily  rotating  reference  frame. 

• Multi- Zone  and  Operlapped  grid  scheme  capability. 

• Computing  time:  ss  lxlO-4  sec/grid  point/iteration 

• Memory  Usage:  ss  45  times  number  of  grid  points  in  words 
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INS3D-UP  Algorithm 


• Central  differencing  for  viscous  fluxes 

• Upwind  differencing  for  convective  fluxes 

^>>  3rd  and  5th  order  flux-difference  splitting  is  used  for  the  right 

hand  side  terms 

• Gauss-Seidel  line  relaxation  relaxation 

• Unlimited  time  step  usage  in  steady-state  formulation. 

• Inflow  and  Outflow  boundaries  based  on  Method  of  Characteristics 

» Inflow  Boundary  : Three  velocity  components  specified 
> Outflow  Boundary  : Static  pressure  specified 

• Quasi-implicit  boundary  conditions  at  zonal  interfaces. 
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Steady-State  Formulation 

i 

• Introduce  artificial  compressibility  term  to  the  continuity  equation 


«E  = _»  (w  + s±  + ?E] 

dr  y d£  dr)  d(  J 


> (3  is  an  artificial  compressibility  constant 

> r is  a pseudo-time  step 

^>5  is  a source  term  as  centrifugal  and  coriolis  forces. 

• Euler  Implicit  time  discretization 

• Solve  system  of  equations  iteratively  in  pseudo-time  until  solution 
converges  to  a steady  state 
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Time- accurate  Formulation 


Discretize  the  time  term  in  momentum  equations  using  second-order 
three-point  backward-difference  formula 


OU  dV  dW 
+ di 1 + d( 


n+1 


= o 


3 qn+1  - 4 qn  + q 
2A  t 


n— 1 


_ _^n+ 1 


Introduce  a pseudo-time  level  and  artificial  compressibility 


- (pn+i?m+i  — pn+1’m)  = — /?V  • gn+1’m+1 
At 


1.5 

A/ 


^n+l,ra+l  _ _^n+l,m+l  _ 


3^n+l,m  _ 4^n  + £n-l 

2Ai 


Iterate  the  equations  in  pseudo-time  for  each  time  step  until 
incompressibility  condition  is  satisfied. 
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Inducer  Computations 

i 

• Grid  size  : 187  x 27  x 35 

• Baldwin-Lomax  algebraic  turbulence  model. 

• Tip  clearance  region  is  included. 

• Computer  time  : 5 - 6 Cray-YMP  hours 

• Multi-zone  computation  (currently  underway) 

• Grid  1 : 63  x 37  x 74  Upstream  of  inducer 

• Grid  2 : 115  x 37  x 48  Inducer  blades 

• Grid  3 : 51  x 37  x 20  Bull-nose  cavity 

• Grid  4 : 51  x 37  x 49  Downstream  of  blades 

• One-equation  Baldwin-Barth  turbulence  model. 


VALIDATION  OF  INS3D-UP 

Pump  Design  Code  Development  for  STME 


K3 

K) 

'*-4 


Rocketdyne  inducer  geometry 


Circumferential  Angle  from  Suction  Side  (Degree) 


Comparison  of  relative  total  velocities 


A BCD 


Schematic  of  the  experimental 
measurement  planes 


Pressure 


Computed  surface  pressure  for 
Rocketdyne  inducer 
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Impeller  Water  Test  Conditions 


• Number  of  Blades 

• Design  Speed,  RPM 

• Design  Flow,  GPM 

• Reynolds  Number,  per  inch 

• Inlet  Tip  Diameter,  inch 

• Inlet  Hub  Diameter,  inch 

• Outflow  Diameter,  inch 


6 full,  6 partial 
6632 
1205 
1.81e+5 
6.0 
3.9 
9.045 


Discharge  Tip  Speed,  fps 


249.5 


NLS  Pump  Impeller 


Surface  Pressure 


Pressure 

-1811.000 
-1711.000 
-1711.000 
-17  1.000 
- 1 ooo 

-1711.000 
-16U.000 
-161 .000 
-16-1.000 
-16: .000 
-16t .000 

— 1 51  .000 
-15* .000 
-15- .000 
-152 .000 
-15C .000 

- 1 -it . uuu 
- 1 it . uuu 

-Mi.  UUU 

- 1 -I2.00U 
-1 40.000 
-138.000 
-136.000 
-134.000 
-132.000 
-130.000 
-128,000 
-126.000 
-124  000 
-122  000 
-120  000 


Relative  Velocity  Vectors  Colored  by  Pressure 


Shroud 


Huh 


Old  Design 


New  Design 
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NLS  Pump  Impeller 

120%  100% 


Particle ir  are  related  near  shroud 


120%  of  Design  Flow 


ORIGINAL  r. 

OF  POOR  QUAUT 


Velocity  Vectors  Colored  by  Pressure 
upstream  of  discharge  plane 
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Velocity  Vectors  Colored  by  Pressure 
upstream  of  discharge  plane 
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-0.3-1250 


Velocity  Vectors  Colored  by  Pressure 
upstream  of  discharge  plane 
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IMPELLER  EXIT  CM  DISTRIBUTION 


100%  of  Design  Row 

120%  of  Design  Row 
80%  of  Design  Row 
-A-  60%  of  Design  Row 


0.2 


0.4 


Relative  X (from  shroud  to  hub) 


BETA 


IMPELLER  EXIT  BETA  DISTRIBUTION 


100%  of  Design  Flow 

120%  of  Design  Flow 
- •-  80%  of  Design  Flow 
-A-  60%  of  Design  Flow 


Relative  X (from  shroud  to  hub) 
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0.25 


CM  @ R/Rtip  =1.0  120%  of  Design  Flow 


0.20 


0.15 


x-0.0  @ shroud 

x-,10 
- - X-.30 

X-.50 

-A-  x-,70 
x-,90 


0.10 


0.05 


0.00 


0.0 


Relative  Angle  (from  suction  side  to  pressure  side) 
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IMPELLER  OVERALL  PARAMETERS 


Design 

Baseline 

Optimized 

Optimized 

Optimized 

Optimized 

Opt. 

Downstream 

no-slip  wall 

no-slip  wall 

slip  b.c. 

slip  b.c. 

slip  b.c. 

slip  b.c. 

Boundary 
Design  Flow 

100% 

100% 

100% 

120% 

80% 

60% 

Efficiency 

.946 

.983 

.943 

.931 

.955 

.963 

Head 

.663 

.636 

.659 

.631 

.667 

.691 

Coefficient 
Relative  Flow 

25.5 

18.3 

19.8 

21.1 

19.2 

17.6 

Angie,  Deg 
Absolute  Flow 

7.01 

6.44 

6.63 

7.31 

4.33 

3.10 

Angle,  Deg 
Meridional 

24.5 

21.1 

19.4 

24.04 

15.52 

11.52 

Velocity 
Flow  Split 

.45/.55 

.47/.53 

.437.57 

.435/.565 

.41 5/585 

.377.63 

Summary 


• The  present  capability  to  compute  a 3-D  flow  through  a complex 
internal  geometry  is  demonstrated.  Advanced  impeller  analysis  showed 
that  overall  parameters  do  not  have  a significant  change  between  80%, 
100%,  and  120%  of  design  flow. 

» 

• Solution  algorithm  was  tested  and  validated  in  model  problems. 
Inducer  computations  indicate  less  than  10  % error  in  velocities.  Tip 
and  wake  regions  show  biggest  discrepancies.  Future  - work  will  be 
focused  on  turbulence  modeling. 

• Solution  procedure  obtained  here  can  be  used  in  the  design  process  of 
pump  components. 


N92-32289  * 

COMPUTATION  OF  THE  FLOW  FIELD  IN  A CENTRIFUGAL  IMPELLER  WITH  SPUTTER  BLADES* 

Frederik  J.  de  Jong,  Sang-Keun  Choi,  T.R.  Govindan  and  Jayant  S.  Sabnis* 

Scientific  Research  Associates,  Inc. 

Glastonbury,  CT 

ABSTRACT 

To  support  the  design  effort  of  the  STME  Fuel  Pump  Stage,  viscous  flow  calculations  have  been 
performed  in  a centrifugal  impeller  with  splitter  blades.  These  calculations  were  carried  out  with 
SRA’s  Navier-Stokes  solver  (MINT),  which  employs  a linearized  block-implicit  ADI  procedure  to 
iteratively  solve  a finite  difference  form  of  the  system  of  conservation  equations  of  mass,  momentum, 
and  energy  in  body-fitted  coordinates.  A computational  grid  was  generated  algebraically  for  the 
"channel'  between  two  main  blades  of  the  impeller  and  extended  both  upstream  of  the  impeller  inlet 
and  downstream  of  the  impeller  exit  so  that  the  appropriate  boundary  conditions  could  be  applied 
(viz.  specified  velocity  profiles  at  the  inflow  boundary,  and  specified  pressure  at  the  outflow 
boundary). 

The  results  of  the  calculations  show  that  although  the  overall  level  of  flow  distortion  near  the 
impeller  exit  is  not  very  large,  there  is  a noticeable  difference  between  the  flow  patterns  in  the  two 
"passages"  (one  passage  between  the  pressure  side  of  the  fun  blade  ar.d  the  suction  side  ot  the 
splitter  biade,  and  the  other  one  between  the  pressure  side  of  the  splitter  blade  and  the  suction  side 
of  the  next  full  blade).  For  example,  the  pressure  distribution  shows  that  the  splitter  blade  is  loaded 
less  heavily  than  the  main  blade.  At  the  same  time,  the  flow  distortion  near  the  suction  side  of  the 
main  blade  is  larger  than  that  near  the  suction  side  of  the  splitter  blade.  A better  understanding  of 
these  results  can  be  obtained  by  studying  "particle  traces"  (streamlines  in  a frame  of  reference  fixed 
to  the  rotating  impeller).  These  traces  show  that  a significant  amount  of  low  momentum  fluid 
(originating  from  the  hub  and  shroud  boundary  layers)  moves  from  the  pressure  side  to  the  suction 
side  in  the  impeller  "channel"  ahead  of  the  splitter  blade,  and  ends  up  in  the  passage  between  the 
pressure  side  of  the  splitter  blade  and  the  suction  side  of  the  full  blade.  This  explains  why  the  region 
of  flow  distortion  in  this  passage  is  larger  than  that  in  the  other  passage,  and  why  the  mass  flow 
through  this  channel  is  less.  The  understanding  of  the  physics  of  impeller  flow  fields  that  results  from 
analyzing  viscous  flow  calculations  such  as  the  one  described  above  is  very  valuable  in  pump  stage 
design. 

* This  work  was  supported  by  NASA  Marshall  Space  Flight  Center  under  Contract  NAS8-38866. 

* Currently  at  United  Technologies  Research  Center,  East  Hartford,  CT 
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PRECEDING  PAGE  BLA.'iK  iW  FILMED 


OBJECTIVES 


DEVELOPMENT  OF  A PREDICTIVE  TOOL  FOR  THE  VISCOUS  FLOW 


IN  CENTRIFUGAL  IMPELLERS 


SUPPORT  OF  DESIGN  EFFORT  OF  STME  FUEL  PUMP  STAGE 


ENHANCEMENT  OF  UNDERSTANDING  OF  PHYSICAL  FLOW  PHENOMENA 


APPROACH 


EQUATIONS  SOLVED: 

REYNOLDS-AVERAGED  NAVIER-STOKES  EQUATIONS  IN 
ROTATING  BODY-FITTED  COORDINATES 

METHOD  OF  SOLUTION: 

LINEARIZED  BLOCK  IMPLICIT  (ADI)  SCHEME 

TURBULENCE  MODEL: 

MIXING  LENGTH  OR  TWO-EQUATION  (K-e) 


Sdantfc 

Research 

Associates 
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BOUNDARY  CONDITIONS 


INFLOW: 

VELOCITY  PROFILES,  STAGNATION  TEMPERATURE 
OUTLFLOW: 

CONSTANT  PRESSURE 
WALLS: 

NO-SLIP,  ADIABATIC 

UPSTREAM  AND  DOWNSTREAM  MAIN  BLADE  EXTENSIONS: 
PERIODICITY 


WATER  TEST  CONDITIONS 


DIMENSIONS 


INLET  TIP  DIAMETER 

6.0 

INCH 

(= 

0.1524 

M) 

INLET  HUB  DIAMETER 

3.9 

INCH 

(= 

0.0991 

M) 

EXIT  TIP  DIAMETER 

D 

9.045 

INCH 

(= 

0.2297 

M) 

FLOW  CONDITIONS 


DESIGN  SPEED 

CO 

6322 

RPM 

(= 

662 

RAD/S) 

TIP  SPEED 

u 

249.5 

FT/S 

(= 

76.05 

M/S) 

DESIGN  FLOW 

Q 

1205 

GPM 

(= 

0.0760 

M3/S) 

AVERAGE  INFLOW  VELOCITY 

23.7 

FT/S 

(= 

722 

M/S) 

FLOW  COEFFICIENT 

♦ 

0.095 

REYNOLDS  NUMBER 

Q/u  D 

5.5  104 

uD/u 

2.9  10* 

Scwtiffc 

Rmmvc/) 

Associate* 
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IMPELLER  GEOMETRY 

MERIDIONAL  SURFACE 


Scientific 

Research 

Associates 


PRESSURE  IN  MERIDIONAL  SURFACE 

BETWEEN  SPLITTER  PRESSURE  SIDE  AND  BLADE  SUCTION  SIDE 
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I 

I 


VELOCITY  MAGNITUDE  IN  MERIDIONAL  SURFACE 

BETWEEN  BLADE  PRESSURE  SIDE  AND  SPLITTER  SUCTION  SIDE 


VELOCITY  MAGNITUDE  IN  MERIDIONAL  SURFACE 

BETWEEN  SPLITTER  PRESSURE  SIDE  AND  BLADE  SUCTION  SIDE 
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STREAMWISE  VELOCITY  IN  MERIDIONAL  SURFACE 

BETWEEN  BLADE  PRESSURE  SIDE  AND  SPLITTER  SUCTION  SIDE 


STREAMWISE  VELOCITY  IN  MERIDIONAL  SURFACE 

BETWEEN  SPLITTER  PRESSURE  SIDE  AND  BLADE  SUCTION  SIDE 


IMPELLER  GEOMETRY 

MID-SPAN  BLADE-TO-BLADE  SURFACE 
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Associates 


PRESSURE 


VELOCITY  MAGNITUDE 


Scientific 
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STREAMWISE  VELOCITY 
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IMPELLER  GEOMETRY 

CROSS-SECTION  NEAR  TRAILING  EDGE 
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VELOCITY  MAGNITUDE 
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STREAMWISE  VELOCITY 
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SUMMARY 


NAVIER-STOKES  PROCEDURE  HAS  BEEN  USED  TO 
SIMULATE  THE  FLOW  FIELD  IN  A CENTRIFUGAL 
IMPELLER  WITH  SPLITTER  BLADES 

ANALYSIS  OF  THE  RESULTS  PROVIDES  UNDERSTANDING 
OF  THE  PHYSICAL  FLOW  PHENOMENA 

• VALUABLE  IN  PUMP  STAGE  DESIGN 

• USEFUL  TO  GUIDE  CFD  CODE  DEVELOPMENT 
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Axsocktas 


• V 


N92-32290 


IMPELLER  TANDEM  BLADE  STUDY  WITH  GRID  EMBEDDING 
FOR  LOCAL  GRID  REFINEMENT 


George  Bache’ 


Aerojet  Propulsion  Division 
Bldg.  2019;  Dept.  5236 
P.O.  Box  13222 
Sacramento  CA  95813-6000 


ABSTRACT 


Flow  non-uniformity  at  the  discharge  of  high  power  density  impellers  can  result  in 
significant  unsteady  interactions  between  impeller  blades  and  downstream 
diffuser  vanes.  These  interactions  result  in  degradation  of  both  performance  and 
pump  reliability.  The  MSFC  Pump  Technology  Team  has  regognized  the 
importance  of  resolving  this  problem  and  has  thus  initiated  the  development  and 
testing  of  a high  head  coefficient  impeller.  One  of  the  primary  goals  of  this 
program  is  to  improve  impel  er  performance  and  discharge  flow  uniformity.  The 
objective  of  the  present  work  is  complimentary.  Flow  uniformity  and  performance 
gains  were  sought  through  the  application  of  a tandem  blade  arrangement.  The 
apDroach  arinptpH  was  tc  numerically  establish  fiuw  characteristics  at  the  impeller 
discharge  for  the  baseiine  MSFC  impeller  and  then  parametrically  evaluate 
tandem  blade  configurations.  A tandem  design  was  sought  that  improves  both 
impeller  performance  and  discharge  uniformity.  The  Navier-Stokes  solver 
AEROVISC  was  used  to  conduct  the  study.  Grid  embedding  is  used  to  resolve 
local  gradients  while  attempting  to  minimize  model  size.  Initial  results  indicate  that 
significant  gains  in  flow  uniformity  can  be  achieved  through  the  tandem  blade 
concept  and  that  blade  clocking  rather  than  slot  location  is  the  primary  driver  for 
flow  uniformity. 
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IMPELLER  TANDEM  BLADE  STUDY  WITH  GRID  EMBEDDING 

FOR  LOCAL  GRID  REFINEMENT 


GEORGE  BACHE’ 
AEROJET  PROPULSION  DIVISION 


10th  WORKSHOP  FOR  COMPUTATIONAL  FLUID  DYNAMIC 
APPLICATION  IN  ROCKET  PROPULSION 


APRIL  28-30,  1992 


Propulsion  Division 


GenCorp 
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IMPELLER  TANDEM  BI.ADE  STUDY  WITH  GRID  EMBEDDING 
FOR  LOCAL  GRID  REFINEMENT 


George  Bache’ 


Aerojet  Propulsion  Division 
Blda.  2019;  Dept.  5236 
P.O.  Box  13222 
Sacramento  CA  95813-6000 


ABSTRACT 


Flow  non-uniformity  at  the  discharge  of  high  power  density  impellers  can  result  in  significant 
unsteady  interactions  between  impeller  bludes  and  downstream  diffuser  vanes.  These  interactions 
result  in  degradation  of  both  performance  and  pump  reliability.  The  MSFC  Pump  Technology  Team 
has  regognized  the  importance  of  resolving  this  problem  and  has  thus  initiated  the  development  and 
testing  of  a high  head  coefficient  impeller.  One  of  the  primary  goals  of  this  program  is  to  improve 
impeller  performance  and  discharge  flow  uniformity.  The  objective  of  the  present  work  is 
complimentary.  Flow  uniformity  and  performance  gains  were  sought  through  the  application  of  a 
tandem  blade  arrangement.  The  approach  adopted  was  to  numerically  establish  flow  characteristics 
at  the  impeller  discharge  for  the  baseline  MSFC  impeller  and  then  parametrically  evaluate  tandem 
blade  configurations.  A tandem  design  was  sought  that  improves  both  impeller  performance  and 
discharge  uniformity.  The  Navier-Stokes  solver  AEROVISC  was  used  to  conduct  the  study.  Grid 
embedding  is  used  to  resolve  local  gradients  while  attempting  to  minimize  model  size.  Initial  results 
indicate  that  significant  gains  in  flow  uniformity  can  be  achieved  through  the  tandem  blade  concept 
and  that  blade  clocking  rather  than  slot  location  is  the  primary  driver  for  flow  uniformity. 
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GenCorp 

AEROJET 


MOTIVATION  FOR  WORK 


Propulsion  Division 


HIGH  POWER  DENSITY  IMPELLERS  USED  IN  ROCKET  PROPUSLION  TURBOPUMPS  PRODUCE 
SIGNIFICANT  UNSTEADY  INTERACTIONS  BETWEEN  IMPELLER  BLADES  AND  DOWNSTREAM 
DIFFUSER  VANES. 


- PERFORMANCE  DEGRADATION 


REDUCED  COMPONENT  LIFE  (RELIABITY) 


GOAL: 


USE  TANDEM  BLADE  CONCEPT  TO  IMPROVE  IMPELLER  PERFORMANCE  AND/OR 
DISCHARGE  FLOW  UNIFORMITY 


APPROACH 


Propulsion  Division 


GenCorp 


- PREDICT  FLOW  IN  MSFC  PUMP  TECHNOLOGY  TEAM  BASELINE  IMPELLER 


- NUMERICALLY  STUDY  SEVERAL  TANDEM  BLADE  CONFIGURATIONS 

ho 

On 


- PARAMETRICALLY  VARY  SLOT  LOCATION  AND  CLOCKING 


- PERFORM  NUMERICAL  PARAMETRICS  WITH  AEROVISC 


JUDGE  IMPROVEMENTS  BASED  ON  PERFORMANCE  AND  DISCHARGE  FLOW  UNIFORMITY 
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Aerovisc  Numerics 


Propulsion  Division 


BenCorp 


• Formulation 

- Reynolds  Stress  Averaged  Navier-Stokes  Equations 

- Cartesian  Primitive  Variable  Approach 

- Strongly  Conservative,  Colocated  Form 

- k-e  and  ARS  Turbulence  Models  With  Wall  Functions 

- 2 Layer  Turbulence  Model 

• Discretization 

- Flux  Element  Based  Finite  Volume  Method 

- General  Non-Orthogonal  Boundary-Fitted  Structured  Grid 

- Advection  Schemes  Have  Two  Components 

- Upwind  Skew  Scheme  (UDS,  MWS,  LPS)  ==>  Transverse  Gradients 

- Physically  Based  Correction  Term  (PAC)  ==>  Streamwise  Gradients 

- Rhie  Type  4th  Order  Pressure  Redistribution 

- Pressure  /Velocity  Coupling  For  I.C.  Flows 


• Solver 

- Vectorized  Gauss-Siedel  or  Incomplete  Cholesky  Base  Solver 

- Multigrid  (Large  Grids)  and  Block  Correction  (High  Aspect  Ratio  Grids) 


GenCorp  Demonstrated  Code  Capabilities  Propulsion  Division 

AEROJET 

• Applicable  Flow  Ranges 

- Incompressible 

- Subsonic,  Transonic  and  Supersonic 

- Laminar,  Turbulent  or  I nviscid 

- Multi-Component 

- 2-Phase  (Solid  Particle  / Gas) 

• Coriolis  and  Centrifugal  Terms  For  Turbomachinery  Applications 

• Conjugate  Heat  Transfer  or  Specified  Wall  Temperature  / Flux 

• Flexible  Geometric  Modeling  Features 

- Arbitrary  Periodicity 

- Multiple  Blocked  Regions 

- Grid  Embedding  or  Attaching 

• Fixed,  Moving  or  Rotating  Walls 

• Variable  Fluid  and  / or  Sol  id  Properties 
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Solid  Model  of 
Eckardt  Compressor 


Propulsion  Division 


The  Krain  Impeller  - Introduction 


• 24  Blades,  22,363  rpm,  4.7 
Design  Pressure  Ratio,  4.0  kg/s 
Total  Mass  Flow 

• Computational  Domain  Includes 
Inlet,  Tip  Region,  and  Diffuser 

• Absolute  Frame  Total  Pressure 
and  Total  Temperature  at  Inlet, 
Mass  Flow  Specified  at  Outlet, 
Log-Law  Used  at  Walls 

• 80,000  Nodes,  Second-Order- 
Accurate  Skew  Upwind  Scheme, 
Coupled  Multigrid 

• Note:  Inlet  Geometry  Was 

Estimated 


IRAD-7 


The  Krain  Impeller  - Comparisons  to  Data 


• Total  Run  Time  About  150  Hours  on  a Personal  Iris  4D/25 

• Measured  Total  Pressure  Ratio  Was  About  4.1,  the  Code  Predicted 

4.26  3.0i 


■f '..y'..-.. i.  ' »■  m,tm% 

0.0  0.5  1.0 

X/SM  — 

Meridional  Velocity  Circumferentially  Averaged 

Components  at  Outlet  Data  Shroud  Static  Pressure 

(Left)  Calculation  (Right) 


IR4D-7 


uhi  i an  r i on  inoucfr 
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TOTAL  RELATIVE  VELOCITY  (FT/SEC) 


INDUCER  CODE  VALIDATION 
PLANE  B,  X-2.725  IN.,  ANGLE=-133.755  DEG 
POSITION  5,  RADIUS-  2.365  IN. 


LEQEND 
• DATA 
8ECA 


RKDN 


LERC 


AERO 


8RA 


CIRCUMFERENTIAL  ANGLE  (DEQ.) 


INDUCEFI  CODE  VALIDATION 
PLANE  B,  X=2.725  IN.,  ANGLE=-133.755  DEG. 
POSITION  5,  RADIUS-  2.365  IN. 
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Propulsion  Division 


• FLUID  HYDROGEN 

• 6 FULL  + 6 PARTIAL 

• N - 30108  RPM 

• TIP  DIAM.  = 14.14  in. 

• TIP  SPEED  = 1857  fps 

• OUTLET  BLADE  ANGLE  = 49.5  deg. 

• SPECIFIC  SPEED  = 1141 

• HEAD  COEFF.  = 0.57? 


MSFC  PUMP  CFD  CONSORTIUM  IMPELLER 
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Propulsion  Division 


MSFC  PUMP  CFD  CONSORTIUM  IMPELLER 
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4ERDJET 

Propulsion  Division 


BASELINE  IMPELLER.  VISCOUS  SOLUTION 


O.  -250  -500  .750 

I HETfll  1 J ''DLLIh  IHLMi 


Utip-1858  fps  , MUM)  Deg  . mass  flow  split  45%  - 55 
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Gencorp 

AEROJET 

Propulsion  Division 


IMPELLEF  DOWNSTREAM  PROFILE 

BASFI  INF  CONFIGURATION  34K  NODES 


r / rtip  = 0.99 
avrg  = 8 deg 


SHROUD 


r / rtip  = 1 .05 
avrg  - 9 deg 


r / rtip  r i.io 
avrg  - 9 deg. 


(>  f Pt  I 
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□encorp 

>)ERDJcT  IMPELLER  DOWNSTREAM  PROFILE 

Propulsion  Division  baseunf  configuration  34k  nodes 


r i rtip  = 1.10 
avrg  = 1389  psi  j 

/... 


wn 


Propulsion 


i\j 

00 


FULL 
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CONFIGURATIO 


ORIGINAL  PA^i 

280  OF  POOR  QUALITY 


□enCdrp 


AEROJET 

Propulsion  Division 


IMPELLER  DISCHARGE  CIRC.  Cm  DISTRIBUTION 
Mass  Averaged  Hub  to  Tip 


IHEim  I I 0 L L T n THETP 


' R/Rtip  - 1.0? 


ORIGINAL  PAGE  IS 

OF  POOR  QUALfTY 
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Propulsion  Division 


IMPFI.  LER  DISCHARGE  CIRC.  Cm  DISTRIBUTION 
Mass  Averaged  Hub  to  Tip 


' 1 • . .--m  l c,.\ 

THFTfli  T 1 TlFI  TD  THFTP  T H F T R i T i / R F I T R T H F T R 


R/Rtip  1.05 
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IMPFLLER  DISCHARGE  CIRC.  Cm  DISTRIBUTION 
Mass  Averaged  Hub  to  Tip 


R/Rtip  -1.10 
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Propulsion  Division 


BASELINE  IMPELLER,  VISCOUS  SOLUTION 


IHJ1SS  flo'7  S|  lit  l!3c'o  • 


Utip  1858  fps  , Mi  i,d  Deg 


5b'\ 


ORIGINAL  PAGE  IS 
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Propulsion  Division 


a/c=0  14  h/s-0  067,  VISCOUS  SOI  U I ION 


iii.i-.s  tlo  v ‘‘|  'll*  - 47%  SVJ 


Utip=1858  fps  . Deg 


Bencdrp 

■dERCJJ^r  BASELINE  vs.  a/c=0.14  ; 1 Deg.  CLOCKED 

Propulsion  Division 
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□enCorp 

AEROJET 

Propulsion  Division 


Cm  DISTRIBUTION  JUST  DOWNSTREAM  OF  SLOT  LOCATION 

R/RTIP  = 0.94 


( BASELINE  VISCOUS) 


( a/oO.14  ; h/s-0.067  ) 


BenCdrp 

AEROJET 


ho 

00 


% BASELINE  HEAD  COEFF.  vs  SLOT  LOCATION 


MASS  FLOW  SPLIT  vs.  SLOT  LOCATION 


partial  ps  — partial  ss 


Propulsion  Division 


% BASEUNE  HEAD  COEFF.  vs  CLOCKING 


MASS  FLOW  SPLIT  vs.  CLOCKING 
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DenCdrp 

AEROJET 


CONCLUSIONS 


Propulsion  Division 


- FLOW  NON-UNIFORMITY  FOR  THE  BASELINE  CASE  WAS  FOUND  TO  BE  SIGNIFICANT 


- SIGNIFICANT  GAINS  IN  UNIFORMITY  CAN  POTENTIALLY  BE  ACHIEVED  WITH  TANDEM  BLADE 
CONCEPT 


- CLOCKING  WAS  SHOWN  TO  BE  THE  PRIMARY  DRIVER  FOR  IMPELLER  DISCHARGE  FLOW 
UNIFORMITY 


- CARE  MUST  BE  TAKEN  IN  THE  DESIGN  PROCESS  TO  ACHIEVE  POSITIVE  PERFORMANCE 
GAINS  RATHER  THAN  LOSSES 


- GRID  EMBEDDING  CAN  SIGNIFICANTLY  REDUCE  MODEL  SIZE  WHILE  NOT  SACRIFICING  FLOW 
GRADIENT  RESOLUTION 
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THREE-DIMENSIONAL  FLOW  FIELDS  INSIDE 
A SHROUDED  INDUCER  AT  DESIGN  AND 
OFF-DESIGN  CONDITIONS  ( CFD  STUDY ) 


C.  HAH,  O.  KWON,  AND  D.  A.  GREENWALD 
NASA  LEWIS  RESEARCH  CENTER 

R.  GARCIA 

NASA  MARSHALL  SPACE  FLIGHT  CENTER 


Three-dimensional  flow  phenomena  in  a shrouded  inducer  have  been  studied  with 
a three-dimensional  Navier-Stokes  method. 

The  details  of  the  three-dimensional  flow  structure  inside  the  inducer  at  design 
and  off-design  conditions  are  analyzed  and  the  results  are  compared  with  some  flow 
visualization  results  obtained  at  the  California  Institute  of  Technology. 
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THREE-DIMENSIONAL  FLOW  FIELDS  INSIDE 
A SHROUDED  INDUCER  AT  DESIGN  AND 
OFF-DESIGN  CONDITIONS  ( CFD  STUDY  ) 


C.  HAH,  O.  KWON,  AND  D.  A.  GREENWALD 
NASA  LEWIS  RESEARCH  CENTER 

R.  GARCIA 

NASA  MARSHALL  SPACE  FLIGHT  CENTER 


OBJECTIVES 


3-D  FLOW  STRUCTURE  AT  LOW  FLOW  COEFFICIENT 

FORMATION  OF  BACKFLOWS 

LATERAL  FORCES  IN  AXIAL  FLOW  INDUCERS 
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OVERALL  APPROACHES 


EXPERIMENTAL  INVESTIGATION 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 

PH.D  THESIS  BY  BH ATTACH ARYY A (PROF.  ACOSTA) 

ASME  PAPER  BY  BHATTACHARYYA,  ACOSTA,  BRENNEN 
AND  CAUGHEY  ON  " BACKFLOW  IN  INDUCER  " 

CFD  INVESTIGATION 

CURRENT  SUBJECT 


INDUCER  10 


SHROUDED  INDUCER  ( 3 BLADES  ) 

BLADE  ANGLE  = 12  DEGREES 

TIP  RADIUS  = 1.9115  INCHES 

TESTED  FLOW  COEFFICIENTS  ( 0.074  & 0.041 ) 


0.030 


N3 

vjO 
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wV3l- 


0.025 
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0.015 

0.010 

0.005 

0.000 
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LATERAL 


CALTECH  SHROUDED  INDUCER 
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COMPUTATIONAL  GRID  ( 40*31*122  ) 
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FLOW  COEFF.  = 0.074 


FLOW  ON  HUB  ( FROM  CALTECH  STUDY  ) 
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VELOCITY  VECTORS  NEAR  SHROUD  ( FLOW  COEFF.=0.074  ) 


300 


VELOCITY  VECTORS  NEAR  HUB  ( FLOW  COEFF.=0.074  ) 


V \VA  \ 


3 920219-16 


OIL  FLOW  ON  SHROUD  ( FLOW  COEFF.=0.074 ) 


3 320213-16:37:38 
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OIL  FLOW  ON  HUB  ( FLOW  COEFF.=0.074 ) 
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OIL  FLOW  ON  SUCTION  SIDE  ( FLOW  COEFF.=0.074 ) 
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OIL  FLOW  ON  PRESSURE  S[DE(  FLOW  COEFF.=0.074 ) 


FLOW  ON  SHROUD  (FROM  CALTECH  STUDY) 
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FLOW  COEFF.  = 0.041 


FLOW  ON  HUB  ( FROM  CALTECH  STUDY ) 
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ELOCITY  VECTORS  NEAR  SHROUD  ( FLOW  COEFF.=0.041 ) 

308 


3 92U22Q- 12:09:461 


VELOCITY  VECTORS  NEAR  HUB  ( FLOW  COEFF.=0.041 ) 
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3 920219-17 
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OIL  FLOW  ON  SHROUD  ( FLOW  COEFF.=0.041  I 
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OIL  FLOW  ON  HUB  ( FLOW  COEFF.=0.041 ) 
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OIL  FLOW  ON  SUCTION  SIDE  ( FLOW  COEFF.=0.041 ) 


920220-12:36: 16 
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OIL  FLOW  ON  PRESSURE  SIDE  ( FLOW  COEFF.=0.041 ) 
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OBSERVATIONS  FROM  CURRENT  EXERCISE 


MAJOR  PHENOMENA  WELL  CALCULATED 
FURTHER  QUANTITATIVE  COMPARISON  NEEDED 


CAVITATION  MODELING  NECESSARY 


N92-32292 


EFFECTS  OF  CURVATURE  AND  ROTATION  ON  TURBULENCE 

IN  THE 

NASA  LOW-SPEED  CENTRIFUGAL  COMPRESSOR  IMPELLER 


Joan  G.  Moore  and  John  Moore 


Mechanical  Engineering  Department 
Virginia  Polytechnic  Institute  and  State  University 

Blacksburg,  Virginia 


The  flow  in  the  NASA  Low-Speed  Impeller  is  affected  by  both 
curvature  and  rotation.  The  flow  curves  due  to 


a)  geometric  curvature 
nrnf  1 pq  TT.SI'idi 


, e.g.  the  curvature  of  the  hub  and  shroud 

onal  luju  uuc  uui  vdum'c  ui  Une 


backswept  impeller  blades,  and 


b)  secondary  flow  vortices,  e.g.  the  tip  leakage  vortex. 


Is  the  turbulence  and  effective  turbulent  viscosity  in  the 
impeller  significantly  affected  by  the  curvature  and  rotation  ? 

Do  these  changes  significantly  affect  the  overall 
three-dimensional  flow  development  ? 

And  do  they  also  impact  on  the  overall  performance  of  the 
impeller  ? 


An  answer  to  these  questions  is  obtained  by  comparing  two 
predictions  of  the  flow  in  the  impeller  - one  with,  and  one 
without  modification  to  the  turbulent  viscosity  due  to  rotation 
and  cu  r vatu  re . 


Some  experimental  and  theoretical  background  for  the  modified 
mixing  length  model  of  turbulent  viscosity  will  also  be  presented. 
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Prediction  of  Flow  in  the  Impeller 
using  a Mixing  Length  Model  for  Turbulent  Viscosity 

p = oL^"du" 
dy 

L = smaller  of  0.41ny" 

0.08  6 

Van  Driest  correction  used  in  0.41y  region 


Throughflow  velocity  vectors 
projected  onto  the  unrolled 
blade-to-blade  plane  at  80%  of 
blade  height. 
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Meridional  view  of  the  velocity 
vectors  at  mid-passage. 


A - 0.74. 


(a)  Near 
impeller 
inlet , 

A - 0.10. 


Cb) 

V = 0.40. 


(d)  At 

impeller 

exit, 

A « 1.0. 


Distributions  of  entropy  on  cross-sectional  planes  through  the  impeller.  Contour  interval  - 2.0  J/*g  K 
P,  pressure  side;  S,  suction  side. 


The  effects  of  tip  leakage  dominate  the  calculated  flow.  The 
calculation  was  made  with  a tip  gap  which  varied  from  2.6%  of 
blade  height  at  the  impeller  inlet  to  4%  at  the  impeller  exit.  The 
secondary  flow  velocity  vectors  in  the  cross-sectional  plane  are 
dominated  by  the  flow  over  the  blade  tip  and  the  resulting  vortex 
in  the  passage.  The  velocity  vectors  in  the  meridional  view  at  mid 
passage  show  the  extent  of  the  backflow  region  near  the  shroud  due 
to  the  tip  leakage.  The  vectors  in  the  blade-to-blade  view  show 
the  penetration  of  the  high  loss/low  velocity  tip  leakage  fluid  at 
80%  of  the  blade  height. 

The  entropy  on  four  cross-sectional  planes  show  the  build  up 
of  the  losses  to  be  dominated  by  the  tip  leakage  flow  with  the 
high  loss  fluid  covering  the  pressure-side/shroud  quarter  of  the 
passage  at  the  impeller  exit. 


Background  of  the  Modification 
of  the  Mixing  Length 
due  to  Curvature  and  Rotation 


MIXING  LENGTH  MODEL 


t 2 

p,  = pL  ou 
dy 


= Equilibrium  form  of  one  equation  k-L  model 


Turbulence  Kinetic  Energy  Equation 


- *.peffVk 


pu  • Vk 

Convection  Diffusion 
Equi 1 ibrium 


rk  ~ "'k 

Production  Dissipation 


= D, 


===> 


= pL2[  ( aUi  + au^  ) aUi  ]1/2 


dXj 


dxi 


ax . 

J 


- pL2"du“ 
dy 


The  mixing  length  turbulence  model  may  be  viewed  as  the 
equilibrium  form  of  a one-equation  turbulence  model.  It  can  be 
derived  by  setting  the  production  of  turbulence  kinetic  energy 
equal  to  its  dissipation. 

The  modification  to  the  turbulent  viscosity  due  to  curvature 
and  rotation  may  be  derived  by  considering  the  Reynolds  stress 
equations.  So  (1978),  derived  the  reduction  by  considering  the 
situation  when  rotation  and  curvature  act  in  a plane.  He  used  the 
equilibrium  form  (production  = dissipation)  of  the  three  normal 
stress  equations  and  one  shear  stress  equation. 
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MIXING  LENGTH  MODEL  WITH  CURVATURE  AND  ROTATION 


= P(LCF)2 


du 

dy 


= Equilibrium  form  of  6 equation  u?u’.  - L model 
— ^ ^ 

Reynolds  Stress  Equations 


pU-V(u^Uj)  - V- (K/a)effV(u^Uj)  = 


Convection 

Equilibrium 


Diffusion 


P. 


ij 


pu 


Di  0 


Production  Dissipation 


D.  . 
1J 


P (L  F)  “du" 
° dy 


F = F(Gradient  Richardson  Number) 


MODIFICATION  FACTOR  FROM  2-D  FLOW  LITERATURE 

F = 1 - B Ri  for  Ri  < 0 

F = 1/(1  + B Ri)  for  Ri  > 0 

Curvature  and  Rotation  Acting  in  a Plane  (So) 

Ri  = ( 2u/R  - 2fl)  (au/ay  + u/R  - 2D) 

( du/dy  - u/R  )^ 

(R  = radius  of  curvature) 
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Y/Delta 


Modified  Mixing  Length  from  Experimental  Data 
2-D  Flow  in  a Curved  Channel  (Gillis  et  al.) 
-pu’v’  = pt',du/dy‘,J 

2 

o - = pL  "du/dy" 


x ~ ^t  = CppLk 


1/2 


Gillis  et  al.  measured  the  six  Reynolds  stresses  and  the 
velocity  profile  on  the  suction  side  of  a curved  channel.  The 
local  turbulent  viscosity  and  the  local  mixing  length  for  either  a 
zero  equation  (Prandtl  mixing  length)  or  one-equation  (k-L) 
turbulence  model  may  then  be  determined  from  their  measurements. 
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L/Delta 


0.04 


0.°2  g0oo 


0.00  

0.0  -0.2 


°*°2  * 


0.00  -f^~r 
0 . 0 


Modified  L Model  using  local  Factor: 


L = L F 
o 


F = 1/(1  + BRi) 


Modified  L Model  using  mean  Factor: 

L = smaller  of  0.41"y" 

_ 6 
0.086  F = 0.08  / F dy 

0 


When  L,  determined  from  the  measurements,  is  compared  to  L 
from  the  model,  to  obtain  B,  it  is  found  that  using  a local 
modification  factor  results  in  the  wrong  shape  for  the  L versus  y 
profile  through  the  boundary  layer.  A mean  factor  applied  to  L in 
the  outer  part  of  the  boundary  layer  gives  better  results. 


PREDICTED  FLOW  IN  THE  IMPELLER 
using  a CURVATURE/ROTATION 
MODIFIED  MIXING  LENGTH  MODEL  FOR  3-D  FLOW 

P.  = pL2,,du" 
dy 

L = smaller  of  0.41"y" 

0.086  F 

Van  Driest  correction  used  in  0. 41y  region 
F = 1 - B Ri  for  Ri  < 0 

F = 1/(1  + B Ri)  for  Ri  > 0 
B = 4 


Ri  = [2£nkjukuiSuj/Sxi  - 2Enlk£kjiUlUinj1 
l£nla£aknulum(£kjiSuo/xi  " Znk>  1 
' [£nkjukui<sui/axj  + S"A’  ]2 


The  mixing  length  model,  with  L modified  by  a mean  factor  for 
curvature  and  rotation,  was  then  used  to  obtain  another  prediction 
of  the  flow  in  the  NASA  low  speed  centrifugal  impeller.  B=4  was 
used  with  a generalized  3-d  form  of  the  Richardson  number  which 
reduces  to  the  correct  2-d  form  in  2-d  situations. 
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—2 

Modification  factor  for  the  turbulent  viscosity,  F , 

3/4ths  of  the  way  through  the  impeller. 

The  classical  2-d  boundary  layer  modifications  can  be  seen 
with  enhancement  due  to  rotation  in  the  pressure  side  boundary 
layer,  reduction  due  to  rotation  in  the  suction  side  boundary 
layer,  and  enhancement  due  to  curvature  on  the  concave  hub  wall. 
In  the  tip  leakage  vortex  there  is  an  enhancement  due  to  rotation 
near  mid-passage  where  the  entropy  gradient  is  in  the  same  sense 
as  in  a pressure  side  boundary  layer.  The  tip  vortex  also  sees  a 
reduction  due  to  curvature  near  mid-height  where  the  entropy 
gradient  is  in  the  same  sense  as  for  a shroud  wall  (convex  wall) 
boundary  layer. 
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A = 0.97 


Secondary  flow  velocity  vectors  and  entropy  contours  on  cross- 
sectional  planes  3/4ths  through,  and  near  the  exit  of  the 
impeller.  Contour  interval  6s  = 2 J/kgK. 


The  character  of  the  flow  in  the  impeller  and  the  mean 
(mass-averaged)  levels  of  entropy  were  essentially  the  same  as  for 
the  first  prediction.  The  shape  and  location  of  the  tip  leakage 
vortex  were  slightly  modified  by  the  increased  turbulence  at 
mid-shroud  and  the  reduced  turbulence  in  the  tip  vortex  near  the 
pressure  side  near  mid-height. 
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Comparison  of  Rotor  Performance  Predictions  at  1920  RPM 
and  30  kg/s  (66  lbm/s)  with  Preliminary  Measurements 


NASA  Design 

Prediction 

Prediction 

NASA 

( Streamline 

Mixing  Length 

Curv/Rot 

Preliminary 

Curvature ) 

Turb.  Model 

Modified  L 

Measurements 

pt2/po 

1. 173 

1.  128 

1.  127 

1.  13 

n . . . _ 

0.934 

0.  922 

0.922 

0.  93 

AH/U2 

0.  607 

0.  469 

0.  467 

0.  50 

M .... 

0.31 

0.  31 

0.31 

reli,t.ip 

M . _ 

0.  20 

0.  272 

reiZ 

M . _ 

0.287 

0.  244 

absz 

Reaction 

0.763 

0.843 

0.846 

Plane  i is  at  the  impeller  inlet;  Plane  2 is  at  the  impeller  exit. 


TOTAL  PRESSURE  RISE 
IN  THE  IMPELLER 

o - NASA  Preliminary 
Measurements 

Predictions 

MEFP,  Mixing  Length 

Turbulence  Model 

- - - MEFP,  Rot/Curv 

Modified  L Model 

- - NASA  Design 

(Streamline  Curvature) 


Ptab/Po  - 1 


Both  MEFP  predictions  compare  well  with  the  preliminary 
measurements  made  by  NASA.  Comparisons  with  preliminary 
measurements  are  included  here  because  NASA  has  decided  to  reduce 
the  tip  clearance  before  making  detailed  measurements  of  flow  in 
the  impeller. 
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CONCLUSIONS 


o A local  modification  factor  due  to  curvature 
inappropriate  in  a mixing  length  model. 

Data  suggests  that  an  average  factor  for  the 
appropriate. 

o Typical  modification  factors  for  the  turbulent  viscosity  in  the 
NASA  Low  Speed  Centrifugal  Impeller  were  2-4. 

o Changes  to  the  predicted  flow  were  small: 

Slight  increase  in  secondary  flow  velocities 
on  the  pressure  side. 

Slight  change  in  shape  and  relocation 
of  the  tip  leakage  vortex. 

o No  change  in  overall  performance. 

o The  NASA  Low  Speed  Centrifugal  Impeller  is  a good  test  case 
for  verification  of  3-d  N-S  solvers  which  include 
mixing  length,  one-equation,  or  two-equation  turbulence  models. 

o The  Upwind  Control  Volume*  discretization  used  in  MEFP  makes 
this  code  a good  vehicle  for  the  evaluation  of  turbulence  models 
in  3-d  flow.  Calculated  vortex  structure  is  sensitive  to  small 
changes  in  the  turbulence  model. 


* The  Upwind  Control  Volume  approach  introduces  no  numerical 

mixing  either  directly  through  second  or  fourth  order  smoothing 
or  indirectly  through  inconsistencies  in  the  discretization  of 
the  convection  term  such  as  upwind  differencing. 


and  rotation  is 
layer  is 
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EFFECTS  OF  CURVATURE  AND  ROTATION  ON  TURBULENCE 


IN  THE 

NASA  LOW-SPEED  CENTRIFUGAL  COMPRESSOR  IMPELLER 


Joan  G.  Moore  and  John  Moore 


Mechanical  Engineering  Department 
Virginia  Polytechnic  Institute  and  State  University 

Blacksburg,  Virginia 


The  flow  in  the  NASA  Low-Speed  Impeller  is  affected  by  both 
curvature  and  rotation.  The  flow  curves  due  to 

a)  geometric  curvature,  e.g.  the  curvature  of  the  hub  and  shroud 

profiles  in  mpriHinnnl  plane  and  the  curvature  of  the 

backswept  impeller  blades,  and 

b)  secondary  flow  vortices,  e.g.  the  tip  leakage  vortex. 

Is  the  turbulence  and  effective  turbulent  viscosity  in  the 
impeller  significantly  affected  by  the  curvature  and  rotation  ? 

Do  these  changes  significantly  affect  the  overall 
three-dimensional  flow  development  ? 

And  do  they  also  impact  on  the  overall  performance  of  the 
impeller  ? 

An  answer  to  these  questions  is  obtained  by  comparing  two 
predictions  of  the  flow  in  the  impeller  - one  with,  and  one 
without  modification  to  the  turbulent  viscosity  due  to  rotation 
and  curvature. 

Some  experimental  and  theoretical  background  for  the  modified 
mixing  length  model  of  turbulent  viscosity  will  also  be  presented. 
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Prediction  of  Flow  in  the  Impeller 
using  a Mixing  Length  Model  for  Turbulent  Viscosity 


= pL  "diT 
* dy 


L = smaller  of  0.41"y" 


0.08  6 


Van  Driest  correction  used  in  0.41y  region 


* \ 


\\\ 

mwi 

rnw^v 


Calculated  secondary  flow 
velocity  vectors  on  cross- 
sectional  plane  three  quarters 
of  the  way  through  the 
impeller.  Note  contour  of  zero 
primary  velocity  near  the 
shroud  wall. 


Ms 

% %y 


W. 


Throughflow  velocity  vectors 
projected  onto  the  unrolled 
blade-to-blade  plane  at  80%  of 
blade  height. 


Meridional  view  of  the  velocity 
vectors  at  mid-passage. 


(c) 

A = 0,74, 


(d)  At 

inpeller 

exit, 

A - 1.0. 


Distributions  of  entropy  on  cross-sectional  plan 
P,  pressure  side;  S,  suction  side. 


es  through  the  impeller.  Contour  interval  * 2.0  J/fcg  K. 


The  effects  of  tip  leakage  dominate  the  calculated  flow.  The 
calculation  was  made  with  a tip  gap  which  varied  from  2.6%  of 
blade  height  at  the  impeller  inlet  to  4%  at  the  impeller  exit.  The 
secondary  flow  velocity  vectors  in  the  cross-sectional  plane  are 
dominated  by  the  flow  over  the  blade  tip  and  the  resulting  vortex 
in  the  passage.  The  velocity  vectors  in  the  meridional  view  at  mid 
passage  show  the  extent  of  the  backflow  region  near  the  shroud  due 
to  the  tip  leakage.  The  vectors  in  the  blade-to-blade  view  show 
the  penetration  of  the  high  loss/low  velocity  tip  leakage  fluid  at 
80%  of  the  blade  height. 

The  entropy  on  four  cross-sectional  planes  show  the  build  up 
of  the  losses  to  be  dominated  by  the  tip  leakage  flow  with  the 
high  loss  fluid  covering  the  pressure-side/shroud  quarter  of  the 
passage  at  the  impeller  exit. 


Background  of  the  Modification 
of  the  Mixing  Length 
due  to  Curvature  and  Rotation 


MIXING  LENGTH  MODEL 


pt 

= PL2"du" 
dy 

= Equilibrium  form  of 

one  equation  k 

-L  model 

Turbulence  Kinetic  Energy 

Equation 

pu*Vk  - 

Pk 

- Dk 

Convection  Diffusion 

Production 

Dissipation 

Equilibrium 

Pk 

= Dk 

===>  "t 

2 

= pL  [ ( + 

au . ) aui  ]1/2 

dx  . 
J 

axf 

pt 

= pL2"du“ 
dy 

The  mixing  length  turbulence  model  may  be  viewed  as  the 
equilibrium  form  of  a one-equation  turbulence  model.  It  can  be 
derived  by  setting  the  production  of  turbulence  kinetic  energy 
equal  to  its  dissipation. 

The  modification  to  the  turbulent  viscosity  due  to  curvature 
and  rotation  may  be  derived  by  considering  the  Reynolds  stress 
equations.  So  (1978),  derived  the  reduction  by  considering  the 
situation  when  rotation  and  curvature  act  in  a plane.  He  used  the 
equilibrium  form  (production  = dissipation)  of  the  three  normal 
stress  equations  and  one  shear  stress  equation. 
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MIXING  LENGTH  MODEL  WITH  CURVATURE  AND  ROTATION 


2 

Pt  = P(L0F) 

dy 


= Equilibrium  form  of  6 equation  u.u*.  - L model 
■■  — ...  3. J - 

Reynolds  Stress  Equations 


pu*7(uJUj)  - V- (|_./.j)effV(u^Uj)  = 
Convection  Diffusion 


Pi4 


D.  . 
ij 


Production  Dissipation 


Equilibrium 


P.  . - D.  . 

U 


P (L  F)2,,du' 
° dy 


F = F(Gradient  Richardson  Number) 


MODIFICATION  FACTOR  FROM  2-D  FLOW  LITERATURE 

F = 1 - B Ri  for  Ri  < 0 

F = 1/(1  + B Ri)  for  Ri  > 0 

Curvature  and  Rotation  Acting  in  a Plane  (So) 

Ri  = (2u/R  - 2fl)  (au/ay  + u/R  - 2D) 

( du/dy  - u/R  )2 

(R  = radius  of  curvature) 
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0.0  0.2  0.4  0.6'  0.8  1.0 

Y/Delta 


Modified  Mixing  Length  from  Experimental  Data 
2-D  Flow  in  a Curved  Channel  (Gillis  et  al. ) 
— pu 5 v ’ = Ht”du/dy‘', 

o - pt  = pL2 "du/dy" 

x - = C pLk1/2 

^ P 


Gillis  et  al.  measured  the  six  Reynolds  stresses  and  the 
velocity  profile  on  the  suction  side  of  a curved  channel.  The 
local  turbulent  viscosity  and  the  local  mixing  length  for  either  a 
zero  equation  (Prandtl  mixing  length)  or  one-equation  (k-L) 
turbulence  model  may  then  be  determined  from  their  measurements. 
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L/DeJta 


Modified  L Model  using  local  Factor: 

L = L F 
o 

F = 1/(1  + BRi) 


Modified  L Model  using  mean  Factor: 
L = smaller  of  0.41"y" 


6 

0. 086  F = 0.  08  / F dy 

0 


When  L,  determined  from  the  measurements,  is  compared  to  L 
from  the  model,  to  obtain  B,  it  is  found  that  using  a local 
modification  factor  results  in  the  wrong  shape  for  the  L versus  y 
profile  through  the  boundary  layer.  A mean  factor  applied  to  L in 
the  outer  part  of  the  boundary  layer  gives  better  results. 


PREDICTED  FLOW  IN  THE  IMPELLER 


using  a CURVATURE/ROTATION 
MODIFIED  MIXING  LENGTH  MODEL  FOR  3-D  FLOW 

p.  = pL^"du" 
dy 

L = smaller  of  0.41‘'y'‘ 

0. 086  F 

Van  Driest  correction  used  in  0.41y  region 
F = 1 - B Ri  for  Ri  < 0 

F = 1/(1  + B Ri)  for  Ri  > 0 
B = 4 


Ei  = t2£nkjukui8V8xi  - 2£nik£kjiuluinj1 
[£nla£aknulum<£kji8uj/Xi  ~ Znk'  1 
7 l£nkoukui(i>ui/8xj  + 8uj/8xi>  )2 


The  mixing  length  model,  with  L modified  by  a mean  factor  for 
curvature  and  rotation,  was  then  used  to  obtain  another  prediction 
of  the  flow  in  the  NASA  low  speed  centrifugal  impeller.  B=4  was 
used  with  a generalized  3-d  form  of  the  Richardson  number  which 
reduces  to  the  correct  2-d  form  in  2-d  situations. 
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Enhancement 
Curvature  only 


Reduction 
Curvature  only 


\\\ 


Net  Enhancement 
Factor 


m\\\ 


Net  Reduction 
Factor 


Modification  factor  for  the  turbulent  viscosity,  F , 

3/4ths  of  the  way  through  the  impeller. 

The  classical  2-d  boundary  layer  modifications  can  be  seen 
with  enhancement  due  to  rotation  in  the  pressure  side  boundary 
layer,  reduction  due  to  rotation  in  the  suction  side  boundary 
layer,  and  enhancement  due  to  curvature  on  the  concave  hub  wall. 
In  the  tip  leakage  vortex  there  is  an  enhancement  due  to  rotation 
near  mid-passage  where  the  entropy  gradient  is  in  the  same  sense 
as  in  a pressure  side  boundary  layer.  The  tip  vortex  also  sees  a 
reduction  due  to  curvature  near  mid-height  where  the  entropy 
gradient  is  in  the  same  sense  as  for  a shroud  wall  (convex  wall) 
boundary  layer. 


A = 0.97 


Secondary  flow  velocity  vectors  and  entropy  contours  on  cross- 
sectional  planes  3/4ths  through,  and  near  the  exit  of  the 
impeller.  Contour  interval  6s  = 2 J/kgK. 


The  character  of  the  flow  in  the  impeller  and  the  mean 
(mass-averaged)  levels  of  entropy  were  essentially  the  same  as  for 
the  first  prediction.  The  shape  and  location  of  the  tip  leakage 
vortex  were  slightly  modified  by  the  increased  turbulence  at 
mid-shroud  and  the  reduced  turbulence  in  the  tip  vortex  near  the 
pressure  side  near  mid-height. 
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Comparison  of  Rotor  Performance  Predictions  at  1920  RPM 
and  30  kg/s  (66  lbm/s)  with  Preliminary  Measurements 


NASA  Design 

Prediction 

Prediction 

NASA 

(Streamline 

Mixing  Length 

Curv/Rot 

Preliminary 

Curvature) 

Turb.  Model 

Modified  L 

Measurements 

pt2/po 

1.173 

1.  128 

1.  127 

1.  13 

nt-ti2 

0.  934 

0.  922 

0.  922 

0.93 

AH/U2 

0.  607 

0.  469 

0.467 

0.  50 

M , . . . 
reli, tip 

0.  31 

0.  31 

0.  31 

M 10 
rel2 

0.  20 

0.272 

M , 0 
abs2 

0.  287 

0.  244 

Reaction 

0.763 

0.843 

0.846 

Plane  i is  at  the  impeller  inlet;  Plane  2 is  at  the  impeller  exit. 


TOTAL  PRESSURE  RISE 
IN  THE  IMPELLER 

o - NASA  Preliminary 
Measurements 


MEFP,  Mixing  Length 

Turbulence  Model 

- - - MEFP,  Rot/Curv 

Modified  L Model 

- - NASA  Design 

( Streaml ine  Curvature ) 


Ptab/Po  - 1 


Both  MEFP  predictions  compare  well  with  the  preliminary 
measurements  made  by  NASA.  Comparisons  with  preliminary 
measurements  are  included  here  because  NASA  has  decided  to  reduce 
the  tip  clearance  before  making  detailed  measurements  of  flow  in 
the  impeller. 
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CONCLUSIONS 


o A local  modification  factor  due  to  curvature  and  rotation  is 
inappropriate  in  a mixing  length  model. 

Data  suggests  that  an  average  factor  for  the  layer  is 
appropriate. 

o Typical  modification  factors  for  the  turbulent  viscosity  in  the 
NASA  Low  Speed  Centrifugal  Impeller  were  2-4. 

o Changes  to  the  predicted  flow  were  small: 

Slight  increase  in  secondary  flow  velocities 
on  the  pressure  side. 

Slight  change  in  shape  and  relocation 
of  the  tip  leakage  vortex. 

o No  change  in  overall  performance. 

o The  NASA  Low  Speed  Centrifugal  Impeller  is  a good  test  case 
for  verification  of  3-d  N-S  solvers  which  include 
mixing  length,  one-equation,  or  two-equation  turbulence  models. 

o The  Upwind  Control  Volume*  discretization  used  in  MEFP  makes 
this  code  a good  vehicle  for  the  evaluation  of  turbulence  models 
in  3-d  flow.  Calculated  vortex  structure  is  sensitive  to  small 
changes  in  the  turbulence  model. 


* The  Upwind  Control  Volume  approach  introduces  no  numerical 

mixing  either  directly  through  second  or  fourth  order  smoothing 
or  indirectly  through  inconsistencies  in  the  discretization  of 
the  convection  term  such  as  upwind  differencing. 
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Computational  Fluid  Dynamic  Design  of  Rocket  Engine  Pump  Components 

by 

Wei-Chung  Chen,  George  H.  Prueger,  Daniel  C.  Chan  and  Anthony  H.  Eastland 
Rocketdyne  Division,  Rockwell  International  Corp. 


Integration  of  CFD  for  design  and  analysis  of  turbomachinery  components  is  needed 
as  the  requirements  of  pump  performance  and  reliability  become  more  stringent  for  the 
new  generation  of  rocket  engine.  A fast  grid  generator,  designed  specially  for  centrifugal 
pump  impeller,  which  allows  a turbomachinery  designer  to  use  CFD  to  optimize  the 
component  design  will  be  presented.  The  CFD  grid  is  directly  generated  from  the 
impeller  blade  G-H  blade  coordinates.  The  grid  points  are  first  generated  on  the 
meridional  plane  with  the  desired  clustering  near  the  end  walls.  This  is  followed  by  the 
marching  of  grid  points  from  the  pressure  side  of  one  blade  to  the  suction  side  of  a 
neighboring  blade.  This  fast  grid  generator  has  been  used  to  optimize  the  consortium 
pump  impeller  design.  A grid  dependency  study  has  been  conducted  for  the  consortium 
pump  impeller.  Two  different  grid  sizes,  one  with  10,000  grid  points  and  one  with 
80,000  grid  points  were  used  for  the  grid  dependency  study.  The  effects  of  grid 
resolution  on  the  turnaround  time,  including  the  grid  generation  and  completion  of  the 
CFD  analysis,  is  discussed.  The  impeller  overall  mass  average  performance  is  compared 
for  different  designs.  Optimum  design  is  achieved  through  systematic  change  of  the 
design  parameters.  In  conclusion,  it  is  demonstrated  that  CFD  can  be  effectively  used 
not  only  for  flow  analysis  but  also  for  design  and  optimization  of  turbomachinery 
components. 
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CFD  DESIGN  OF  ROCKET  ENGINE  COMPONENT 

by 

Wei-Chung  Chen,  George  H.  Prueger 
Daniel  C.  Chan,  Anthony  H.  Eastland 

Rocketdyne  Division 
Rockwell  International  Corporation 

Presented  at  NASA  Marshall  Space  Flight  Center 
Tenth  Workshop  for  Computational  Fluid  Dynamic  Applications 

in  Rocket  Propulsion 
April  28-30, 1 992 


Rockwell  International 

Rockstdyn*  Division 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

TYPICAL  TURBOPUMP  LAYOUT 


. KEY  COMPONENTS 

. INDUCER,  STATOR,  IMPELLER,  VOLUTE,  TURBINE,  BEARING 
. IMPELLER  SERVES  AS  THE  HEART  OF  PUMP 

. ADVANCED  IMPELLER  DESIGNS  REQUIRE  HIGH  HEAD  COEFFICIENTS 

. RELIABILITY  AND  COST  REQUIREMENTS  LIMIT  MAXIMUM  ALLOWABLE  TIP 
SPEED 

. LOW  COST  REQUIRES  MINIMUM  NUMBER  OF  STAGES 

. HIGH  PUMP  HEAD  COEFFICIENTS  INCREASE  FLOW  TURNING  AND  DIFFUSION 
. INCREASE  EXIT  FLOW  NON-UNIFORMITY 

. CFD  INCORPORATED  FOR  PUMP  IMPELLER  DESIGN 

. IDENTIFY  FLOW  PROBLEMS  INSIDE  IMPELLER  PASSAGE  AND  DISCHARGE 
. OPTIMIZE  IMPELLER  CONFIGURATION  DURING  DESIGN  PROCESS 


Rockwell  International 

Rockstdyns  Division 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

CONSTRUCTICN  OF  IMPELLER  BLADE 


• GENERATE  BLADE  MEANLINE  G-H  CURVE  ALONG  EACH  STREAMLINE 

• REDUCE  3-D  PROBLEM  TO  2-D  COORDINATE  SYSTEM 

. EACH  STREAMLINE  INDEPENDENTLY  GENERATED  TO  MATCH  FLOW  FIELD 
. CHANGE  BLADE  WRAP  ANGLE  TO  CONTROL  SOLIDITY 
. STACK  BLADE  TO  ACHIEVE  OPTI VIUM  IMPELLER  PERFORMANCE 

• CONSTRUCT  BLADE  SURFACE  FROM  MEANLINE  COORDINATES 

• SURFACE  COORDINATE  GENERATED  ACCORDING  TO  BLADE  THICKNESS 
DISTRIBUTION 

• FLEXIBILITY  OF  USING  MEANLINE:,  PRESSURE  SIDE,  SUCTION  SIDE  OR 
HYBRID  FAIRING 

. SURFACE  INFORMATION  DIRECTLY  USED  FOR 

• BLADE  LAYOUT  IN  CATIA  MODEL 

• CFD  GRID  GENERATION 


Rockwell  International 

Rockstdyns  Division 


IMPELLER  MEANLINE  BLADE  GENERATION 


CONSORTIUM  PUMP  IMPELLER 
STREAMLINE  DEFINITION  ( ITERATION  22  ) 


/,  INCH 


CONSORTIUM  PUMP 

IMPELLER  BLADE  G-H  CURVE  ( ITERATION  22  ) 


(;,  inch 


Rockwell  International 

Rockstdyn*  Division 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

FAST  GRID  GENERATOR 


• USE  BLADE  SURFACE  PRESSURE  AND  SUCTION  SURFACE  G-H  INFORMATION 
. BOTH  FULL  AND  PARTIAL  BLADES 

• GENERATE  GRID  IN  MERIDIONAL  PLANE 

. DETERMINE  L.E.  TO  T.E.  STATION  NUMBER 
. SELECT  HUB  TO  TIP  STREAMLINE  NUMBER 

• INTERPOLATE  BLADE  SURFACE  COORDINATES  TO  GRID  MESH  POINTS 
. REQUIRE  SURFACE  INTERPOLATION 

. CREATE  3-D  GRID  POINTS  BETWEEN  TWO  BLADE  SURFACES 

• EXTEND  GRID  POINTS  OUTWARDS 

. ACCORDING  TO  INLET  AND  OUTLET  BLADE  ANGLE  DISTRIBUTION 

• H-GRID  ALGEBRAICALLY  GENERATED 

• POISSON  GRID  SMOOTHER  INCORPORATED 

Rockwell  International 

Rockstdyns  Division 


9V£ 

RADIAL  DISTANCE 


DEVELOPMENT  OF  10K  IMPELLER  GRID 


COMPUTATIONAL  G1?1I) 


COM F5UTATION A L GRID 


MERIDIONAL  DISTANCE 


* 


Rockwell  International 

Rockstdyne  Division 


KAUIAL  DISTANCE 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

BOUNDARY  AND  INITIAL  CONDITIONS 

. BOUNDARY  CONDITIONS  AT  INLET 

. MERIDIONAL  AND  TANGENTIAL  VELOCITY  PRESCRIBED  FROM 
MEASUREMENT  OR  PREDICTION 

. BOUNDARY  CONDITION  AT  OUTLET 

. VELOCITY  EXTRAPOLATED  FROM  INTERIOR  POINTS 
. MASS  AND  ANGULAR  MOMENTUM  CONSERVED 

. ALONG  BLADE  SURFACES  AND  END  WALLS 
. NO-SLIP  BOUNDARY  CONDITION  IMPOSED 
. SLIP  END  WALL  EFFECTS  STUDIED 

. PERIODICITY  APPLIED  AT  INLET  AND  OUTLET  IN  BLADE-TO-BLADE  DIRECTION 
. INITIAL  CONDITIONS 

. UNIFORM  VELOCITY  ASSUMED  AT  EACH  STREAMWISE  STATION 
. VELOCITY  BASED  ON  1-D  PREDICTION 

. FLOW  DIRECTION  IS  ALIGNED  WITH  LOCAL  GRID  ORIENTATION 

* Rockwell  International 

Rocketdyne  Division 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

ACCURACY  OF  CFD  RESULTS 


. FLOW  SOLVER  REACT3D 

• ROCKETDYNE  ELLIPTIC  ANALYSIS  CODE  FOR  TURBOMACHINERY  USED 
. VALIDATED  FOR  INDUCER  AND  TURBINE  PERFORMANCE 

. COMPARE  TO  ONE-DIMENTIONAL  PROGRAM  FOR  IMPELLER  PERFORMANCE 
. IMPELLER  HYDRO  EFFICIENCY  AGREE  WITH  1-D  PROGRAM  ( 94.5%  ) 

. EULER  HEAD  8%  HIGHER  THAN  THAT  OF  1-D  PREDICTION 

. EULER  HEAD  DISCREPENCY  ATTRIBUTED  TO  DIFFERENCES  IN  DEVIATION 
ANGLE  AND  BLOCKAGE 

. 1-D  MODEL  EXTRAPOLATED  TO  HIGH  HEAD  COEFFICIENT 

• UNCERTAINTY  IN  SOLIDITY  CALCULATION  WITH  PARTIAL  BLADE 
. SIMPLIFIED  BLADE  TRAILING  EDGE  MODELING  IN  REACT3D 

. WATER  TEST  PLANNED  TO  RESOLVE  THESE  ISSUES 


Rockwell  International 

Rockstdyns  Division 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

DISCUSSION  OF  GRID  DEPENDENCY 


• 2 GRID  SIZE  USED  FOR  IMPELLER  ANALYSIS 

. 10K  (59X14X1 1 , STREAMWISE  X BLADE-TO-BLADE  X HUB-TO-TIP) 

. 80K  (119X30X23,  STREAMWISE  X BLADE-TO-BLADE  X HUB-TO-TIP) 

. GRID  GENERATED  ON  APOLLO  WORKSTATION 
. 2 CPU  MINUTES  FOR  10K  CASE 
. 10  CPU  MINUTES  FOR  80K  CASE 
. CFD  ANALYSIS  ALSO  CARRIED  OUT  ON  APOLLO  WORKSTATION 
. 4 CPU  HOURS  FORI  OK  CASE 
. 15  CPU  HOURS  FOR  80K  CASE 
. COMPARISON  OF  CFD  RESULTS  BETWEEN  10K  AND  80K 

. CONSISTENT  RESULTS  FOR  BOTH  BASELINE  AND  OPTIMUM  DESIGN 
. EFFICIENCY  AND  PUMP  HEAD  WITHIN  2% 

. TRENDS  IN  EVALUATION  CRITERIA  SIMILAR 
. USE  10K  FOR  IMPELLER  DESIGN  OPTIMIZATION 
. USE  80K  FOR  FINAL  OPTIMUM  DESIGN  FLOW  ANALYSIS 


Rockwell  International 

Rockstdyns  Division 
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COMPARISON  OF  CFD  SOLUTIONS 


CASE 

HEAD  (FT) 

EFFICIENCY 

FLOW  SPLIT* 

IMPELLER  DISCHARGE 
EC  1 EC  2 

CROSSOVER 
EC  1 

INLET 
EC  2 

1 0K 

BASELINE 

1256. G 

94.2% 

•<•/-  6% 

3.26  DEG 

0.142 

4.35  DEG 

0.0809 

BOK 

BASELINE 

1256.5 

94.1  % 

+/-  0.6% 

5.40  DEG 

0.148 

7.67  DEG 

0.0657 

1 0K 

OPTIMUM 

1250.0 

94.2% 

+ /-  2.5% 

2.00  DEG 

0.121 

2.99G  DEG 

0.0G03 

BOK 

OPTIMUM 

1272.9 

94.0% 

-»/-  5% 

2.36  DEG 

0.142 

5.516  DEG 

0.0536 

• HIGH  ON  FULL  BLADE  PRESSURE  SURFACE  10  PAR  HAL  BLADE  SUCTION  SURFACE 


Rockwell  International 

Rocketdyne  Division 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

DISCUSSION  OF  END  WALL  EFFECTS 


. 3 CASES  USED  FOR  PARAMETRIC  STUDY 

. CASE  1 : ROTATING  WALL  WITH  NO  SLIP  AT  IMPELLER  UPSTREAM  AND 
DOWNSTREAM 

• CASE  2:  SAME  AS  CASE  1 EXCEPT  SLIP  CONDITIONS  AT  DOWNSTREAM 

. CASE  3:  SAME  AS  CASE  1 EXCEPT  STATIONARY  SHROUD  AT  UPSTREAM 

. SMALL  CHANGE  FOR  ALL  PARAMETERS  EVALUATED 

• DIFFERENCE  OF  HEAD  AND  EFFICIENCY  WITHIN  1 % 

• SMALL  CHANGE  OF  IMPELLER  DISCHARGE  CM  AND  CU  DISTRIBUTION 

. FOURTH  CASE  WITH  STATIONARY  DOWNSTREAM  WALL  DID  NOT  FULLY 
CONVERGE 

. UNCONVERGED  RESULTS  ( RESIDUALS  10-2  ) SHOWED  EXTENSIVE 
RECIRCULATION  AT  DISCHARGE  HUB  AND  SHROUD 


Rockwell  International 

Rockstdyns  Division 


NON-DIMENSIONAL  CM  (CM/UTIP) 


0.12 


0.10  b 


0.08 


0.06  b 


0.04  b 


0.02  b 


0.00  >- 

2.5 


IMPELLER  EXIT  CM  DISTRIBUTION 
OPTIMIZED  GEOMETRY  WITH  VARIOUS 
BOUNDARY  CONDITIONS 


SHROUD 


-N. 


LEGEND 

Q ROT *1  INC  WAUS  W/NO  SIJI’  CONDITION  At  DISCHARGE 

0 ROtATIMG  WALLS  w/sur  CONDITION  AT  DISCHARGE 

V ROIAIIHO  WALLS  W/NO  SLID  CONDI  I ION  A1  DISCHARGE 
STATIONARY  SIIIIOUD  INI-ET  WALL 


.5.0 

X,  INCH 


Rockwell  International 
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vse 

NON-DIMENSIONAL  CL'  (CU/UTIP) 


IMI’KI.I.KK  EXIT  Hi  DISTRIBUTION 
(M’TIMIZKU  GEOMETRY  WITH  VARIOUS 
BOUNDARY  CONDITIONS 


LEGEND 


Q ROTATING  HALLS  W/NO  SLIP  CONDITION  AT  DISCHARGE 

£ ROTATING  WALLS  1/SLIP  CONDITION  AT  DISCHARGE 

V ROTATING  WALLS  W/NO  SUP  CONDITION  AT  DISCHARGE 
STATIONARY  SHROUD  INLET  WALL 


11U1J 


/ 


.>.0 

X,  INCH 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

DESIGN  PARAMETRIC  CHANGES 


. IMPELLER  DISCHARGE  HUB-TO-TIP  WIDTH 
. RANGE  FROM  0.88  TO  1 .20  INCHES 

. IMPELLER  DISCHARGE  BLADE  ANGLE  FROM  TANGENTIAL 
. RANGE  FROM  38  TO  54  DEGREES 

. IMPELLER  AXIAL  LENGTH  FROM  2.38  TO  3.0  INCHES 

. BLADE  TOTAL  WRAP  ANGLES 
. 65  TO  90  DEGREES  FOR  TIP 
. 65TO  115  DEGREES  FOR  HUB 

. WRAP  DIFFERENCE  BETWEEN  HUB  AND  TIP 
. RANGE  FROM  0 TO  5.0  DEGREES 

. ALL  CHANGES  ACHIEVE  REQUIRED  HEAD  WITH  CONSTANT  RPM  AND  DIAMETER 


Rockwell  International 

Rockstdyfls  04  vision 
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CONSORTIUM  IMPELLER  CONFIGURATION  STUDIED 


KEY  PARAMETER 

DISCHARGE 

DISCHARGE 

BLADE 

ANGLE 

AXIAL 

TOTAL  WRAP 

DISCHARGE 

WRAP 

DIFFERENCE 

CFD 

GRID 

POINT 

CASES 

B2  WIDTH 

LENGTH 

TIP(HUB) 

BASELINE 

1.00 

47.2 

2.50 

65  (85) 

5 

10000 

BASELINE 

1.00 

47.2 

2.50 

65  (85) 

5 

80000 

CASE  1 

1.00 

47.2 

2.50 

65  (85) 

0 

10000 

CASE  2 

0.88 

54.0 

2.38 

55  (78) 

3 

10000 

CASE  3 

1.12 

39.0 

2.62 

77  (100) 

2.7 

10000 

CASE  4 

1.20 

34.4 

2.70 

90  (115) 

5.0 

10000 

CASE  5 

1.00 

47.2 

3.00 

72  (92) 

5.0 

10000 

OPTIMUM 

1.12 

38.0 

2.82 

83  (105) 

2.6 

10000 

OPTIMUM 

1.12 

38.0 

2.82 

83  (105) 

2.6 

80000 

Rockwell  International 
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REACT  USED  TO  OPTIMIZE  HIGH  PERFORMANCE  IMPELLER 


BASELINE  GEOMETRY  00K  GRID 
IMPELLER  DISCHARGE  RADIAL  VELOCITY  (EPS) 


PRESSURE  ( psf ) 


FLOW  RECIRC  JLATION 
ELIMINATED  AT  IMPELLER 
DISCHARGE 


OPTIMUM  GEOMETRY  BOK  GRID 
PELLER  DISCHARGE  RADIAL  VELOCITY  (EPS) 


e 


Rockwell  International 

Rockstdyne  Division 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

CONCLUSIONS 


. A FAST  GRID  GENERATOR  HAS  BEEN  DEVELOPED 

. USED  FOR  DESIGN  OPTIMIZATION  OF  CONSORTIUM  IMPELLER 
. USED  FOR  SSME  HPFTP  IMPELLER  FOR  CFD  CODE  VALIDATION 
. SUCCESFULLY  APPLIED  FOR  INDUCER 

. CFD  INCORPORATED  INTO  PUMP  DESIGN  PROCESS 

. IMPELLER  PERFORMANCE  EVALUATION  CRITERIA  DEVELOPED 
. OPTIMUM  DESIGN  ACHIEVED  THROUGH  CFD  ANALYSIS 
. TURNAROUND  TIME  ACCEPTABLE  FOR  DESIGN  PROCESS 


Rockwell  International 
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CFD  DESIGN  OF  ROCKET  ENGINE  PUMP  COMPONENT 

CONCLUSIONS  (continued) 


. CFD  GRID  DEPENDENCY 

. SMALL  GRID  SIZE  ACCEPTABLE  FOR  PARAMETRIC  STUDIES 
. FINE  GRID  REQUIRED  ( 300K  ) TO  FINALIZE  GRID  DEPENDENCY  STUDY 

. END  WALL  BOUNDARY  EFFECTS 

. CHANGE  OF  IMPELLER  HEAD  AND  EFFICIENCY  WITHIN  1 % 

. SMALL  CHANGE  OF  IMPELLER  DISCHARGE  CM  AND  CU  DISTRIBUTION 

. TEST  DATA  REQUIRED 

. CONFIRM  HIGH  HEAD  COEFFICIENT  IMPELLER  PERFORMANCE 
. VALIDATE  CFD  RESULTS  FOR  CENTRIFUGAL  PUMP  APPLICATION 


Rockwell  International 

Rockstdyno  Division 
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SSME  HPOTP  IMPELLER  BACKCAVITY  CFD  ANALYSIS 

W.W.  HSU  and  S.J.  LIN 

Rockwell  International  Corp.,  Rocketdyne  Division 
6633  Canoga  Avenue,  MS  IA34 
Canoga  Park,  California  91303 


The  ball  bearings  behind  the  SSME  HPOTP  preburner  pump  have  a history  of 
premature  wear  requiring  their  replacement.  Extensive  tests  have  been 
conducted  in  an  attempt  to  identify  the  operating  factors  that  contribute  to  the 
wear.  It  has  been  conjectured  that  the  coolant  inflow  velocity  swirl  pattern  can 
aid  bearing  operation  by  matching  ball  orbit  speed  and  thus  affect  bearing  life. 
However,  control  of  the  velocity  distribution  up  to  now  could  only  be  achieved 
by  trial  and  error  following  hardware  testing.  Observation  of  hardware  from 
recent  flight  and  development  operation  led  to  the  hypothesis  that  certain 
assemblies  with  more  extensive  grinding  patterns  on  the  backwall  of  the 
impeller  for  rotor  balancing  correlated  with  improved  bearing  wear. 

To  analytically  evaluate  the  effect  of  cavity  configuration  on  the  flowfield,  3-D 
.CFD  analyses  of  various  geometries  was  successfully  executed  using 
REACT3D.  Height  of  the  anti-vortex  ribs  on  the  stationary  wall  was  varied,  as 
was  the  configuration  of  the  rotating  wall,  from  smooth  to  simulations  of  various 
grindout  patterns.  The  results  obtained  indicate  the  effects  of  the  various 
geometries  and  provide  valuable  guidelines  for  cavity  modification  to  optimi7P 
bearing  cooling. 
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SSME  HPOTP  PREBURNER  IMPELLER  BACKCAVITY 

CFD  ANALYSIS 


W.W.  HSU,  S.J.  LIN 

ROCKWELL  INERNATIONAL,  ROCKETDYNE  DIVISION 

APRIL  1992 


Rockwell  International 

Rocketdyne  Division 
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BACKGROUND 


SSME  HPOTP  BALL  BEARINGS  #1  AND  #2  BEHIND  PREBURNER 
IMPELLER  HAVE  HISTORY  OR  WEAR  AND  PREMATURE  REPLACEMENT 

EXTENSIVE  ENGINE  AND  SUBCOMPONENT  TESTS  HAVE  IDENTIFIED 
VARIOUS  OPERATING  FACTORS  THAT  AFFECT  WEAR 

LUBRICATION  CAGE  COATINGS  SUCH  AS  FEP,  BRAYCOTE 

AXIAL  PRE-LOAD  CORRECT  SPRING  STIFFNESS  TO 

MAINTAIN  PRE-LOAD 

MATERIALS  SILICON  NITRIDE  BALLS,  PLATINGS 

COOLANT  FLOW  MAINTAIN  ADEQUATE  VAPOR  MARGIN 

RECENT  TESTS  CORRELATED  MORE  EXTENSIVE  ROTOR  BALANCING 
GRINDOUTS  ON  IMPELLER  REAR  FACE  WITH  REDUCED  BEARING  WEAR 


Rockwell  International 

Rockstdyns  Division 


GEN  TUftSOPUMP 
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INCENTIVE  FOR  CFD  ANALYSIS 


COOLANT  INFLOW  SWIRL  DISTRIBUTION  CONJECTURED  TO  AFFECT 
BEARING  WEAR 


MATCHING  BALL  ORBIT  SPEED  REDUCES  INFLOW  RESISTANCE 
AND  DRAG  TORQUE  ON  BALLS  AND  CAGE 


SUGGESTS  DESIGN  CHANGES  TO  UPSTREAM  CAVITY  COULD  REDUCE 
BEARING  WEAR 


OPTIMIZE  MAGNITUDE  AND  RADIAL  DISTRIBUTION  OF  INLET 
SWIRL  VELOCITY 


ADJUST  HEIGHT  OF  ANTI- VORTEX  RIBS  AND/OR  SIZE  OF 
IMPELLER  GRINDOUTS  TO  ACHIEVE  DESIRED  DISTRIBUTION 


Rockwell  International 
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CFD  ANALYSIS  OBJECTIVES 


DEFINE  VARIATION  OF  INLET  SWIRL  VELOCITY  WITH  GEOMETRY 

EFFECT  OF  ANTI-VORTEX  RIBS  AND  IMPELLER  GRINDOUTS 
INDIVIDUALLY  AND  TOGETHER 

UNDERSTAND  FLOW  WELL  ENOUGH  TO  SUGGEST  DESIGN  CHANGES 

NARROW  CAVITY  WITH  HIGH  WALL  TANGENTIAL  VELOCITY 
RIBS  ON  STATIONARY  WALL 
HIGH  VELOCITY  JET  AT  INLET 


Rockwell  International 

Rocketdyne  Division 


CFD  MODELING 


REACT3D  STEADY  NAVIER-STOKES  ANALYSIS,  SINGLE  ZONE 

GRINDOUTS  ANALYZED  IN  ROTATING  FRAME  OF  REFERENCE 

RIBS  ANALYZED  IN  STATIONARY  FRAME  OF  REFERENCE 

QUASI-STEADY  APPROACH  PROPOSED  FOR  GRINDOUT-RIB 
COMBINATIONS 

FLOWFIELD  SIMPLIFICATIONS 

SMALL  LEAKAGE  PATH  PARALLEL  TO  BEARINGS  IGNORED 

EFFECTS  OF  ROTATING  CAGE  AND  ROLLING  BALLS  NOT 
SIMULATED 

GRID  EXTENDED  DOWNSTREAM  BEYOND  BEARING  INLET 

AIDS  CONVERGENCE  WITH  POTENTIAL  BACKFLOW 


Rockwell  International 

Rocketdyne  Division 


GEOMETRIES  INVESTIGATED 


FOUR  BASIC  GEOMETRIES  TO  BE  INVESTIGATED  (3  COMPLETE) 

RIBBED  STATIONARY  WALL  / SMOOTH  ROTATING  WALL 

SMOOTH  STATIONARY  WALL  / SMOOTH  ROTATING  WALL 

SMOOTH  STATIONARY  WALL  / ROTATING  WALL  WITH  GRINDOUTS 

RIBBED  STATIONARY  WALL  / ROTATING  WALL  WITH  GRINDOUTS 

GRINDOUTS  SIMULATED  WITH  SMOOTH  INDENTATIONS  ON  ROTATING 
WALL 

RADIAL  HEIGHT  MATCHED  TO  AVERAGE  OBSERVED 
TOTAL  CIRCUMFERENTIAL  EXTENT  180  DEG. 

GEOMETRY  VARIATIONS 

RIB  HEIGHTS  100%,  50%,  25%  OF  CURRENT  DESIGN 
GRINDOUTS  WITH  4 LOBES  0.05"  DEEP,  6 LOBES  0.1"  DEEP 


Rockwell  International 

Rocketdyns  Division 


CFD  ANALYSIS  CONDITIONS 


GRID  SIZE  MERIDIONAL  DIRECTION  I = 75  TO  85 

CIRCUMFERENTIAL  DIRECTION  J =15  TO  30 
NORMAL  DIRECTION  K = 1 5 TO  20 

GEOMETRY  IMPELLER  ROTATION  29141  RPM 

HUB  SEAL  RADIUS  1 .67  INCH 
SEAL  GAP  0.005  INCH 

HUB  INNER  RADIUS  1.06  INCH 
CAVITY  AXIAL  WIDTH  0.113  INCH 

RIB  HEIGHT  0.05  INCH 

FLUID  LOX,  DENSITY  57.4  LB/CUB.  FT. 

FLOWRATE  10.8  LB/SEC 

JET  VELOCITY  AXIAL  245  FT/SEC  (SEAL  EXIT  CHAMFERED) 

TANGENTIAL  50%  WHEEL  SPEED  (PHASE  1 ) 

30%  WHEEL  SPEED  (PHASE  2) 
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RESULTS 


REACT3D  RESULTS  QUALITATIVELY  AS  EXPECTED 

GRINDOUTS  SIGNIFIICANTLY  INCREASE  INLET  TANGENTIAL 
VELOCITY 

RIBS  SIGNIFICANTLY  DECREASE  INLET  TANGENTIAL  VELOCITY 

STRONG  VORTEX  MOTIONS  DRIVEN  BY  JET  FROM  DAMPING 
SEAL  AND  ROTATING  WALLS 

NO  SIGNIFICANT  IMPACT  ON  AXIAL  VELOCITY 

FLOW  SEPARATIONS  OFF  INNER  RACE  MAY  SIGNIFICANTLY 
REDUCE  THROUGHFLOW 

LESS  SIGNIFICANT  SEPARATION  AT  BOTTOM  OF  RIBS 

RIB  DEPTH  CAN  BE  USED  TO  CONTROL  INLET  VELOCITY 

JET  INLET  TANGENTIAL  VELOCITY  HAS  LITTLE  EFFECT  ON  BEARING 
INLET  VELOCITY  PROFILE 
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Smooth  Kolalin^  Wall  with  Rib 

Tangential  Velocity  (ft/s) 
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Smooth  Rotating  Wall  without  Rib 

Tangential  Velocity  (ft/s) 
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(•round  Rotating  Wall  without  Kill 

Tangential  Velocity  (ft/s) 
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Smooth  Rotating  Wall  with  Rih 

Axial  Velocity  (t't/s) 
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Smooth  Rotating  Wall  without  Rib 

Axial  Velocity  (li/s) 
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( i round  Rolalin^  Wall  without  Rib 

Axial  Velocity  (It/s) 
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Smooth  Rotating  Wall  with  Rih 
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Ground  Rotating  Wall  without  Rib 
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RADIUS  - INCH 


HPOTP  PREBURNER  BACKCAVITY  (AT  INNER-RACE  INLET) 

U = 50%  ROTATIONAL  SPEED  U,r  = 30%  ROTATIONAL  SPEED 

a W/  RIB,  FULL  RIB  HEIGHT  a W/  RIB,  FULL  RIB  HEIGHT 

O W /RIB,  HALF  RIB  HEIGHT 
v W/  RIB,  QUARTER  RIB  HEIGHT 
O W/O  RIB, 

• W,/0  RIB,  HUB  GRIND  4 LOBES/0.05  DEEP 
■ W/O  RIB,  HUB  GRIND  6 LOBES/O.IO  DEEP 


50  -100  -50  0 50  100 


TANGENTIAL  VELOCITY,  DIFFERENTIAL  ( Cu  0 43  u)  FPS  | 04/23/92 
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CONCLUSIONS  / FUTURE  WORK 


REACT3D  CFD  RESULTS  PROVIDE  INSIGHT  INTO  COMPLEX  FLOW 
PHYSICS 


ADVANCED  POST-PROCESSING  ESSENTIAL  TO 
UNDERSTANDING  FLOW 


RESULTS  CAN  BE  USED  TO  DIRECT  REDESIGN  EFFORTS 


CASE  COMBINING  RIBS  AND  GRINDOUTS  IN  WORK 


Rockwell  International 

Rocketdyne  Division 
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NLS  CLUTCHING  BEARING  CAVITY  FLOW  ANALYSIS 

Ken  Tran,  Daniel  C.  Chan  and  Armen  Darian 
Rocketdyne  Division  - Rockwell  International 
6633  Canoga  Ave.,  Canoga  Park  - CA  91303 

In  turbopumps  with  hydrostatic  bearings,  clutching  bearings  are  one 
technique  that  can  be  used  to  control  the  transient  axial  thrust.  At  steady 
state  operation,  the  clutching  bearing  inner  race  is  decoupled  from  the 
rotating  shaft  and  spins  at  a speed  which  is  determined  by  the  fluid  driving 
forces  in  the  bearing  cavity  and  the  ball  bearing  resistance.  The  life  of  the 
clutching  bearing  depends  on  the  speed  of  rotation  of  the  inner  race; 
therefore,  it  is  important  to  predict  the  latter  with  accuracy. 

Cavity  flow  analysis  is  difficult  due  the  complicated  nature  of  the 
geometry,  which  often  results  in  a totally  skewed  mesh.  A quick  study  of  a 
simple  cavity  flow  was  performed  to  gain  insight  into  important  parameters. 
It  was  concluded  that  the  multi-domain  ( or  multi-zone)  approach,  the  double 
precision  code,  the  initial  condition  and  a good  combination  of  relaxation 
factors  are  the  4 essential  features  in  the  search  for  a quick  converged 
solution.  The  multi-domain  approach  enables  the  user  to  divide  the  model 
into  small  blocks  which  are  gridded  separately;  therefore  insuring  the 
creation  of  a reasonable  mesh.  The  double  precision  code  solves  the  problem 
of  various  scales  in  different  regions  of  the  flow  and  the  good  initial  guess  in 
conjunction  with  a good  selection  of  relaxation  factors  helps  reduce  the 
computational  time. 

A flow  model  of  the  NLS  clutching  bearing  cavity  was  built  for  2-D 
axisymmetric  viscous  analyses.  From  the  CFD  output,  the  tangential  force 
exerted  on  the  surfaces  of  the  inner  race  was  integrated  to  calculate  the 
driving  torque  which,  in  conjunction  with  the  resistance  torque,  was  used  to 
predict  the  operating  speed  of  the  inner  race. 

In  order  to  further  reduce  the  inner  race  rotation,  the  swirling  flow  at  the 
cavity  inlet  was  partially  re-directed  to  generate  an  opposing  torque.  Thirty 
six  slanted  slots  was  incorporated  into  the  anti-vortex  rib  to  achieve  this 
goal.  A 3-D  flow  analysis  performed  on  this  configuration  indicated  a drastic 
reduction  of  the  driving  torque  and  inner  race  RPM. 
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NLS  CLUTCHING  BEARING  FLOW 

ANALYSIS 


By  Ken  Tran,  Daniel  C.  Chan  and  Armen  Darian 
Rocketdyne  Division,  Rockwell  International 


Rockwell  International 

Rocketdyne  Division 
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CONTENT 


• Clutching  bearing  description 

• Objective  and  Methodology 

• Two-Dimensional  analysis 

• Three-Dimensional  analysis 

• Concluding  remarks 


Rockwell  International 

Rocketdyne  Division 


392 


CLUTCHING  BEARING  FUNCTION 


• During  transient  operation: 

-Bearing  element  is  used  to  control  transient  axial  thrust 

• At  steady  state  operation: 

- Bearing  is  decoupled  from  the  rotating  shaft 

- Balance  piston  takes  over  the  control  of  the  axial  thrust 

- Inner  race  is  induced  to  rotate  by  fluid  friction 


Rockwell  International 

Rocketdyne  Division 
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NLS  CLUTCHING  BEARING  CONFIGURATION 


Rockwell  International 

Rocketdyne  Division 


395 


OBJECTIVE 


. Determine  life  of  the  clutching  bearing 

• Predict  the  inner  race  RPM 

. Calculate  the  torque  acting  on  the  bearing  inner  race 
faces 

• Investigate  design  features  to  minimize  inner  race 
speed 


Rockwell  International 

Rockotdyno  Division 
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APPROACH 


• Estimate  resistance  torque  as  a function  of  inner  race  RPM 
from  bearing  mechanical  characteristics 

• Determine  driving  torque  in  function  of  inner  race  RPM 

• Intersection  of  2 torque  curves  determines  inner  race  speed 

• Simplify  the  flow  geometry 

- Labys  are  not  modeled 

- Flow  through  the  balls  not  included 


Rockwell  International 

Rocketdyne  Division 


CFD  METHODOLOGY 


• Reynolds  averaged  Navier-Stokes  solver  (REACT): 

- Control  volume,  pressure  correction  method 

- Two-equation  k-s  turbulence  model 

• Validation: 

- Daily  and  Nece  cavity 


Rockwell  International 
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COMPUTATIONAL  MODEL 


• Following  elements  are  important  to  achieve  converged 
solution: 

- Multi-domain  grid:  better  mesh  and  control  of  Y+ 

- Double  precision  code 

- Good  initial  condition 

- Relaxation  factors:  Taguchi  parametrics  for  simple  case 
performed  on  u,v,w  (0.35 , 0.5,  0.8)and  p (0.1,  0.15, 
0.2)  relaxation  factor 


Rockwell  International 
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COMPUTATIONAL  MODEL  (cont.) 

• Grid  resolution: 

- Y+  ~ 50-500 

- 4205  grid  points  for  2-D  model 

- 95366  grid  points  for  3-D  model:  meridional  grid 
identical  to  2-D  mesh 

• Boundary  conditions: 

- Swirling  jet  imposed  at  one  circumferencial  node  line 

The  jet  represents  the  exit  condition  of  the 
hydrostatic  bearing 


Rockwell  International 
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COMPUTATIONAL  MESH 


SINGLE  ZONE  MULTI-ZONE 
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TAGUCHI  ANALYSIS  OF  RELAXATION 

FACTORS 


2 1 1 

1 st  2nd  3rd 

LEVEL 
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COMPUTATIONAL  GRID 
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2-D  CFD  RESULTS 


• Streamlines  show  the  jet  diffusion  is  slow 

• Swirl  flow  is  still  present  near  the  inner  race  front  face 

• Radial  pressure  distribution  on  the  inner  race  front  face  is 
relatively  uniform  except  at  stagnation  point: 

-1-D  model  can  be  used  to  estimate  axial  load 

• Predicted  axial  load  is  independent  of  inner  race  speed 
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2-D  CFD  RESULTS  (cont.) 


• Integrated  torque  acting  on  the  inner  race  indicates  that: 

- At  zero  RPM,  the  inner  race  front  face  contributes  over 
50%  to  the  total  torque 

- At  10000.  RPM,  this  contribution  is  only  25% 

• Resistance  and  driving  torque  characteristics  determine 
inner  race  RPM: 

- Predicted  RPM  is  about  9000 


Rockwell  International 
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MERIDIONAL  VELOCITY 
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TORQUE  VERSUS  INNER  RACE  RPM 


4> 

O 
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PRESSURE  (PSI) 


PRESSURE  DISTRIBUTION  ON  INNER  RACE 

FRONT  FACE 


RADIUS  (IN.) 


Rockwell  International 

Rocketdyne  Division 


409 


METHODS  TO  REDUCE  INNER  RACE 

SPEED 


• Increase  axial  load  to  raise  resistance  torque: 

-Tighter  laby  clearance:  unacceptable  due  to  high 
assembly  cost 

• Reduce  the  effect  of  the  swirl: 

- Anti- vortex  ribs:  limited  result  because  of  small 
contribution  of  the  front  face  torque 

- Redirecting  the  jet  against  direction  of  rotation  to  offset 
driving  torque 

• Rockwell  International 

Rocketdyn*  Division 
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3-D  CFD  ANALYSIS 


. Feature:  36  Slots  on  rib  used  to  redirect  jet  ( 20  deg.  from 
axial) 

• Results: 

- Driving  torque  reduced  significantly 

- Inner  race  speed  lowered  to  7000  RPM 

- 12  Hrs  of  YMP  Cray  CPU  time  per  case 

• Slots  can  be  modified  to  increase  effectiveness 


Rockwell  International 
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TORQUE  VERSUS  INNER  RACE  SPEED 
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CONCLUDING  REMARKS 

• Multi-domain  methodology  is  needed  for  complicated 
geometry 

• REACT  is  a mature  code: 

- Enables  the  exploration  of  new  design  concepts 

- Helps  determine  optimum  configurations 

• Graphics  post-processing  is  an  important  tool: 

- Assists  users  in  the  understanding  of  the  flow  field 

- Highlights  flow  features 


Rockwell  International 
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CFD  Analysis  to  Optimize 
a Design  Modification  of  BSMT 

Mark  Ratcliff  and  Ram  Awa 
CFD  Research  Corporation,  Huntsville,  Alabama 
Robert  Williams 

NASA  MSFC,  Huntsville,  Alabama 

The  Bearings,  Seals  and  Material  Tester  (BSMT)  is  a test  article  being  used  at  MSFC  to 
evaluate  the  performance  of  conventional  rolling  contact  bearings.  Pressure  differentials 
between  the  BSMT  inlet  and  exit  cavities  are  found  to  cause  large  parasitic  axial  loads  on 
the  bearing-carrier  walls.  These  parasitic  loads,  besides  being  detrimental  to  the  life  of 
bearings,  make  the  testing  and  evaluation  of  bearing  performance  very  difficult,  and  need 
to  be  eliminated  if  at  all  possible. 

CFDRC  is  currently  under  contract  to  MSFC  to  perform  a detailed  analysis  of  the  flow  fields 
inside  the  BSMT  cavities  and  manifolds.  The  objectives  of  this  study  ate  to  estimate  the 
hydrodynamic  loads  on  the  bearings  and  to  recommend  feasible  design  modifications  for 
BSMT  to  eliminate  the  parasitic  loads. 

Three-dimensional  computational  analyses  of  inlet  and  exit  cavities  in  their  baseline 
configuration  were  performed  with  REFLEQS  which  is  an  advanced  finite-volume  Navier- 
Stokes  code.  Computations  were  performed  with  and  without  a 1/4  inch  diameter 
temperature  probe  included  in  each  of  the  cavities.  The  results  of  the  analyses  indicate  that 
the  'temperature  probes  substantially  alter  the  flow  field  and  reduce  the  pressure  drop/rise 
in  the  cavities.  The  overall  pressure  drop  across  the  tester  compares  quite  well  with  the 
measurements. 

One  of  the  potential  design  modifications  to  reduce  the  parasitic  loads  on  the  bearings  is  to 
place  baffles  in  the  inlet  cavities  to  isolate  the  coolant  flow  from  the  slinger  wall.  Three- 
dimensional  analyses  of  the  inlet  cavities  with  the  baffle  were  performed  to  assess  the  effect 
of  baffles  on  the  axial  load.  The  baffle  length  was  varied  as  a parameter.  Results  suggest 
that  axial  loading  should  be  reduced  considerably  with  the  baffle  extended  inward  to  the 
radius  of  the  outer  race. 

Thermal  analyses  of  the  inlet  cavities  were  performed  to  determine  the  temperature  rise  due 
to  viscous  dissipation.  The  deflection  of  the  baffle  due  to  the  hydrodynamic  pressure  load 
was  also  determined  by  performing  structural  analysis.  The  analyses  suggest  that  the 
temperature  rise  and  the  baffle  deflection  are  not  of  much  concern.  Therefore  the 
considered  design  modification  seems  feasible  and  should  be  investigated  further  from 
structural,  manufacturing,  and  test  assembly  considerations. 
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OUTLINE 


CFMC 


• Overview 

• Modeling  of  Inlet  and  Exit  Cavities 

• Results  of  3D  Analyses 

• Summary  and  Concluding  Remarks 
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Bearings,  Seals  & Material  Tester 


CFMC 


Exit  Ports 


Radial  Loader 


BSMT  INSTRUMENTATION 


Upstream 

Temperature  Probe 
(Bearing  Inlet) 


Downstream 
Temperature  Probe 
(Bearing  Exit) 
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STATEMENT  OF  PROBLEM 


CFMC 


• Problem 

- Tester  Designed  with  Improper  Cavity  Pressure 
Distribution 

- Pressure  Differentials  Between  BSMT  Inlet  and 
Exit  Cavities  Cause  Large  Parasitic  Axial  Loads 


• Objective  of  Design  Modification 
- Reduce  Axial  Loads  in  BSMT 
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NUMERICAL  CODE 


REactive  FLow  EQation  Solver  (REFLEQS) 


CFD^C 


• Density  Averaged  N-S  Equations 

• Finite  Volume 

• Pressure-Based  Algorithm  (SIMPLEC) 

• Incompressible  and  Compressible  Flows 

• Cartesian,  Axisymmetric  and  BFC  Options 

• Turbulence  and  Combustion  Models 


MODELING  OF  BSMT  CAVITIES 

LOX  at  -279°F  (100°K)  and  600  PSI  (4x10  6 Pa) 


CrMC 


Assumptions: 

• Incompressible 

• Isothermal 

• Constant  Properties 

• Single  Phase 


427 


MODELING  OF  INLET  CAVITIES 
Mass  Flux  Matched 


CFDRC 


Rotating 


INLET  CAVITY  NO.  1 


INLET  CAVITY  NO.  2 


» t 


MODELED  GEOMETRY 
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MODELING  OF  EXIT  CAVITY 


CFDRC 


Actual  Geometry 


Modeled  Geometry 
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BOUNDARY  CONDITIONS 


CFDRC 


Specified  Mass  Flow 


INLET  CAVITY  NO.  1 INLET  CAVITY  NO.  2 
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CAVITY  NO.  1 
COMPUTATIONAL  GRID 
30  x 30  x 34 


r0- PLANE 


Xr- PLANE 
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INLET  & EXIT  LOCATIONS 


CFD=rc 


Tester  Inlet 
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SUMMARY  OF  3D  ANALYSES 
Baseline  with  Probes 


CFDRC 


n 

Case 

r> 

Static  Pressure  at  (psi) 

Net  A P across 
Tester  (psi) 

Bearing 
Inlet  Swirl 
% cage 

Net  Load 
Ibf 

Tester 

Inlet 

Bearing 

Inlet 

Bearing 

Exit 

Tester 

Exit 

Computed 

Measured 

speed 

Cavity  No.  1 
with  probe 

600 

570 

555 

578 

22 

24 

42 

945 

Cavity  No.  2 
with  probe 

600 

556 

541 

564 

36 

38 

60 

1,495 
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POTENTIAL  DESIGN  MODIFICATION 


CFD=TC 


Problem: 


Cause: 


Solution: 


Large  axial  (parasitic)  loads 


CAVITY  NO.  1 (66%  BAFFLE) 


STATIC  PRESSURE  ON  CARRIER  WALL 


CAVITY  NO.  1 (100%  BAFFLE) 
STATIC  PRESSURE  ON  CARRIER  WALL 
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SUMMARY  OF  3D  ANALYSES 


CFMC 


Case 

Static  Pressure  at  (psi) 

Net  AP 
across 
Tester  (psi) 

Bearing 
Inlet  Swirl 
% cage 

Net 
Axial 
Load  on 

Tester 

Inlet 

Bearing 

Inlet 

Bearing 

Exit 

Tester 

Exit 

Measured 

speed 

Bearings 

Ibf 

Drive-side 

Baseline 

600 

570 

555 

578 

22 

24 

42 

945 

66%  Baffli 

0 

m 

585 

570 

593 

7 

27 

345 

100%  Baffle 

600 

591 

576 

599 

1 

21 

105 

1 

Baseline 

o 

556 

541 

564 

36 

38 

60 

1,495 

66%  Baffle 

600 

586 

571 

594 

6 

41 

265 

1 00  % Baf  fle 

HI 

591 

576 

599 

1 

38 

75 
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BSMT  THERMAL  ANALYSIS 


CFMC 


1 st  Law  of  Thermodynamics 


Assuming  adiabatic  walls 
w = rh[(h  + KE)out  - (h  + KE)jn] 

h - h. 

Bulk  AT  = — 111 

cp 

Solving  Enthalpy  Equation  # Peak  AT 


Inlet  Cavity  No.  1 (Drive-Side) 


Baseline 

24.9 

6.72 

17.5 

33%  baffle 

27.5 

7.4 

66%  baffle 

25.5 

7.0 

24.3 

100%  baffle 

22.9 

6.5 

out 


STRUCTURAL  ANALYSIS 
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N92-32297 


Combustion  Instability  Analysis  For  Liquid  Propellant  Rocket  Engines 

\ .M.  Kim,  C.P.  Chen,  and  J.P.  Ziebarth 

University  of  Alabama  in  Huntsville 
Huntsville,  Alabama  35S99 


Abstract 

The  multi-dimensional  numerical  model  has  been  developed  to  analyze  the 
nonlinear  combustion  instabilities  in  liquid-fueled  engines.  The  present  pressure- 
based  approach  can  handle  the  implicit  pressure-velocity  coupling  in  a non-iterative 
way.  The  additional  scalar  conservation  equations  for  the  chemical  species,  the 
energy,  and  the  turbulent  transport  quantities  can  be  handled  by  the  same  predictor- 
corrector  sequences.  This  method  is  time-accurate  and  it  can  be  applicable  to  the 
all-speed,  transient,  multi-phase,  and  reacting  flows. 

Special  emphasis  is  given  to  the  acoustic/vaporization  interaction  which  may 
act  as  the  crucial  rate-controlling  mechanism  in  the  liquid-fueled  rocket  engines. 
The  subcritical  vaporization  is  modeled  to  account  for  the  effects  of  variable  ther- 
mophysical properties,  non-unitary  Lewis  number  in  the  gas-film,  the  Stefan  flow 
effect,  and  the  effect  of  transient  liquid  heating.  The  test  cases  include  the  one-  di- 
mensional fast  transient  non-reacting  and  reacting  flows,  and  the  multi-dimensional 
combustion  instabilities  encountered  in  the  liquid-fueled  rocket  thrust  chamber.  The 
present  numerical  model  successfully  demonstrated  the  capability  to  simulate  the 
fast  transient  sprav-combusting  flows  in  terms  of  the  limiting-cycle  amplitude  phe- 
nomena, correspondence  between  combustion  and  acoustics,  and  the  steep-fronted 
wave  cz  flame  propagation.  The  investigated  parameters  include  the  spray  initial 
conditions,  air-fuel  mixture  ratios,  and  the  engine  geometry.  Stable  and  unstable 
operating  conditions  are  found  for  the  liquid-fueled  combustors.  Under  certain  con- 
ditions, the  limiting  cycle  behavior  of  the  combusting  flowfields  is  obtained.  The 
numerical  results  indicate  that  the  spray  vaporization  processes  play  an  important 
role  in  releasing  thermal  energy  and  driving  the  combustion  instability. 
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MOTIVATION 


• To  predict  the  nonlinear  instability  phenomena  in  liquid-fueled  rocket  en- 
gines. 

• To  gain  deeper  understanding  of  the  effects  of  the  design  parameters. 

• To  get  the  detailed  information  about  driving  mechanism  of  combustion  in- 
stabilities influenced  by  the  physical  processes  such  as  atomization,  vaporiza- 
tion , and  drop  breakup  Sc  collision. 

• To  develop  an  efficient,  accurate,  and  stable  numerical  model(pressure- 
based)  for  fast  transient  spray-combusting  flows. 
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APPROACH 


• Eulerian-Lagrangian  Formulation 

- Pressure-based  method 

- Applicable  to  all  speed  flows 

- Non-iterative  for  transient  calculations 

• Stochastic  Particle  Tracking  Technique 

- Delta  function  stochastic  separated  flow(SSF)  model 

- Stochastic  dispersion  width  transport (SDWT)  model 

• Equilibrium,  Non-equilibrium,  PDF  Combustion  Models. 

• Infinitive  &,  Effective  Conductivity  Vaporization  Model 

• Second-Order  Upwind  Scheme 
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ISSUES 


• Physical  processes  involved  in  the  driving  mechanism  of  combustion  instabil- 
ity. 

• Correlation  between  the  vaporization  response  characteristics  and  the  oscil- 
lating flowfield. 

• Prediction  capabilities  for  the  limiting  cycle  and  the  triggered  instability. 

• Effects  of  operating  conditions,  combustor  geometry,  and  stabilization  de- 
vices. 

• Validation  of  numerical  model  for  nonlinear  chamber  wave  phenomea. 
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SHOCK  TUBE  PROBLEM 


0.0  0.5  1.0 


FLAME  PROPAGATION  IN 


CLOSED  TUBE 


o o o o o 


OXYGEN  MASS  FRACTION  FUEL  MASS  FRACTION 

0.00000  0.11584  0.23167  0.000000  0.033162  0.066325 


I . QOms 

2.. 00ms 

3. . 0.0  ms 
4 , 00ms 
5 . 00ms 


X ( m ) 


Premixed  flame  propagation^  = 800 K,  D = 9.0x10  4m2 
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OXYGEN  MASS  FRACTION  FUEL  NASS  FRACTION 

0.00000  0.11542  0.23085  0.012259  0.041431  0.070603 


1 . 00ms 
2.0i 
3,01.... 

4 . 00ms 
5 


0-00  0.01  0.02  0.03  0.04  0. 

X Cm) 


Spray  flame  propagation^,  = 800 K,D  = 9.0x10 ~4m2/s,rk 


Spray  flame  propagation^  = 600 K,  D = 9.0x10 


0-01  0.02  0.03  0.04 

X ( m) 
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0.013489 


MRSS  FRACTION 
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l 
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Table  1.  Dimensions  of  Three  Liquid-Fuel  Rocket  Engines. 


Engine 

Li(m) 

L2(m) 

R(m) 

R|(m) 

©(deg) 

El 

0.3534 

0.2094 

0.2266 

0.1309 

65.44 

E2 

0.1767 

0.2094 

0.2266 

0.1309 

65.44 

E3 

0.1767 

0.1047 

0.2266 

0.1309 

47.57 
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(Atomization  + Vaporization  + Turbulent  Mixing  + Chemical 

Reaction  + etc) 

I 

Tc  (Overall  Combustion  Time  Scale) 

IF  tc  ~ Ta  (characteristic  acoustic  time  scale)  ► Combustion  Instability 


TIMEfms)  TIME (ms)  TIME  I ms) 


Figure  2 Pressure  oscillations  for  three  perturbation  levels; 
(rjktO=100/im,  </>=1.3,  X / — Y/i2=0.5,  E\) 


PRESSURE (atm) 

191.050  194.753  198. 457  202.160  205.863 


PRESSURE! atm) 


r*i0  = 100  /im 


16.0  27.0  36.0  15.0 

TIMEIms) 


0.0  6.0  16.0  27.0 

TIMEUbI 


16.0  2J 

TIMEIms ) 


Figure  4 Pressure  oscillations  for  three  initial  spray  conditions; 
(<£=1.3,  Cpert=  0.02,  X/L}=  0.0,  Y/R= 0.5,  E ,) 


♦ V ’-,11  • -f*  . . j 


100/im 


af.', 
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Figure  5 Spray  parcel  distribution  for  three  initial  spray  conditions 
($=1.3,  Cpert=  0.02,  Ei) 


| Mrii 


Figure  6 Contours  of  temperature  for  three  initial  spray  conditions; 
(0 — 1.3,  Cpert—  0.02,  X f L\—0.Q.  Ei) 


TlMEUs) 


Figure  7 Flow  oscillations  of  pressure,  velocity,  and  temperature; 

(r*,o=120pm,  0=1.3,  Cpert=  0.02,  X/Lx  = 0.0,  Y/R= 0.5,  Fi) 


1.0  6.0 
TIME (ms) 


< 0 6.0  8.0  10.0 
T I ME ( ms  ) 


4.0  6.0 

TIME (ms) 


Figure  8 Limit-cycle  flow  oscillations  of  pressure,  velocity,  and  temperat 
(r*,o= 120pm,  0=1.3,  X/Li=  0.0,  Y/R= 0.5,  Ei) 


Figure  9 Spray  parcel  distribution  and  temperature  contours; 
(rfc,o=120/xm,  <£=1.3,  Cpc rt=0.02,  Ei) 
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SUMMARIES 


Successful  predictions  for  the  unsteady  non-reacting,  flame-propagating,  and 
spray-combusting  flows. 

Variations  in  the  droplet  size,  the  combustor  length,  and  the  nozzle  length  of 
converging  section  have  a significant  effect  on  the  combustion  instability. 

Extension  to  transverse  mode  instability  analysis  and  incorporation  of  phys- 
ical submodels  dominantly  involved  in  the  driving  mechanism  of  combustion 
instability. 

Validation  of  numerical  model  for  linear  and  nonlinear  chamber  wave  phe- 
nomena. 
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Inverse  Design  of  a Proper  Number,  Shapes,  Sizes,  and  Locations  of 

Coolant  Flow  Passages 


George  S.  Dulikravich 

Department  of  Aerospace  Engineering,  233  Hammond  Building 
The  Pennsylvania  State  University,  University  Park,  PA  16802,  USA 


During  the  past  several  years  we  have  developed  an  inverse  method  that  allows  a 
thermal  cooling  system  designer  to  determine  proper  sizes,  shapes,  and  locations 
of  coolant  passages  (holes)  in,  say,  an  internally  cooled  turbine  blade,  a scram  jet 
strut,  a rocket  chamber  wall,  etc.  Using  this  method  the  designer  can  enforce  a 
desired  heat  flux  distribution  on  the  hot  outer  surface  of  the  object,  while 
simultaneously  enforcing  desired  temperature  distributions  on  the  same  hot  outer 
surface  as  well  as  on  the  cooled  interior  surfaces  of  each  of  the  coolant  passages. 
This  Constitutes  an  over-specified  problem  which  is  solved  by  allowing  the 
number,  sizes,  locations  and  shapes  of  the  holes  to  adjust  iteratively  until  the  final 
internally  cooled  configuration  satisfies  the  over-specified  surface  thermal 
conditions  and  the  governing  equation  for  the  steady  temperature  field. 


The  problem  is  solved  by  minimizing  an  error  function  expressing  the  difference 
between  the  specified  and  the  computed  hot  surface  heat  fluxes.  The  computed 

outer  surface  heat  flux  qout*P  will  not  be  the  same  as  the  specified  outer  surface 


heat  flux,  qSoUetc.  A properly  scaled  L-2  norm  of  the  difference  between  the 


specified  outer  surface  heat  flux,  q^ftc,  and  the  computed  outer  surface  heat  flux, 

qrp,  is  then  minimized  by  iteratively  changing  the  sizes,  shapes,  and  locations  of 

coolant  passages.  Starting  with  a large  number  of  guessed  holes,  all  unnecessary 
coolant  passages  are  efficiently  eliminated  when  their  sizes  reduce  below  a 
prespecified  minimal  allowable  value.  The  minimization  has  been  performed 
automatically  using  a standard  optimization  algorithm  of  Davidon-Fletcher- 
Powell.  Local  minimas  in  the  optimization  process  were  successfully  avoided  by 
changing  the  formulation  for  the  objective  function  whenever  the  local  minimas 
were  detected.  The  temperature  field  analysis  was  performed  using  our  highly 
accurate  boundary  integral  element  code  with  linearly  varying  temperature  along 
straight  surface  panels.  Examples  of  the  inverse  design  applied  to  internally 
cooled  turbine  blades  and  scram  jet  struts  (coated  and  non-coated)  having  circular 
and  non-circular  coolant  flow  passages  will  be  shown. 
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1.  Mathematical  model 


Steady  heat  conduction  in  internally  cooled  objects  is  modeled  as  a boundary 
value  problem  for  Laplace's  equation  over  a multiply-connected  domain. 
Assumptions  are:  - temperature  field  is  steady 

- solid  material  of  the  blade  is  thermally  isotropic. 

- thermal  expansion  is  neglected 
Governing  equation  is  Laplace’s  equation: 

V2T  = o|  (1) 


2.  Objectives 

Determine:  - the  exact  number  of  the  holes, 

- radii  of  the  holes, 

- locations  of  the  holes, 

such  that  relative  error  between  specified  and  computed  heat  fluxes  at  the  outer 
boundary  is  minimized. 

3.  Boundary  Conditions  - 111  Posed  Boundary  Value  Problem 

Both,  Dirichlet  and  Neumann  boundary  conditions  are  specified  on  the 
outer  boundary.  Such  an  overspecified  problem  can  be  solved  by  inverse  (design) 
approach.  The  problem  is  solvable  since  the  domain  is  multi-connected:  positions, 
shapes  and  dimensions  of  the  holes  will  provide  additional  degrees  of  freedom. 

4.  Constraints 

Besides  minimizing  the  heat  flux  error,  optimized  shape  has  to  satisfy  these 
constraints:  - minimum  distance  between  holes, 

- minimum  distance  between  holes  and  the  outer  boundary 
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5.  Objective  Functions 

Two  different  definitions  of  objective  function  were  used.  The  difference 

between  the  specified  and  the  heat  flux  and  heat  flux  obtained  by  the  current 

design  can  be  computed  as  a global  error: 

N 


(2) 


or  as  a local  error  in  heat  flux  at  each  node: 

V (qjc-qjr)2 

F2(x)  = Zj 


j=l 


(tf)- 


(3> 


Two  constraints  were  incorporated  into  the  objective  function  using  a barrier 
function 


M 


B(g(x),wb)  = ^ X 
° i=l 


N2 

X 


ds 


M 

X 


j=l  (^Dj-ds-rij  k=l  ^-dh-ri-rkjj 


(4) 


The  composite  objective  function  can  have  two  forms: 

Fi(g(x)»wb)  * Fj(x)  + B(g(x),wb),  1 = 1*2  (5) 

depending  whether  global  or  local  objective  function  is  used  for  its  evaluation. 
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2.  The  Optimization  Procedure 

The  optimization  procedure  consists  of  the  following  steps: 

(1)  Specify  shape  of  the  outer  surface  and  coating  of  the  turbine  blade. 

(2)  Specify  desired  temperature  Tjr  values  on  the  outer  and  inner  surfaces. 

(3)  Specify  desired  heat  flux  qjr  values  on  the  outer  surface. 

(4)  Specify  manufacturing  constraints: 

(i)  minimum  distance  ds  between  holes  and  the  outer  surface, 

(ii)  minimum  distance  d*1  between  any  two  neighboring  holes. 

(5)  Specify  initial  guess  for  the  number  of  holes,  M,  their  dimensions,  q,  and 

locations  of  the  centers  of  the  holes,  xj  and  yi.  Thus,  there  will  be  3xM 

design  variables  if  we  limit  ourselves  to  circular  holes  only. 

(6)  Using  the  Boundary  Element  Method,  the  Laplace’s  equation  for  a given 
domain  and  temperature  boundary  conditions  is  solved  and  heat  fluxes  at 
the  outer  boundary  are  computed.  The  Laplace's  equation  is  solved  3xM 
times,  ones  for  each  perturbed  design  variable  to  compute  the  gradient. 

(7)  Determine  relative  error  between  specified  and  computed  heat  fluxes  and 
evaluate  the  objective  function.  At  the  same  time  the  barierr  function  has  to 
be  evaluated  to  determine  the  composite  objective  function  Fj. 

(8)  Davidon-Powel-Fletcher  technique  is  used  to  find  the  new  values  of  design 
variables  repeating  the  optimization  procedure  from  the  step  (6)  until  the 
corresponding  composite  objective  function  F is  less  than  a prespecified 
value  If  the  dimension  of  a hole  becomes  less  than  a prespecified  value,  the 
hole  is  eliminated  from  further  optimization.  If  the  optimization  procedure 
stalls  in  a local  minimum  the  objective  function  formulation  is  changed 
from  Eq.  2 to  Eq.  3 while  continuing  with  optimization  from  the  step  (6). 
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Fig.  2b  Linear  element 


Pig.  1 Geometry  and  manufacturing  constraints 
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Fig.  1.1  Initial  configuration  (three  holes)  and  final  configuration  (one  large  centrally 
located  hole  and  two  dots  marked  with  arrows)  corresponding  to  a solution  with 


0.1%  integrated  flux  error. 


Fig.  1.2  Integrated  heat  flux  error  (L2  norm)  convergence  history  during  the 
optimization. 
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Fig.  1. 


Fig.  1. 


Coated  disk  problem  with  initially  ten  holes.  Convergence  history  shows  Five 
holes  are  reduced  to  zero.  Hole  elimination  method  was  not  used. 


Convergence  history  for  a circular  domain  with  initially  ten  holes.  Hole 
elimination  method  was  not  used.  Minimization  process  terminates  in  a local 
minimum. 


4 


Fig.  1.7  Coated  disk  problem  with  initially  ten  holes.  Hole  elimination  method  was  used 
together  with  objective  function  switching. 


Fig.  1.8  Convergence  history  for  a circular  domain  with  initially  ten  holes  when  hole 
elimination  method  is  applied  together  with  objective  function  switching. 
Discontinuities  represent  changing  of  the  objective  function. 
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NODE  NUMBER 


Fig.  1.9  Initial  ( ),  intermediate  (-•-•-)  and  final  ( ) heat  flux 

distribution  through  the  outer  boundary  for  a cylinder  with  ten  holes  initially. 
Hole  elimination  method  was  not  used. 


0 5 10  15  20  25  30  35 

NODE  NUMBER 


Fig.  1.10  Initial  ( ),  after  5 cycles  ( ) , after  10  cycles  ( ) and  final 

( ) heat  flux  distribution  through  the  outer  boundary  for  a cylinder  with 

ten  holes  initially.  Hole  elimination  method  and  objective  function  switching 
was  used. 
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Fig.  1.11  Total  CPU  time  (IBM  3090)  vs.  number  of  iterations  for  a circular  cylinder  with 
ten  holes  initially.  Hole  elimination  method  was  not  applied.  Total  number  of 
analysis  code  calls  (BEM  code)  was  1428. 


Fig.  1.12  Total  CPU  time  (IBM  3090)  vs.  number  of  iterations  for  a circular  cylinder  with 
ten  holes  initially  when  hole  elimination  method  was  applied  together  with 
objective  function  switching. 
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Fig.  1.13  A five-hole  coated  turbine  blade  from  which  thermal  boundary  conditions  were 
used  represents  an  actual  solution  for  the  case  of  the  turbine  blade  with  ten  holes 
initially 


Fig.  1.14  Initial  guess  for  a coated  turbine  blade  configuration  with  ten  holes  using 
thermal  boundary  conditions  from  the  five-hole  configuration 
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Fig.  1.15  Optimized  solution  for  initial  configuration  with  ten  holes.  Number  of  holes  is 
minimized  to  six,  where  the  sixth  hole  continues  to  shrink.  Hole  elimination 
method  was  used  together  with  objective  function  switching. 


SO  — r — 


0 5 >0  15  20  25  30 


NUMBER  OF  OPTIMIZATION  CYCLES  NO.  OF  BEM  CALLS  3775 


Fig.  1.16  Convergence  histoiy  for  a coated  turbine  blade  with  ten  holes  initially  when 
hole  elimination  method  was  applied  together  with  objective  function  switching. 
Discontinuities  represent  changing  of  the  objective  function. 
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2000 


Fig.  1.17  Initial  ( ),  after  5 cycles  ( ) , after  10  cycles  ( ) and  final 

( ) heat  flux  distribution  through  the  outer  boundary  for  a turbine  blade 

with  initially  ten  holes.  Hole  elimination  method  was  used  together  with 
objective  function  switching. 


Fig.  1.18  Total  CPU  time  (IBM  3090)  vs.  number  of  iterations  for  a turbine  blade  with 
ten  holes  initially.  Hole  elimination  method  was  applied  together  with  objective 
function  switching. 
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Fig.  1 Discretized  Ceramically  Coated  Scram  Jet  Combustor  strut 
With  prescribed  Temperatures  and  Outer  Surface  Heat  Flux 
chord  length  of  the  strut  : 19. 
maximum  thickness  of  the  strut  : 5. 


Fig.  2 Iteration  sequence  of  case  1 (norm  error  = 0.554  %) 


ERROR  PERCENTAGE 


Fig.  2.1  Initial  configuration  (an  off-center  inclined  almost  rectangular  hole) 
and  optimized  configuration  (one  large  centrally  located  hole)  for  one- 
hole  coated  disk  with  intermediate  hole  shapes. 


ITERATION 


Fig.  2.2  Convergence  history  of  the  coated  disk  with  one-hole  configuration. 
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Fig.  2.3  Initial  configuration  consisting  of  three  different  holes  (solid  line) 
and  their  inermediate  shapes  during  the  first  64  optimization  cycles  for  a 
coated  disk. 


\/v-  \v. 

\ ''  V\ 
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Fig.  2.4  Inermediate  shapes  of  the  three  different  holes  during  the 
optimization  cycles  65-121  for  a coated  disk. 


Fig.  2.5  Initial  configuration  (three  circular  holes)  and  optimized 
configuration  (ellipse,  rectangle,  and  a square)  for  a coated  turbine  blade 
with  intermediate  hole  shapes  (dotted). 


Fig.  2.6  Initial  configuration  (three  unequal  almost  rectangular  holes)  and 
optimized  configuration  (three  diffeently  sized,  positioned  almost 
rectangular  partially  constrained  holes)  for  the  coated  turbine  blade 
airfoil  with  intermediate  hole  shapes. 
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NUMERICAL  ANALYSIS  OF  THE  HOT-GAS-SIDE  AND  COOLANT-SIDE 
HEAT  TRANSFER  IN  LIQUID  ROCKET  ENGINE  COMBUSTORS 


T.S.  WANG 

Computational  Fluid  Dynamics  Branch 
NASA  - Marshall  Space  Flight  Center 
Marshall  Space  Flight  Center,  AL  35812 

V.  LUONG 
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Marshall  Space  Flight  Center,  AL  35812 


Abstract 


The  objectives  of  this  paper  are  to  develop  computational 
methods  to  predict  the  hot-gas-side  and  coolant-side  heat  transfer, 
and  to  use  these  methods  in  parametric  studies  to  recommend 
optimized  design  of  the  coolant  channels  for  regeneratively  cooled 
liquid  rocket  engine  combustors.  An  integrated  numerical  model 
which  incorporates  computational  fluid  dynamics  (CFD)  for  the  hot- 
gas  thermal  environment,  and  thermal  analysis  for  the  coolant 
channels,  was  developed.  The  model  was  validated  by  comparing 
predicted  heat  fluxes  with  those  of  hot-firing  test  and  industrial 
design  methods.  Parametric  studies  were  performed  to  find  a 
strategy  for  optimized  combustion  chamber  coolant  channel  design. 
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OBJECTIVES 


* TO  DEVELOP  COMPUTATIONAL  METHODS  FOR  THE 
PREDICTION  OF  THE  COUPLED  HOT-GAS-SIDE  AND 
COOLANT-SIDE  HEAT  TRANSFER  IN  A LIQUID  ROCKET 
ENGINE  COMBUSTOR 


* TO  PERFORM  PARAMETRIC  STUDIES  TO  RECOMMEND 
OPTIMIZED  DESIGN  OF  THE  COOLANT  CHANNELS  FOR 
REGENERATIVELY  COOLED  LIQUID  ROCKET  ENGINE 
COMBUSTORS 
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THE  AERO-THERMAL  MODEL 


★ CFD  MODEL  FOR  HOT-GAS-SIDE  ENVIRONMENT 

- AXISYMMETRIC  MCC  FLOWFIELD 

- FULLY  VISCOUS  FLOW 

- SHOCK  CAPTURING 

- SEVEN  SPECIES  EQUILIBRIUM  CHEMISTRY 

■k  SINDA  THERMAL  MODEL  FOR  LINER,  RIB,  AND  JACKET 

- THREE-DIMENSIONAL 

- VARIABLE  WALL  THICKNESS,  CHANNEL  DIMENSIONS 
AND  NUMBER  OF  CHANNELS 

- RADIATION  CORRECTED 

- WALL  TEMPERATURE  AND  THERMAL  GRADIENT 

★ SINDA  HYDRAULIC  MODEL  FOR  COOLANT  FLOW 

- COOLANT  TEMPERATURE  AND  PRESSURE  DROP 
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AERO-THERMAL  MODEL 


FLOW  CONDITIONS 
ASSUME  Tw 


r 

TOLERANCE  CHECk| 

TEST  CASES 


40K  CALORIMETER  THRUST  CHAMBER  TEST 
VALIDATION 


BASELINE  STANDARD  THROAT  SSME  MCC 
COMPARISON 


LARGE  THROAT  AMCC  DESIGN  PARAMETRIC 
STUDIES 


SCHEMATIC  OF  40K  TEST  CONFIGURATION 
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SUMMARY 


* AN  INTEGRATED  CFD/THERMAL  MODEL  HAS  BEEN  DEVELOPED  TO 
PREDICT  THE  HOT-GAS-SIDE  AND  COOLANT  SIDE  HEAT  TRANSFER 
FOR  LIQUID  ROCKET  COMBUSTION  CHAMBER 

★ MODEL  VALIDATED  FOR  40K  CALORIMETER  THRUST  CHAMBER 
TEST 

if  MODEL  COMPARED  FOR  BASELINE  STANDARD  THROAT  SSME  MCC 
HEAT  TRANSFER 

if  PERFORMED  LARGE  THROAT  AMCC  DESIGN  PARAMETRIC  STUDIES 
- INCREASED  ASPECT  RATIO  AND  NUMBER  OF  CHANNELS  REDUCE 
THE  WALL  TEMPERATURE  AND  THERMAL  GRADIENT 
• REDUCED  WALL  THICKNESS  REDUCES  THE  SURFACE  WALL 
TEMPERATURE 
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ABSTRACT 


| 


The  development  of  an  efficient  and  robust  grid  optimization  algorithm  is  presented.  This 
algorithm  is  developed  by  combining  the  best  characteristics  of  algebraic,  elliptic  and 
hyperbolic  grid  generation  techniques  (Ref.  1-3).  This  development  is  based  on  the 
following  observations  and  evaluations: 

Algebraic  systems  are  fast  and  economical. 

Precise  spacing  control  (well  distributed  grid)  is  always  achieved  with  algebraic  systems. 

Grid  generation  by  elliptic  system  is  always  smooth. 

Algebraic  system  may  cause  grids  to  everlap,  however,  elliptic  system  resist  grid  line  overlapping. 

Weighted  transfinite  interpolation  method  blended  with  Bezier,  B-spline  curves/surfaces  can  produce 
well-distributed,  orthogonai(at  Boundaries)  and  smooth  grids  (not  in  all  cases,  but  most  all). 

The  control  functions  can  be  formulated  to  achieve  boundary  orthogonality  and  spacing  control  (near 
solid  boundary  surface)  by  elliptic  generation  system. 

The  control  functions  can  be  formulated  to  accomplish  field  orthogonality  in  a given  computational 
direction  (h,  x,  or  z ) and  spacingcontrol  by  elliptic  generation  system  by  iteratively  updating  various 
terms  in  tne  generation  system.  This  is  very  time  consuming  especially  m three  dimensional  problems. 

Algebraic  systems  require  a high  degree  of  understanding  and  visual  user  interaction.  However,  elliptic 
systems  can  be  readily  adaptable  for  generalization.  This  is  extremely  useful  in  grid  adaptation. 

The  hyperbolic  system  preserves  the  orthogonality  at  the  solid  boundary  and  the  point  distribution  in 
the  field.  However,  its  applicability  is  restricted  to  external  flows  where  the  accurate  geometrical  shape 
of  the  outer  boundaries/surfaces  are  not  important  as  long  as  their  location  is  a certain  distance  away 
from  the  body.  Also  in  three  dimensional  applications  of  hyperbolic  system  the  grid  auality  is  directly 
influenced  by  the  characteristics  of  the  surfaces  associated  with  the  computational  domain. 


Computational  examples  representing  practical  internal  flow  configurations  are 
presented  to  demonstrate  the  success  of  this  algorithm. 
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Grid  Methods 


Direct  (Algebraic) 

• Fast  and  Economical 

• Precise  Spacing  Control 

• Propogation  of  Slope  Discontinuities 

• Interactive  User  Interface 

• Possible  Overlapping 

- can  be  avoided  1 

• High  Degree  of  Understanding  and 
Visual  User  Interaction 

• Orthogonality  and  Smoothness 

• Transfinite  : Lagrange,  Hermite,  Bezier, 

B-Splines,  NURBS 


Indirect  (PDES) 

• Time  Consuming 

• Distribution  Loss  I 

• Inherent  Smoothness 

• Iterative  Background  Crunching 

• Resistant  to  Grid  Line  Overlapping 

• Readily  Adaptable  for  Generalization 

• Competitive  Enhancement  of 
Smoothness,  Orthogonality, 
and  Concentration 

• Elliptic  Hyperbolic 
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APPROACH 


Objective:  Accomplish  othogonality  - smoothness 
without  any  distribution  loss. 


Work  hard  with  Algebraic 

- Precise  Spacing  Control  (Grid  Spacings,  Areas,  Volume) 

- Inexpensive  and  Fast 

- Interior  Bezier  Curve/Surface  Specification  for  Sub-blocks 

- Weighted  Transfinite  Lagrange  and  Hermite  Interpolation 

- Precise  Spacing  Control  (Grid  Spacings,  Areas,  Volume) 


Use  elliptic  for  a quick  fix 

- Smart  Forcing  Functions 

- 3-5  Iterations  (maximum) 
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Transfinite  Interpolation 
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HERMITE  TRANSFINITE 
INTERPOLATION 


Slope  Evaluation: 


Going  in  % direction 


*Vrn  - 
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Going  in  r|  direction 
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A Two  Dimensional  Elliptic 
Grid  System 

Q22(r^-^)+  gn(rim-lFrT1)-2  gi2r$n=0 

r=(x,y)  physical  space 

computational  space 

gn-A  • r^=xf+y| 
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922=.^  • It,=X^+y^ 

(j),'F  control  functions 
J 


Control  Functions 
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Cell  Area  Approach  (I) 

□ Evaluation  of  rn 
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Cell  Area  Approach  (II) 

□ Evaluation  of  r , 
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Cell  Area  Approach  (III) 


□ Evaluation  of  r 


-S.TP1 


(r^  • r^)ri  -(xgXq+ygr^ii-O 

H x rJL  =(x^-x^)ri=vt, 
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Grid  Spacing  Approach 

□ Evaluation  of  r,, 

(r^  • it i)ti  = + - 0 

VY*=922 

□ Evaluation  of  rt 

(r§  • i^) % =■  (j^Xq  + yy^H  - 0 

(rn  * rti)^  = + y!)^  = (922)^ 


□ Evaluation  of 


(ri;  • rn)ri  = (xy*,  + yyrn)ii  = 0 
(fn  • ft))1!  = (x*  + y^)ti  = (g22>Ti 
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Similarly  in  3D 

Working  on  £ = Constant  Surface 


r$  • r^  =0 

• rg  = 0 OR 

rj;  .•  (r-r,  X r8)  = V 


r^  • rn  = 0 
• rg  =0 
*V  = gii 
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Control  Functions  Using 

Metric  Term  cont'd 
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Similarly  In  3D 


. Id...  9n  v. 

4>  = 5-  -r-  ( ta  <- — — )) 

2 04  Q22  933 


* 2 dq  ' 'gugss" 


0=  \ -f  < In  ( 933  )) 
2 <fc  gii  g22 


J 


TtTT 


BEOMETW 


SSiii 


5SOO 


GRID 


GRIO 


GEOMETRY 


523 


525 


mmm 

mumm 


pmM 

■mbs 


526 


1 1 5 5! S5S  SjjJftlsS 
■ 


m 


§ 

■■I 


mm 


529 


530 


hi 


■■■I 


m 

iiiiib 

lailSfi 


minis 

IIIIUII 

■iiiiaMMi 


II  Si 

R 


■V.v 


m 

m 


gaS^P 


H88IIMB 


BW! 


B— —— 


532 


533 


534 


535 


536 


537 


540 


r 


FUTURE 
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FUTURE 


• Surface  Grid  Optimization  Using  NURBS 
Evaluation  and  Elliptic  System  Applications 
to  the  Parametric  Space 

• Full  3D  Applications  in  a Multiblock 
Environment 
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Enhancements  to  the  GRIDGEN  Structured  Grid  Generation 
System  for  Internal  and  External  Flow  Applications 


John  P.  Steinbrenner  and  John  R.  Chawner 
MDA  Engineering,  Inc. 

Arlington,  Texas 


GRIDGEN  is  a government  domain  software  package  for  interactive  generation  of  multiple 
block  grids  around  general  configurations.  Though  it  has  freely  available  since  1989,  it  has 
not  been  widely  embraced  by  the  internal  flow  community  due  to  a misconception  that  it 
was  designed  for  external  flow  use  only.  In  reality  GRIDGEN  has  always  worked  for  internal 
flow  applications,  and  GRIDGEN  ongoing  enhancements  are  increasing  the  quality  of  and 
efficiency  with  which  grids  for  external  and  internal  flow  problems  may  be  constructed. 

The  software  consists  of  four  codes  used  to  perform  the  four  steps  of  the  grid  generation 
process.  GRIDBLOCK  is  first  used  to  decompose  the  flow  domain  into  a collection  of 
component  blocks  and  then  to  establish  interblock  connections  and  flow  solver  boundary 
conditions.  GRIDGEN2D  is  then  used  to  generate  surface  grids  on  the  outer  shell  of  each 
component  block.  GRIDGEN3D  generates  grid  points  on  the  interior  of  each  block,  and 
finally  GRIDVUE3D  is  used  to  inspect  the  resulting  multiple  block  grid.  Three  of  these  codes 
(GRIDBLOCK,  GRIDGEN2D,  and  GRIDVUE3D)  are  highly  interactive  and  graphical  in 
nature,  and  currently  run  on  Silicon  Graphics,  Inc.  and  IBM  RS/6000  workstations.  The 
lone  batch  code  (GRIDGEN3D)  may  be  run  on  any  of  several  Unix  based  platforms. 

The  ease  of  flow  domain  decomposition  using  GRIDBLOCK  has  been  improved  through 
incorporation  of  edge  point  distribution  commands  and  a new  intermediate  construction 
entity  know  as  a domain.  Grid  point  dimensions  and  distributions  are  now  assigned  to  block 
boundary  curves  (connectors)  before  block  construction.  From  here,  block  subsurfaces  are 
defined  by  domains,  which  are  simply  a loop  of  connectors  that  represent  the  perimeter  of 
the  surface.  The  bounding  connectors  of  the  domain  and  the  grid  point  distributions  on  the 
connectors  provide  sufficient  data  for  the  automatic  initialization  of  surface  grid  points,  which 
may  be  later  refined  as  necessary  in  the  GRIDGEN2D  code.  Blocks  are  then  constructed 
by  grouping  domains  into  faces,  and  then  by  grouping  6 faces  into  a block.  Grouping  takes 
place  in  a point-and-click  environment,  and  the  reorientations  of  domains  and  faces  needed 
to  fit  these  components  into  the  developing  block  is  maintained  automatically  within  the 
code,  so  that  block  construction  may  proceed  in  an  intuitive  manner.  Further,  block  to  block 
interfaces  are  determined  automatically  on  the  domain  level,  and  domains  without  interblock 
connections  may  be  assigned  flow  solver  boundary  conditions  in  a graphical  interface. 

Surface  grid  generation  in  GRIDGEN2D  is  being  improved  with  the  addition  of  higher 
order  surface  definitions  (NURBS  and  parametric  surfaces  input  in  IGES  format  and  bicubic 
surfaces  input  in  PATRAN  Neutral  File  format)  and  double  precision  mathematics.  In 
addition,  two  types  of  automation  have  been  added  to  GRIDGEN2D  that  reduce  the  learning 
curve  slope  for  new  users  and  eliminate  work  for  experienced  users. 

Volume  grid  generation  using  GRIDGEN3D  has  been  improved  via  the  addition  of  an 
advanced  hybrid  control  function  formulation  that  provides  both  orthogonality  and  clustering 
control  at  the  block  faces  and  clustering  control  on  the  block  interior. 
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Overview  of  the  GRIDGEN  System 


• GRIDGEN  is  a series  of  four  codes 
for  the  generation  of  3D,  multiple 
block,  structured  grids. 

• GRIDBLOCK  (interactive): 
domain  decomposition. 

• GRIDGEN2D  (interactive): 

3D  surface  grid  generation. 

• GRIDGEN3D  (batch): 
volume  grid  generation. 

• GRIDVUE3D  (interactive): 
volume  grid  visualization. 
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Overview  of  the  GRIDGEN  System 

• The  interactive  codes  have  been  written  using  IRIS  GL  and  cur- 
rently run  only  on  Silicon  Graphics  4D  and  IBM  RS/6000  work- 
stations. 

• The  interactive  codes  also  require  24-bit  planes  and  Z-buffer. 

• GRIDGEN  documentation  consists  of  an  official  Air  Force  manual 
and  several  AIAA  and  AGARD  papers. 

— “The  GRIDGEN  3D  Multiple  Block  Grid  Generation  System”,  Vols.  I and  II, 
WRDC-TR-90-3022,  Flight  Dynamics  Lab.,  Wright- Patterson  AFB,  July  1990. 

— “Enhancements  to  the  GRIDGEN  System  for  Increased  User  Efficiency  and  Grid 
Quality”,  AIAA  paper  no.  92-0662,  AIAA  30th  Aerospace  Sciences  Meeting, 
January  1992. 

— “A  Structured  Approach  to  Interactive  Multiple  Block  Grid  Generation”,  from 
AGARD-CP-644  “Applications  of  Mesh  Generation  to  Complex  3-D  Configura- 
tions”, 1989. 
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Overview  of  the  GRIDGEN  System 

• Version  6,  USAF 

— Developed  for  USAF  at  Wright-Patterson  AFB,  1987-1990. 

— Technical  Supervision:  Dr.  Donald  W.  Kinsey. 

— Software  Distribution:  Lt.  John  Seo  (513)  255-2481,  to  U.S. 
government  agencies  and  U.S.  industry. 

• Version  8,  NASA  Langley 

- Currently  being  developed  for  NASA  Langley  Research  Cen- 
ter, 1991-1992. 

— Technical  Supervision:  Dr.  Robert  E.  Smith. 

— Software  Distribution:  Dr.  Jamshid  Abolhassani  (804)  864- 
5776  (Sept.  1992). 

• Version  9,  ? (currently  being  negotiated) 
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Databases  (Version  6) 


The  user  provides  GRID  GEN  with  a geometric 
description  (a  database)  of  the  configuration.  ^ 

A database  consists  of  a collection  of  \ \ _ 

patches  called  networks. 

Each  network  is  a 2D  array  of 
coordinate  data  on  the  configuration. 

The  networks  are  not  the  same  as 
the  surface  grid. 

The  database  may  be  obtained  from  a 

CAD  system,  an  external  user  program,  or  GRIDGEN2D. 


MDA  Engineering,  Inc. 
Arlington,  TX 


551 


GRIDBLOCK  (Version  6) 


• GRIDBLOCK  is  used  to  decompose  the  domain  surrounding  the 
database  into  blocks. 

• The  user  interactively  draws  3D  curves  (connectors)  that  define 
the  edges  of  each  block. 

- Straight  Line 

- Circular /Elliptical  Arc 

- Piecewise  Cubic 

- Line  on  Database 

• Connectors  may  be  drawn  in  any 
order  and  in  any  direction. 
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GRIDBLOCK  (Version  6) 


• The  user  groups  connectors 
into  blocks. 

• Blocks  may  contain  up  to 

12  edges;  singularities  are  allowed. 

• The  user  interactively  specifies 
computational  (£,  ?7,  C)  coordinate 
axes  and  number  of  points 

in  each  block. 
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GRIDBLOCK  (Version  6) 


• The  user  specifies  interblock 
connections  and  flow  solver  BCs. 

• Interblock  connections  must  be  set. 

- GRIDGEN2D  can  then  ensure 
consistency  between  blocks. 

- GRIDGEN3D  can  then  provide 
slope  continuity  across  interfaces. 

• The  user  may  set  TEAM  (USAF  Euler  solver)  flow  BC’s. 

- TEAM  restrictions  on  connections  are  checked  to  be  sure  grid 
is  compatible. 

- GRIDGEN3D  writes  the  connection  and  BC  data  in  TEAM 
format. 
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GRIDGEN2D  (Version  6) 

• GRIDGEN2D  is  used  to  generate  the  surface  grids  on  the  six 
faces  of  each  block  in  the  system. 

• It  may  also  be  used  to  generate  single  block  or  single  surface  grids 
without  running  GRIDBLOCK  first. 

• For  each  face  of  each  block... 

— Distribute  points  on  each  of  the  four  edges. 

— Initialize  surface  points  using  algebraic  methods. 

— Refine  surface  points  using  elliptic  PDE  methods. 
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GRIDGEN2D  (Version  6) 

• GRIDGEN2D  edge  point  distribution. 

- GRIDBLOCK  connectors  are  used  to  define  edge  shape  or  a 
new  edge  shape  may  be  drawn  interactively. 

- The  edge  may  be  divided  into  subedges  for  more  control  of 
point  distribution. 

— Grid  points  are  distributed  using... 

* Two-sided  tanh  (Vinokur)  stretching. 

* One-sided  sinh  and  tanh  stretching. 

* One-sided  geometric  progression. 

* Equal  spacing. 

* Copy  spacing  from  elsewhere  in  grid. 

* Cluster  to  edge  curvature. 
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GRIDGEN2D  (Version  6) 


• GRIDGEN2D  algebraic  methods. 

— Standard  TFI  with  computational  LaGrange  BF 

— Standard  TFI  with  arclength  based  LaGrange  BF 

— Ortho  TFI  with  computational  Hermite  BF 


- Polar  TFI 

- Re-distribution  methods. 

- Parametric  methods  to  fit 
the  grid  to  the  database. 
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GRIDGEN2D  (Version  6) 

• GRIDGEN2D’s  elliptic  PDE  methods. 

- Poisson’s  Equation  on  3D  surfaces  solved  using  pointwise  SOR 
with  Ehrlich’s  optimal  relaxation  factor. 

- Six  hybrid  control  function  formulations. 

— Five  solver  types. 

— Five  edge  BC  types 
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GRIDGEN2D  (Version  6) 

• GRIDGEN2D  elliptic  PDE  methods  cont’ 


• Hybrid  control 

functions  combine  j 

background  and 
foreground  control 
functions. 

• Background  control 
functions  tend  to 
influence  interior 
grid  points,  e.g. 

LaPlace,  Thomas  and 
Middlecoff,  or  Fixed  Grid 

• Foreground  control 
functions  tend  to 
influence  grid  points 
points  near  the  edges, 
e.g.  Sorenson. 

• (b)  Thomas  and 
Middlecoff 

• (c)  Sorenson 


(d)  Thomas  and 
Middlecoff  plus  Sorenson 


A y*  .. 

X As*' 

xAA'cf 


M 
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GRIDGEN2D  (Version  6) 

• GRIDGEN2D  elliptic  PDE  methods  cont’. 

— Three  conventional  solvers: 

* Solve  for  x,  y and  leave  z as  is. 

* Solve  for  x,  y , z. 

* Solve  for  x,  y and  interpolate  z 
from  the  database. 

— Two  parametric  solvers: 

* Solve  for  x , y,  z in  terms  of 
the  current  surface  shape. 

* Solve  for  x,  y,  z in  terms  of 
a database  network. 
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GRIDGEN2D  (Version  6) 

• GRIDGEN2D  tools 

— A face  may  be  divided  into  subfaces. 

* This  allows  the  shape  of 
and  distribution  of  points  on 
grid  lines  on  the  face  interior 
to  be  explicitly  set  by  the  user. 

- 8600  lines  of  help  text  may  be  accessed 
via  a browser. 

- There  is  a utility  to  graphically  move  any  point. 
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GRIDGEN3D  and  GRIDVUE3D  (Version  6) 

• GRIDGEN3D  is  a batch  code  written  for  a Cray  X/MP  running 
the  UNICOS  operating  system 

• It  may  easily  be  modified  for  other  UNIX  hardware. 

• Algebraic  methods  include  standard  TFI  with  computational  or 
arclength  based  LaGrange  blending  functions. 

• Elliptic  PDE  methods  solve  the  3D  Poisson  equations  using  point- 
wise  SOR  with  optimal  relaxation  factors. 

• LaPlace,  Thomas  and  Middlecoff,  Fixed  Grid,  or  Sorenson  control 
functions  are  available. 

• GRIDVUE3D  is  used  to  visualize  volume  grids  written  in  either 
GRIDGEN  or  PLOT3D  format. 
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Summary  of  GRID  GEN  enhancements  (Version  8) 

• Double  precision 

• Add  edge  point  grid  generation  to  GRIDBLOCK. 

- New  connector  shapes:  Cubic  on  surface  and  read  from  file. 

- New  distribution  function:  Monotonic  Quadratic  Rational  Spline 
(MQRS),  allows  a smooth  variation  of  grid  point  spacing  along 
the  connector  with  explicit  control  over  grid  point  locations  on 
the  interior. 

— Improved  editing  capability:  shape  or  number  of  points  can  be 
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Summary  of  GRIDGEN  enhancements  (Version  8) 

• Add  the  domain  entity  to  GRIDBLOCK. 

— Connectors  grouped  into  surfaces  called  domains.  Then  do- 
mains are  grouped  into  blocks. 

— A domain  may  be  a whole  face  or  only  a subface. 

- Point  to  point  interblock  connections  will  be  determined  au- 
tomatically. 

- Flow  BCs  will  be  set  by  graphically  picking  the  domain;  no 
more  typing  indices. 

- Algebraic  surface  grid  generation  will  be  performed  automat- 
ically; the  GRIDGEN2D  workload  is  drastically  reduced. 

- Changes  in  number  of  points  on  a single  connector  will  be 
propagated  semi-automatically  throughout  the  grid. 
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Summary  of  GRIDGEN  enhancements  (Version  8) 


• Add  the  domain  entity  to  GRIDBLOCK  cont’. 

— A domains  is  an  entity  between  the  connector  and  block  enti- 
ties in  the  GRIDBLOCK  hierarchy. 

— They  may  represent  a region  of  a single  flow  solver  BC  or  an 
interblock  connection. 

— The  user  creates  a domain  by  interactively  picking  the  indi- 
vidual connectors  in  a closed  loop. 

— Domains  must  be  computationally  rectangular. 


16  16  16  16 


This  domain  maps  into  This  domain  does  not  map 
a 45x31  region.  into  a rectangular  region. 
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Summary  of  GRIDGEN  enhancements  (Version  8) 

• Add  the  domain  entity  to  GRIDBLOCK  cont’. 

— Blocks  are  now  defined  by  six  faces,  rather  than  twelve  edges. 

- Faces  are  defined  by  at  least  one  domain. 

— Blocks  and  faces  are  checked  for  a consistent  number  of  points 
during  construction. 

- An  example  of  a face  consisting  of  four  non-unique  domains. 


Domain  B 
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Summary  of  GRIDGEN  enhancements  (Version  8) 


• Change  Number  of  Points  Utility 

— Low  level  changes  to  an  existing  grid  will  be  propagated  semi- 
automatically  throughout  the  entire  blocking  system. 

- Rather  than  edit  a journal  file,  the  code  will  do  most  of  the 
work  and  prompt  the  user  for  any  changes. 


- After  a change  in  the  number  of  points  on  the  indicated  connec- 
tor, the  code  would  ask  the  user  to  apportion  the  new  points 
across  the  affected  connectors. 
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Summary  of  GRIDGEN  enhancements  (Version  8) 

• GRIDGEN3D  upgrades. 

- Hybrid  control  functions  added,  for  improved  grid  quality. 

— Background  CF  control  on  interior  of  block,  e.g.  LaPlace, 
Thomas  & Middlecoff 

— Foreground  CF  control  near  faces,  e.g.  Sorenson 

- Hybrid  = Background  + Foreground 

* Compute  on  block  interior,  on  faces. 

* Calculate  = $/-$&  on  faces. 

* Interpolate  $a  from  faces  into  interior  using  exponentially 
decaying  blending  functions. 

* Sum  $ = $a  + $b 
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Summary  of  GRIDGEN  enhancements  (Version  8) 

• GRIDGEN3D  upgrades. 

- Grid  sequencing  added:  faster  convergence  rate  in  the  PDE 
solver. 

— Robustness  improved:  one  sided  differencing  based  on  the  sign 
of  the  control  function. 

— Efficiency  improved:  I/O  of  temporary  files  changed  to  reduce 
overhead. 

— Grid  quality:  several  quality  measures  are  written  to  a file  for 
visualization  using  PLOT3D  or  FAST. 
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Summary  of  GRIDGEN  enhancements  (Version  8) 

• GRIDGEN2D  Customization. 

— Goal:  reduce  the  effort  required  to  use  GRIDGEN2D. 

- Method:  eliminate  seldom-used  buttons  and  text  prompts. 

— Benefit:  fewer  keystrokes,  less  confusion. 

— Implementation:  verbosity  setting  and  preferencing. 

* Terse  verbosity  hides  obscure  prompts  from  the  user  (meant 
for  novices). 

* Preferencing  allows  the  user  to  pre-select  certain  options 
such  as  control  function  type  (meant  for  experts). 

• Double  precision  GRIDGEN2D. 
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Summary  of  GRIDGEN  enhancements  (Version  8) 

• Standardized  higher  order  surface  models  (databases). 

• PATRAN  Neutral  File. 

- Bicubics. 

• DT-IGES,  a simplified  form  of  the  IGES  standard. 

— Parametric  surface. 

- Rational  B-spline  surface. 

— Implementation  will  use  the  Navy  David  Taylor  Research  Cen- 
ter DT_NURBS  Library  of  surface  geometry  routines. 
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Continuing  GRIDGEN  improvement 

• Merge  GRIDBLOCK  and  GRIDGEN2D. 

- Less  user  confusion  (What  do  I do  in  which  code?) 

- Improved  usability  through  a single  GRIDBLOCK  style  GUI. 

— More  maintainability  through  elimination  of  duplicate  func- 
tionality. 

• Support  higher  order  surface  definitions  (databases)  in  standard 
file  formats  (e.g.,  NASA-IGES). 

• Increase  user  base  by  porting  to  other  hardware. 
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Conclusions 

• GRID  GEN  currently  provides  at  no  cost  a practical  and  well- 
tested  structured  grid  generation  capability. 

• Improvements  are  currently  being  added  to  the  government  do- 
main version  of  GRIDGEN. 

• GRIDGEN  will  remain  in  the  government  domain  well  into  the 
future. 

• MDA  Engineering  is  committed  to  supporting  GRIDGEN. 
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ABSTRACT 


Pregress  realized  in  the  development  of  a Computer  Aided  Grid  Interface  (CAGI) 
software  system  in  integrating  CAD/CAM  geometric  system  output  and/or  IGES 
files,  geometry  manipulations  associated  with  grid  generation  and  robust  grid 
generation  methodologies  is  presented.  CAGI  is  being  developed  in  a modular  fashion 
and  will  offer  fast,  efficient  and  economical  response  to  geomfetry/grid  preparation 
allowing  ability  to  upgrade  basic  geometry  in  a step-by-step  fashion  interactively  and 
under  permanent  visual  control  along  with  minimizing  the  differences  between  the 
actual  hardware  surface  descriptions  and  corresponding  numerical  analog. 


The  computer  code  GENIE  (Ref.  1-3)  is  used  as  basis.  The  Non-Uniform  Rational 
B-Splines  (NURBS)  representation  of  sculptured  surfaces  is  utilized  for  surface  grid 
redistribution.  The  computer  aided  analysis  system,  PATRAN,  is  adapted  as  a 
CAD/CAM  system.  The  progress  realized  in  NURBS  surface  grid  generation,  the 
development  of  IGES  transformer,  and  geometry  adaption  using  PATRAN  will  be 
presented  along  with  their  applicability  to  grid  generation  associated  with  rock 


propulsion  applications. 
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GEOMETRY  DEFINITIONS 

# Analytic 

# Drawings 

# Discretized  Points 

# Combo  (Combination  of  AboveThree) 

# CAD/CAM  Output 

# IGES 

$ CAD  Output  igi  MASAJ  Surface 

T Software  Data 

# Scale  Model 


SCULPTURED  SURFACES 

AP!^^ 

$ Fits  Given  Information 
$ Smoothness 

# Shape  Fidelity 

>!c  Parametorf^ftepresentation 

# Local  vs  Global  Schemes 

# Interactive  Design 

# Interactive  Viewing 


NGRB  Curves 


A NGRB  curve  c(u)  is  a piecewise  rational  curve  of  the  form 


<5(U)  - £ mdtN^u)  / £ WiNi,k(u ) , u e uk.u  Um*i 
1*0  («0 


defined  by 

• an  order  k (k  equalling  the  degree  of  the  polynomials  -1), 

• a set  of  3D  control  points,  { do,  , dm 

• a set  of  real  weights,  {<*> ©m}, 


* 
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NQRB  Curves  Gont'd 


• a set  of  real  knots  {“o.—.um*  * I u<s  ui»  1. 1 - ?•••(»*  + *■  • i)  V 

• B- spline  basis  functions  /■«,*(<*).  u « L,  u,.  1 1 . where 


«<.*(«)  - - nu-  ,(U)  ♦ u‘,>'u  n„t.k.  i(u) 

UUfc.l-Uf  V Ut+k-Ut+i 


Iiilm  f u < u*+ ! ; 

I 0,  otherwise. 


I - 0...  m, 


• and  curve  segments  c<u),  u€  ui(ui+i,  i - (k  - 1)  ...m 
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NGRB  Surfaces 


n m 


x(u,  v)  - 


X X <®) 

Jm0l»0  


n m 

1 1 


JmO  1-0 


G>W  ^ *(tt)ty  K0) 


1 

. 

, u« 

Uk  • 1 • ♦ 1 

0|-1.  »/*♦  1 

» 


defined  by 

• two  orders  fc  and  l (equalling  the  degree  of  the  polynomials  -1) 

• a set  of  3 D control  points  {<*0.0 <Vn  j- 
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NGRB  Surfaces  cont’d 

a set  of  real  u-knots  + * i ut  z ut+  j,  i - o...(m  + k - i)}, 

a set  of  real  u-knots  {oo»....on*  1 1 ojs:  dj*  ij  - o...(n  + 1 - 1)  J, 


B-spline  basis  functions  n^u),  u e uit  ut  «.  kj,  i - o...m 
^jc(u)  defined  as  for  the  curve  case, 


• B-spline  basis  functions  HjA°)>  0 « [°j*  i]./»  0"  ^ 

tyi(o)  defined  as  for  the  cure  case,  and 


• surface  segments  u uu  ijt  i - (*  - 

•0  elty,  + ll  J - (I  - !)...«• 


iit 
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£3  pto&d  g 


6E0nETRV 

Original  Geowet ry 
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QC 

UJ 

o 

/: 

UJ 

cc 


>o 

G <-0 

— s 
C-G 
c/o  — 

— a 
G1— 


> 

GOO 

—1—1 


a*-  log 

era.  — x 

13  0-1 


! 


a 

x 

UJ 


“TCSi"^  — LD 


oo 

3 uJ 
O M3 
— 1 LUX  — 

-crxji 

CXXJUJ 

C33I3 


X 


> 
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cnur^cNj 
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HELP 


GIPSI 


Graphically  Interactive  PATRAN  Structured  Grid  Interface 


INHUT  UIRECTIUE  UR  "ENU" 

SET.  riENU.  ON 

"MENU”  IS  SOL  ON  ILHS  ON  ). 
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IGES  Entitles:  (IGES  V.5.0) 


Circular  Arc  (type  100) 
Composite  Curve  (type  102) 
Conic  Arc  (type  104) 

• Parabola  (form  1) 

• Ellipse  (form  2) 

• Hyperbola  (form  3) 

• General  Equation  (form  0) 

Copious  Data 


Plane  (type  108) 

Line  (type  110) 

Parametric  Spline  Curve  (type  112) 

Parametric  Spline  Surface  (type  114) 

Point  (type  116) 

Ruled  Surface  (type  118) 


• Equal  Relative  Arc  Length  (form  0) 

• Equal  Relative  Parametric  Value  (form  1 ) 


• Center  line  (form  20-21)  Surface  of  Revolution  (type  120) 

• Section  (form  31-38)  Tabulated  Cylinder  (type  122) 

• Witness  line  (form  40)  Transformation  Matrix  (type  124) 


• Orthogonal  Matrix  (det  = 1 ) (form  0) 

right  handed  system 

• Orthogonal  Matrix  (det  ■ -1)  (form  1) 

left  handed  system 


Q[ 


IGES  Entitles:  (IGES  V.5.0)  cont. 


Rational  B-Spline  Curve  (type  1 26) 

• General  Parameters  (form  0) 

• Line  (form  1) 

• Circular  Arc  (form  2) 

• Elliptical  Arc  (form  3) 

• Parabolic  Arc  (form  4) 

• Hyperbolic  Arc  (form  5) 


Rational  B-Spline  Surface  (type  128) 

• General  (form  0) 

• Plane  (form  1) 

• Right  Circular  Cylinder  (form  2) 

• Cone  (form  3) 

• Sphere  (form  4) 

• Torus  (form  5) 

• Surface  of  Revolution  (form  6) 

• Tabulated  Cylinder  (form  7) 

• Ruled  Surface  (form  8) 

• General  Quadric  Surface  (form  9) 
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IGES  Entitles:  (IGES  V.5.0)  cont. 

Rational  B-Spline  CurveOffset  Curve  (type  130) 
Offset  Surface  (type  140) 

Boundary  Entity  (type  141) 

(set  of  curves  lying  on  surface) 


Curve  on  a Parametric  Surface 
Boundary  Surface 
Trimmed  Parametric  Surface 
User  Defined  Surface  Data  Form 


0 


1 


GEO'lETRV  2 

1 ■«  l 'i  m i 


flNfljLYS  IS 


HODEL  3 . fl^flLVZE  -I,  RESULTS  . 

"'M.l  ■ I ' t J ' • I 


"aaaaaaaaaaaaaaaaaaa aaaaaa aaaaaa aaaa aaaaaa a a 


**  # OF  POINTS  CONVERT  : 21 
**  # OF  LINES  CONVERT  : 0 
**  # OF  CIRARC  CONVERT  : 0 
**  # OF  PCCURVE  CONVERT  : 17 


******************************************* 


******************************************* 


**  INQUIRE 

A* 

AA 

AA 

AA 

AA 

AA AAA AAA AAA 


FUNCTION  FOR  THE  CONVERT  IGES  ** 

1 ^ DHTIITC!  ** 


1>. 

POINTS 

A X 

2> . 

LINES 

A A 

3> . 

CIRARC 

A A 

4>  . 

PARAMETER  CURVES 

A X 

5>  . 

QUIT  INQUIRE 

* A 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa* 


FIRST  , USE  PAT RAN  2 . 5 TO  CRESTE  THE  PC  CURVES 
(THE  LEFT  UPPER  WINDOW  ) . 

OUTPUT  THE  PC  CURVES  TO  IGES  FILES  . 

USE  THE  CAGI -IGES  TRANSLATOR  TO  CONVERT  THE  IGES  FILE 
BY  RUNNING  CAGI  AND  READ  THE  M PAT RAN . IGS . 1"  AS  INPUT  . 
LEFT  DOWN  WINDOW  SHOW  THERE  ARE  21  POINTS  AND  17  PC 
CURVES  BEEN  CONVERTED  . 

PLOT  THE  PC  CURVES  ON  THE  RIGHT  UPPER  WINDOW  . 
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ccmflNO  sfliE  nisnov  nmuF  limwl  uifl  sfi/shol 


INPUT  DIRECTIVE  OR  "E,'«D" 
E 

■ H 'H  t ' ■ •*  ■ ' 




PLEASE  KEYIN  THE  IGXS  FILE  NEME  .... 

body wing . igs 

AAAAAAAAA A A A A A A A A A A A A A A A A A AAAAAAAAAAAAAAAAA 


aa  | OF  POINTS  CONVERT  : 19 
**  # OF  LINES  CONVERT  : 6 
**  # OF  CIRARC  CONVERT  : 0 
aa  | OF  PCCURVE  CONVERT  : 6 
aa  | OF  PC  SUR  CONVERT  : 10 
**  f OF  BODY  OF  ROV.  0 
**  # or  CONIC  ARC  0 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

**  INQUIRE  FUNCTION  FOR  THE  CONVERT  ICES  ** 
aa  1>.  POINTS  ** 


1>.  FIRST  , USE  PATRAN  2 . 5 TO  CREATE  THE  PC  SURFACE 
(THE  LEFT  UPPER  WINDOW  ) . 

2> . TRANSLATE  THE  PARAMETER  SURFACES  TO  ICES  FILES  . 

3> . USE  THE  CAGI-IGZS  TRANSLATOR  TO  CONVERT  THE  ICES  FILE 
BY  RUNNING  CAGI  AND  READ  THE  "BODYWING. IGS"  AS  INPUT 
4> . LEFT  DOWN  WINDOW  SHOW  THERE  ARE  19  POINTS  AND  6 PC 
CURVES  10  SURFACE  AND  6 LINES  BEEN  CONVERTED  . 

5> . PLOT  THE  PC  SURFACE  ON  THE  RIGHT  UPPER  WINDOW  . 

D 
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HELP 


geometry 


21*20  ** 

20x20 

: 6RID  24 

20x20 

* 6RI0  25 

20X20 

GRID  26 

20x20 

GRID  27 

20x20 

GRID  28 

20x20 

GRID  29 

20X20 

GRIO  30 

20X20 

GRIO  31 

20X20 

GRID  32 

20x20 

GRIO  33 

20X20 

6RID  34 

20x20 

6RI0  35 

20X20 

GRIO  36 

20x20 

6R10  37 

609 


610 
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2DF0IL.  1GES 

, , 1 6HPATCAD  IGES  FILE, lOHpatran . igs, 
20HPDA/PATRAN,  Rel  2.4  , 3H4 .0, ,,,,,, 1 . 0, 
, , 13H910313 .235011, , , , , 6, 0; 


1« 4HINCH, 


SUUUUUU1 

G0000001 

G0000002 

G0000003 


116, 

116, 

116, 

116, 

116, 

116, 

116, 

116, 

110, 


116 

1 

2 

2 

0 

0 

0 

OOOOOOOOODOOOOOOl 

116 

0 

0 

1 

0 

0 

OPQINT 

1D0000002 

116 

2 

2 

2 

0 

0 

0 

OOOOOOOOODOOOOOO3 

116 

0 

0 

1 

0 

0 

OPOINT 

2D0000004 

116 

3 

2 

2 

0 

0 

0 

00000000000000005 

116 

0 

0 

1 

0 

0 

OPOINT 

3D0000006 

116 

4 

2 

2 

0 

0 

0 

OOOOOOOOODOOOOOO? 

116 

0 

0 

1 

0 

0 

OPOINT 

400000008 

116 

5 

2 

2 

0 

0 

0 

00000000000000009 

116 

0 

0 

1 

0 

0 

OPOINT 

5D0000010 

116 

6 

2 

2 

0 

0 

0 

OOOOOOOOOOOOOOOll 

116 

0 

0 

1 

0 

0 

OPOINT 

600000012 

116 

7 

2 

2 

0 

0 

0 

00000000000000013 

116 

0 

0 

1 

0 

0 

OPOINT 

700000014 

116 

8 

2 

2 

0 

0 

0 

00000000000000015 

116 

0 

0 

1 

0 

0 

OPOINT 

600000016 

110 

9 

2 

1 

0 

0 

0 

00000000000000017 

110 

0 

0 

2 

0 

0 

OLINE 

100000018 

124 

11 

2 

1 

0 

0 

0 

000010000D0000019 

124 

0 

0 

4 

0 

0 

OMATRIX 

200000020 

100 

15 

2 

1 

0 

0 

• 19 

00000000000000021 

100 

0 

0 

3 

0 

0 

OCIRC  UC  2D0000022 

110 

18 

2 

1 

0 

0 

0 

00000000000000023 

110 

0 

0 

2 

0 

0 

OLINE 

3D0000024 

110 

20 

2 

1 

0 

0 

0 

00000000000000025 

110 

0 

0 

2 

0 

0 

OLINE 

4D0000026 

110 

22 

2 

1 

0 

0 

0 

00000000000000027 

110 

0 

0 

2 

0 

0 

OLINE 

5D0000028 

110 

24 

2 

1 

0 

0 

0 

00000000000000029 

110 

0 

0 

2 

0 

0 

OLINE 

600000030 

124 

26 

2 

1 

0 

0 

0 

000010000D0000Q31 

124 

0 

0 

4 

0 

0 

OMATRIX 

8D0000032 

100 

30 

2 

1 

0 

0 

31 

00000000000000033 

100 

0 

0 

3 

0 

0 

OCIRC  ARC  800000034 

110 

33 

2 

1 

0 

0 

0 

00000000000000035 

110 

0 

0 

2 

0 

0 

OLINE 

9D0000036 

114 

35 

2 

1 

0 

0 

0 

00000000000000037 

114 

0 

0 

20 

0 

0 

OSPLSRF 

100000030 

114 

55 

2 

1 

0 

0 

0 

00000004000000039 

114 

0 

0 

20 

0 

0 

OSPLSRF 

200000040 

114 

75 

2 

1 

0 

0 

0 

00000000000000041 

114 

0 

0 

20 

0 

‘o 

OSPLSRF 

600000042 

O.OOOOOOOOOE+OO, 
-10.0000000 
0.000000000E+00, 
10.0000000 
20.0000000  , 
20.0000000 
10.0000000 
1.00000000 
-10.0000000 
O.OOOOOOOOOE+OO, 

124,  -1.00000000 
O.OOOOOOOOOE+OO, 
O.OOOOOOOOOE+OO, 
O.OOOOOOOOOE+OO, 


O.OOOOOOOOOE+OO, 
O.OOOOOOOOOE+OO, 
10.0000000  , 
10.0000000 
10.0000000 
O.OOOOOOOOOE+OO, 
O.OOOOOOOOOE+OO, 
1.00000000 


O.OOOOOOOOOE+OO; 

O.OOOOOOOOOE+OO; 

O.OOOOOOOOOE+OO; 

O.OOOOOOOOOE+OO; 

O.OOOOOOOOOE+OO; 

O.OOOOOOOOOE+OO; 

O.OOOOOOOOOE+OO; 

O.OOOOOOOOOE+OO; 

O.OOOOOOOOOE+OO, 


, O.OOOOOOOOOE+OO, 

Q.OOOOOOOOOE+OO,  O.OOOOOOOOOE+OO; 

V,  O.OOOOOOOOOE+OO,  O.OOOOOOOOOE+OO, 
O.OOOOOOOOOE+OO,  1.00000000  , 

O.OOOOOOOOOE+OO,  O.OOOOOOOOOE+OO, 
-1.00000000  , O.OOOOOOOOOE+OO; 

100,  O.OOOOOOOOOE+OO,  O.OOOOOOOOOE+OO,  O.OOOOOOOOOE+OO, 
10.0000000  , O.OOOOOOOOOE+OO,  O.OOOOOOOOOE+OO, 

10.0000000 

110,  O.OOOOOOOOOE+OO,  10.0000000  , O.OOOOOOOOOE+OO, 


1P0000001 

3P0000002 

5P0000003 

7P0000004 

9P0000005 

11P0000006 

13P0000007 

15P0000008 

17P0000009 

17P0000010 

19P0000011 

19P0000012 

19P0000013 

19P0000014 

21P0000015 

21P0000016 

21P0000017 

23P0000018 


- - v 

Vpoint',  point"  2,  rl— 10.00000  \ 0 . 0000000e+00,  0 . 0000000e+00  $ 

$ ' point' , point-  3,  rl-0 . 0000000e+00,  10.00000  , 0 . 0000000e+00  $ 

S'point',  point-  4,  rl-  10.00000  , 10.00000  , 0 . 0000000e+00  $ 

$' point',  point-  5,  rl-  20.00000  , 10.00000  , 0 . OOQOOOOe+OO  $ 

$' point' / point-  6,  rl-  20.00000  , 0 . OOOOOOOe+CK),  0 . OOOOQOOe+OO  $ 

$' point',,  point-  7,  rl-  10;  00000  , 0. 0000000e+00,  0 . ©&&&Q©0«+00  $ 

$'point',  point-  8,  rl-  1.000000  , 1.000000  , 0 . QOOOOOOe+OO  $ 

$'line',  points-  10,  rl— 10.00000  ,0.0000000e+00,0.0000000e+00, 

r2— 0 . 0000000e+00,  0 . 0000000e+00, 0 . 0000000e+00,  coreout-  1 $ 

$'oonicur',  type-' circle' ,points-20,  angle-  0.00,  90 . 00, radius-  10.00, 

coreout-  2$ 

$'trans', corein-  2,  rl-  0.00,  0.00,  0.00,r2-  0.00,  0.00,  0.00, 

cosines— 1.0,  0.0,  0.0,  0.0,  1.0,  0.0,  0.0,  0.0, -1.0, coreout-  2$ 
$'line' , points-  10,  rl-0.0000000e+00,  10.00000  ,0.0000000e+00, 

r2-  10.00000  , 10.00000  , 0 .0000000e+00,  coreout-  3 $ 

$' line' , points-  10,  rl-  10.00000  , 10.00000  ,0.0000000e+00, 

r2-  20.00000  , 10.00000  , 0 .0000000e+00,  coreout-  4 $ 

$'line' , points-  10,  rl-  20.00000  ,0.0000000e+00,0.0000000e+00, 

r2-  20.00000  , 10.00000  ,0.0000000e+00,  coreout-  "5  $ 

$'line',  points-  10,  rl-  10.0Q000  ,0.0000000e+00,0.0000000e+00, 

r2-  20.00000  , 0.0000000e+00,0 .0000000e+00,  coreout-  6 $ 

$'conicur',  type-' circle' , point s-20,  euigle-  0.00,  90. 00, radius-  1.00, 

coreout-  7$ 

$'trans', corein-  7,  rl-  -1.00,  0.00,  0.00,r2-  0.00,  0.00,  0.00, 

cosines— 1.0,  0.0,  0.0,  0.0,  1.0,  0.0,  0.0,  0.0, -1.0, coreout-  7$ 
$'line' , points-  10,  rl-  1.000000  , 1.000000  ,0.0000000e+00, 

r2-  10.00000  , 0 . 0000000e+00, 0 . 0000000e+00,  coreout-  8 $ 

$' combine' , corein-1,-8, form-'plot3d' , fileout-4$ 

$'end'$ 


GENIE+  + 


• Simple  Minded 

• Portable 

• Modular 

• Journal  File  Execution  Control 

• FORTRAN  SC 

• Extensive  Error  Checking 

• Easy  Customization 

• (4-)  Genie  Shortcomings 

- Spacings 

- Geometry/Grid  Manipulations 

- On  Line  Storage 

- B Spline  - NURBS  Applications 

- Geometry  Interface 

- Interactive  Visualization 
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USING  ADAPTIVE  GRID  IN  MODELING  ROCKET  NOZZLE  FLOW 

i 

By  Alan  S.  Chow  * and  Kang-Ren  Jin** 

♦NASA/Performance  Analysis  Branch,  Marshall  Space  Flight  Center,  AL  35812. 

**Department  of  Civil  Engineering,  Mississippi  State  University,  Mississippi  State,  MS  39762. 

ABSTRACT 

I The  mechanical  behavior  of  a rocket  motor  internal  flow  field  results  in  a system  of 

nonlinear  partial  differential  equations  which  cannot  be  solved  analytically.  However,  this  system 
of  equations  called  the  Navier-Stokes  equations  can  be  solved  numerically.  The  accuracy  and  the 
convergence  of  the  solution  of  the  system  of  equations  will  depend  largely  on  how  precisely  the 
sharp  gradients  in  the  domain  of  interest  can  be  resolved.  With  the  advances  in  computer 
technology,  more  sophisticated  algorithms  are  available  to  improve  the  accuracy  and  convergence 
of  the  solutions.  An  adaptive  grid  generation  is  one  of  the  schemes  which  can  be  incorporated 
into  the  algorithm  to  enhance  the  capability  of  numerical  modeling.  It  is  equivalent  to  putting 
intelligence  into  the  algorithm  to  optimize  the  use  of  computer  memory.  With  this  scheme,  the 
finite  difference  domain  of  the  flow  field  called  the  grid  does  neither  have  to  be  very  fine  nor 
strategically  placed  at  the  location  of  sharp  gradients.  The  grid  is  self  adapting  as  the  solution 
evolves.  This  scheme  significantly  improve  the  methodology  of  solving  flow  problems  in  rocket 
nozzle  by  taking  the  refinement  part  of  grid  generation  out  of  the  hands  of  computational  fluid 
dynamics  (CFD)  specialists  and  place  it  into  the  computer  algorithm  itself. 

PRECEDING  PAGE  BLANK  NOT  FILMED 
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OBJECTIVE 


• To  develope  a user-friendly  solution- 
adaptive  grid  generator  that  will  simplify 
grid  generation  process  so  that  a 'perfect’ 
grid  can  be  generated  everytime  without 
the  intervention  of  CFD  experts. 
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dU  M JF  ac 

dt  ax  a y a z 


Os 

00 


m 


u- 


p 

pu 

pv 


E « 


P***  +P~rm 
P*>-T*y 

pUW  - Tst 


pw 

E, 


(Et  + p)u  - UTjtx  - vrXy  - *"*«  + Qx 


pv 

pw 

PUV  - Txy 

pttw  — rm 

PU5  + p — Tyy 

c* 

pVW-Tyg 

PVW-Tyg 

pw*  +p-  Tm 

(Ef  + P)V  ~ UTXy  -VTyy—  WT  yg  + Q y^ 

(Et+  p)w-UTm  -vTyg  — wrm  +<7f_ 

MSFC 
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dS  dF,  1 dGj 
dt  + dX,  ~ Re  dXj 

where  S is  a vector  containing  the 
conservation  variables, 


The  Fj  vectors  represent  the  inviscld 
flux  vectors, 

r ***> 

Fi  ~ pu.Uj  + P,6„ 

■ (E  + Pt)uj  . 

and  Gj  vectors  are  the  viscous  flux 
vectors 


0 


u*rj*  - J 
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GRID  GENERATION  METHODS 

• Complex  variables  (Conformal  Mapping) 

• Algebraic 

• Partial  Differential  Equations  (PDE) 
-Elliptic 

-Hyperbolic 
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■ constant 

W(x)  : weight  function 

x(£):  point  distribution 

xttW  + xtW(  - 0 


xtt  + 

P - _fil 
x( 


1.2,3) 


W - 1 + |Vp| 
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The  Adaptive  Grid  at  Time  Step  - 500 


for  Inviscid  Flow. 
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The  Adaptive  Grid  at  Tima  Step 


1000  for  Inviacid  Flow. 
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The  Adaptive  Grid  at  Time  Step  « 2000  for  Inviecod  Flow 
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WEIGHT  FUNCTION 
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The  Weight  Function  Distribution  st  Time  Step  ~ 300  for 
Viscous  Flow. 


631 


MSFC 

GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  ALABAMA 


NASA 

NATIONAL  AERONAUTICS 

And 

SPACE  ADMINISTRATION 


MSFC 

GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE.  ALABAMA 


NASA 

NATIONAL  AERONAUTICS 
And 

SPACE  ADMINISTRATION 


633 


MSFC 

GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  ALABAMA 


9I>91 


NASA 

NATIONAL  AERONAUTICS 

And 

SPACE  ADMINISTRATION 


634 


MSFC 

GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE.  ALABAMA 


NASA 

NATIONAL  AERONAUTICS 

AND 

SPACE  ADMINISTRATION 


MSFC 

GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  ALABAMA 


O' 

u> 

Ul 


NASA 

NATIONAL  AERONAUTICS 
And 

SPACE  ADMINISTRATION 


GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  ALABAMA 


NATIONAL  AERONAUTICS 
AND 

SPACE  ADMINISTRATION 


637 


MSFC 

GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  ALABAMA 


COWTOU*  urn s 
0.04000 
a.OGOCD 
0.00000 
0. 100GD 
0.12000 
a.  14003 
0.16000 
0. I HOOD 
0.20000 
0.22000 
0.24003 
0.26000 
0.26000 
0.30000 
0.32000 
0.34000 
0.36000 
0.36000 
0.40000 
0.42000 
0.44000 
0.46000 
0.46000 
0.50000 
0.52000 

0. 54000 
0,56000 
osKm 
0.60000 
0.62000 
0.64000 
0.86000 
0.66000 
0.7OUOO 


PRESSURE 


X 


INVISCID  NOZZLE 


NASA 

NATIONAL  AERONAUTICS 
And 

SPACE  ADMINISTRATION 


FLOW 


MSFC 

GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 


NASA 

NATIONAL  AERONAUTICS 
AND 

SPACE  ADMINISTRATION 


2.0  3-0 

INVISCID  NOZZLE  FLOW 


MSFC 

GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  ALABAMA 


NASA 

NATIONAL  AERONAUTICS 
AND 

SPACE  ADMINISTRATION 


CONTOUR  LEVELS 
0.10000 
0.13000 
0. 14000 

8.16000 
.10000 
0.30003 
0.23000 
0.3«XZ3 
0.20000 
0.28000 
0.30000 
0.32000 
0.34000 
0.36000 
0.38000 
0.40000 
0.42000 
0.44000 
0.46000 
0.48000 
0.50000 
0.52000 
0.54000 
0.56000 
0.56000 
0 80000 
0.63000 
0.64000 
0.60000 
0.68000 
0.70000 


PRESSURE 


1.00.10’ 
101 .51 


VISCOUS  NOZZLE  FLOW 


M 3.0  4.0 

VISCOUS  NOZZLE  FLOW 


641 


MSFC 

GEORGE  C.  MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  ALABAMA 


NASA 

NATIONAL  AERONAUTICS 
AND 

SPACE  ADMINISTRATION 


WEIGHT  FUNCTION 


IUImM  ( Mlllt 


The  Weight  Function  Distribution  at  Tiae  Step  » 5500  for 
Viscous  Flow. 
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The  Weight  Function  Distribution  at  Time  Step  - 1000  for 
Viscous  Flow. 
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CONCLUSIONS  & RECOMMENDATIONS 

/2-D  Solution-Adaptive  Grid  Generator  has 
been  completed  and  demostrated. 

V It  simplifies  grid  generation  and  makes 
better  use  of  computer  and  human  resources. 

y Should  be  verified  on  various  CFD  codes  to 
ensure  robustness. 

/ Continue  development  in  3-D  and  Time-accurate 
Solution-Adaptive  Grid  Generator. 
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Abstract  for  the  Tenth  CFD  Working  Group  Meeting: 

Complex  Three-Dimensional  Internal  Flows 
in  the  ASRM  and  RSRM  Aft  End  Segments 

Presented  By:  Dr.  Edward  J.  Reske 

Dr.  Dana  F.  Billings 
Ms.  Joni  W.  Cornelison 

Results  from  computational  fluid  dynamic  analyses  for  complex 
three-dimensional  internal  flows  in  the  Advanced  Solid  Rocket  Motor  (ASRM) 
and  Redesigned  Solid  Rocket  Motor  (RSRM)  are  presented.  In  particular,  a 
parametric  study  for  the  case  of  a gimballed  nozzle  in  these  motors  at 
various  burn  times  and  gimbal  angles  is  presented.  The  resultant 
pressure  fields  are  used  to  determine  the  location  of  the  center  of 
pressure  and  hinge  moments  due  to  the  interal  flow  for  these  geometries. 
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ASRM  AFT  END  FLOW  ANALYSIS 


ASRM  AFT  END  ANALYSIS 


\ 


•Objective 

* Characterize  flow  environment  in  aft  end  of  ASRM 

• Purpose 

Hinge  moments  due  to  internal  flow  for  a gimballed  nozzle 
> Nozzle  performance 

- Heat  transfer  for  insulation  sizing 

•Approach 

-Axisymmetric  analyses 

CM'lNT  (48K  and  J24K  grid  points) 

--  FDNS  (14K  grid  points) 

- Three-dimensional  gimballed  nozzle  analysis 

- FDNS3D  (14K  X 26  planes  - 366K  grid  points) 

• Results 

- Axisymmetric  analysis  complete 

- 3-D  gimballed  nozzle  analyses  nearing  completion 

— - 2L J 
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Greed  Outline  = RSRM  at  19  seconds 
Rad  Outline  * RSRfl  at  19  seconds 


: t>?0xR4  GRin  1 

>*0x64  GRID  2 

: 220 k6 4 CRID  3 

220x64  GRID  4 
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COMPLEX  3D  FLOWS  IN  THE  ASRM 
AND  RSRM  AFT  END  SEGMENTS 


HINGE  MOMENTS  AND  LOADS 

MOTOR  BURN  GIMBAL  HINGE  NORMAL  AXIAL  CENTER  OF  PRESSURE 


TIME 

sec. 

ANGLE 

deg. 

MOMENT 

in.-lb. 

LOAD 

lb. 

LOAD 

lb. 

wrt  pivot  wrt  throat 
in. 

ASRM 

19 

4 

627  K 

39.6  K 

2.52  M 

-15.9 

+1.7 

ASRM 

19 

8 

1.28  M 

80.8  K 

2.51  M 

-16.0 

+1.6 

RSRM 

19 

4 

730  K 

49.5  K 

3.18  M 

-14.7 

+2.9 

ASRM 

115 

4 

150  K 

6.9  K 

1.25  M 

-21.7 

-4.1 

ASRM 

115 

8 

546  K 

12.8  K 

1.25  M 

-42.6 

-25.0 

RSRM 

114 

4 

165  K 

14.7  K 

0.96  M 

-11.2 

+6.4 

Note:  All  the  above  hinge  moments  are  non-restoring  torques.  The  axial  load  acts 
along  the  axis  of  symmetry  of  the  nozzle,  and  the  normal  load  acts  in  the 
direction  perpendicular  to  this  axis,  with  both  components  acting  in  a 
direction  away  from  the  motor.  The  center  of  pressure  is  determined  by 
finding  a point  on  the  nozzle  axis  of  symmetry  where  the  torque  vanishes. 

A negative  value  indicates  that  it  is  upstream  of  the  reference  point,  whereas 
a positive  value  indicates  that  it  is  downstream  of  the  reference  point. 
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A simple  model  for  calculating  the  hinge  moment  about  a 
pivot  point  that  is  shifted  relative  to  the  nominal  location: 

Tz  = TDZ+yFM-  xFa 


Pz  ' N A 

where,  using  a body-fixed  coordinate  system, 
Tz  = torque  about  the  new  pivot  point; 


T] pz=  tor9ue  about  the  nominal  P'vot  P°inb 
y = axial  displacement  of  the  new  pivot  point; 


x = normal  displacement  of  the  new  pivot  point; 


F,  = normal  load; 
F = axial  load. 
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An  Analysis  of  the  Flow  Field  in  the  Region  of  the  ASRM  Field  Joints 

Richard  A.  Dill,  ERC  Incorporated 
Harold  R.  Whitesides,  ERC  Incorporated 

Abstract 

The  flow  field  in  the  region  of  a solid  rocket  motor  field  joint  is  very  important  since  fluid 
1 dynamic  and  mechanical  propellant  stresses  can  couple  to  cause  a motor  failure  at  a joint.  This  paper 

presents  an  examination  of  the  flow  field  in  the  region  of  the  ASRM  field  joints.  The  analyses  were 
performed  as  a first  step  in  assessing  the  design  of  the  ASRM  forward  and  aft  field  joints  in  order  to 
assure  the  proper  operation  of  the  motor  prior  to  further  development  or  test  firing. 

I 

The  analyses  discussed  are  a first  step  in  the  process  of  a full  analysis  of  the  ASRM  field  joints. 
The  first  step  involves  the  analysis  of  both  the  forward  and  aft  motor  field  joints  at  the  0 and  19  second 
motor  burn  back  times.  The  zero  second  bum  back  time  has  the  potential  for  causing  the  greatest 
| possible  fluid  dynamic  induced  stresses  at  the  joints.  This  is  because  the  port  flow  Mach  number  and 

j dynamic  pressure  decrease  as  the  motor  bums,  thus  reducing  the  stresses  at  the  joints.  Initial  analyses 

| have  also  been  performed  on  the  inhibitor  stub  left  protruding  into  the  port  flow  field  at  the  field  joint 

caused  by  propellant  bum  back  at  the  19  second  bum  back  time.  The  analyses  discussed  are  for  non- 
deformed  propellant  grains.  Analyses  of  the  field  joints  deformed  from  cure  shrinkage,  thermal  cool 
down  and  gravity  loading  will  be  included  at  a later  time.  Also  a coupled  fluid  dynamic/mechanical 
stress  analysis  will  be  performed  in  conjunction  with  NASA/MSFC  mechanical  stress  analysts  in  order  to 
| assess  any  adverse  dynamic  mechanical  effects  of  the  flow  field  on  the  propellant  grain  shape. 

( The  analyses  presented  in  this  paper  have  been  performed  by  employing  a two-dimensional  axi- 

symmetric  assumption.  Fluent/BFC,  a three  dimensional  full  Navier-Stokes  flow  field  code,  has  been 
used  to  make  the  numerical  calculations.  This  code  utilizes  a staggered  grid  formulation  along  with  the 
SIMPLER  numerical  algorithm.  Wall  functions  are  used  to  determine  the  character  of  the  laminar  sub- 
layer flow  and  a standard  k—c  turbulence  model  is  used  to  close  the  fluid  dynamic  equations. 

The  analyses  performed  to  this  date  verify  that  the  ASRM  field  joint  design  operates  properly. 
The  fluid  dynamic  stresses  at  the  field  joints  are  small  due  to  the  inherent  design  of  the  field  joints.  A 
problem  observed  in  some  other  solid  rocket  motors  is  that  large  fluid  dynamic  stresses  are  generated  at 
the  motor  joint  on  the  downstream  propellant  grain  due  to  forward  facing  step  geometries.  The  design  of 
the  ASRM  field  joints  are  such  that  this  is  not  a problem  as  shown  by  the  analyses.  Also,  the  analyses  of 
the  inhibitor  stub  left  protruding  into  the  port  flow  from  normal  propellant  burn  back  show  that  more 
information  is  necessary  to  complete  these  analyses.  These  analyses  were  performed  as  parametric 
analyses  in  relation  to  the  height  of  the  inhibitor  stub  left  protruding  into  the  motor  port.  A better 
estimate  of  the  amount  of  the  inhibitor  stub  remaining  at  later  bum  times  must  be  determined  since  the 
height  which  the  inhibitor  stub  protrudes  into  the  port  flow  drastically  affects  the  fluid  dynamic  induced 
stresses  on  the  propellant  grain  at  the  field  joints. 
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OBJECTIVES 


1)  DETERMINE  SLOT/PORT  FLOW  INTERACTIONS  FOR  ASRM  FWD  AND  AFT  FIELD 
JOINT  DESIGNS 

2)  PERFORM  PRELIMINARY  CFD  ANALYSES  OF  INITIAL  GRAIN  CONFIGURATIONS  AT 
THE  FORWARD  AND  AFT  FIELD  JOINTS  TO  DETERMINE  PROPELLANT  GRAIN 
PRESSURE  LOADS  AND  IDENTIFY  POTENTIAL  EARLY  DESIGN  PROBLEMS 


ERC,  Inc.  4/29/92 


CFD  METHODOLOGY 


- GOVERNING  EQUATIONS  ARE  THE  3-D  ENSEMBLE-AVERAGED  NAVIER  STOKES  EQUATIONS  IN 
CONSERVATION  FORM 

- CLOSURE  OF  THE  EQUATIONS  BY  THE  STANDARD  TWO-EQUATION  k-e  MODEL  OF 
TURBULENCE 

- WALL  FUNCTIONS  USED  TO  DETERMINE  NEAR  WALL  GRADIENTS 

- DISCRETIZATION  METHOD 

• GOVERNING  EQUATIONS  ARE  WRITTEN  IN  COMPONENT  FORM  USING  CONTRAVARIANT 
VELOCITY  COMPONENTS 

• THIS  ALLOWS  THE  USE  OF  A BOUNDARY  FITTED  CURVILINEAR  COORDINATE  SYSTEM 

• NUMERICAL  METHOD  IS  FINITE  VOLUME  BASED 

• STAGGERED  GRID  STORAGE  SYSTEM  IS  USED 

• CONVECTION  AND  DIFFUSION  FLUXES  ARE  APPROXIMATED  USING  A POWER-LAW  SCHEME 

• TIME  DERIVATIVES  ARE  CALCULATED  USING  A FULLY  IMPLICIT  FIRST  ORDER  SCHEME 

- PRESSURE-VELOCITY  COUPLING  IS  ACCOMPLISHED  BY  USING  THE  SIMPLER  ALGORITHM 

- SOLVER  USES  LINEARIZED  BLOCK  IMPLICIT  SCHEME 

ERC,  Inc.  4/29/92 


667 


OVERVIEW  OF  THE  GENERAL  SLOT  GEOMETRY 


Propellant  Surface 
Inhibitor  Surface 
Symmetry  Boundary 
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ASRM  FIELD  JOINT  CONFIGURATIONS  ANALYZED 


• 0 SECOND  BURN  TIME  MOTOR  CONFIGURATION 

- FORWARD  SLOT 

- AFT  SLOT 

• 19  SECOND  BURN  TIME  MOTOR  CONFIGURATION 

- FORWARD  SLOT 

INHIBITOR  STUB  HEIGHT,  3.9  INCHES 
INHIBITOR  STUB  HEIGHT,  0 INCHES 

- AFT  SLOT 

INHIBITOR  STUB  HEIGHT,  3.9  INCHES 
INHIBITOR  STUB  HEIGHT,  0 INCHES 


ERC,  Inc.  4/29/92 
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AFT  SLOT 

INLET  STATIC  PRESSURE 

821.9  psia 

AVERAGE  PORT  VELOCITY 

177.5  ft/s 

STAGNATION  TEMPERATURE  AT  THE  INLET 

6345  °R 

RATIO  OF  SPECIFIC  HEATS 

1.128 

PROPELLANT  INJECTION  VELOCITY 

13.467  ft/s 

MOLECULAR  WEIGHT 

29.489 

MASS  FLOW  RATE  (INLET) 

8682  lbm/s 

CFD  CALCULATED  MASS  FLOW  RATE  (INLET) 

CADUIADn  Cl  AT 

8562  lbm/s 

FORWARD  SLOT 
INLET  STATIC  PRESSURE 

855.9  psia 

AVERAGE  PORT  VELOCITY 

877.5  ft/s 

STAGNATION  TEMPERATURE  AT  THE  INLET 

6345  °R 

RATIO  OF  SPECIFIC  HEATS 

1.128 

PROPELLANT  INJECTION  VELOCITY 

9.9837  ft/s 

MOLECULAR  WEIGHT 

29.489 

MASS  FLOW  RATE  (INLET) 

6178  lbm/s 

CFD  CALCULATED  MASS  FLOW  RATE  (INLET) 

6103  .bm/s 

ERC,  Inc.  4/29/92 
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ASRM  Fwd  Slot  Undeformed  Grain  Port  Pressures  On  the 
Surface  and  at  the  Motor  Centerline 


870 
860 
850 
840 
830 

Pressure  (psia)  820 
810 
800 
790 
780 
770 

20  21  22  23  24  25  26  27  28  29 

Axial  Distance  (inches) 

° Wall  ♦ Centerline 
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ASRM  Aft  Slot  Undeformed  Grain  Port  Pressures  On  the 
Surface  and  at  the  Motor  Centerline 


Pressure  (psia) 


Axial  Distance  (inches) 


ASRM  MOTOR  FIELD  JOINT  BOUNDARY  CONDITIONS 
19  SECOND  BURN  TIME  CONFIGURATION 


AFT  SLOT 


INLET  STATIC  PRESSURE  861.3  psla 

AVERAGE  PORT  VELOCITY  746.25  ft/s 

STAGNATION  TEMPERATURE  AT  THE  INLET  6317.6  °R 

RATIO  OF  SPECIFIC  HEATS  1 .1 28 

PROPELLANT  INJECTION  VELOCITY  1 0.1 35  ft/s 

MOLECULAR  WEIGHT  29.295 

MASS  FLOW  RATE  (INLET)  8846  Ibm/s 

CFD  CALCULATED  MASS  FLOW  RATE  (INLET)  8824  Ibm/s 

FORWARD  SLOT 

INLET  STATIC  PRESSURE  886.0  psla 

AVERAGE  PORT  VELOCITY  521 .2  ft/s 

STAGNATION  TEMPERATURE  AT  THE  INLET  6317.6  °R 

RATIO  OF  SPECIFIC  HEATS  1.128 

PROPELLANT  INJECTION  VELOCITY  9.956  ft/s 

MOLECULAR  WEIGHT  29.295 

MASS  FLOW  RATE  (INLET)  5963  Ibm/s 

CFD  CALCULATED  MASS  FLOW  RATE  (INLET)  5944  .bm/s 


ERC,  Inc.  4/29/92 


679 


Caaaaai”“l<ln"1 

SBI^B^I^SSSSSSSSSbbbB  ■ i ■ ■ ■ 

■SSSSSSSSaaiHMiHiKi'111 

■■■■■■■■■■IllllllllUlllia 

■■■■■■■■■■1111111111111 M 


MB— 


pSS^jfgifj 

wS 

tan 


mini 


flSRH  AFT  SLOT  AT  19  SEC0N0S/3. 87"  1 

Inhibitor  Hmlaht I 

FLUENT/BFC  V3.  0P 

Finite-Difference  Grid 

2D  Domo i n 

68 


681 


KEY 


Mini  nun  * 
-3.48E*01 

-3.48E.0l 

-1.25E+01 

9. BSE  *00 

3. 2 2 C *0 1 

5.  ^ 5£  *0 1 

7.68E*01 

9.91E.01 

1.21E.02 

1.44E.02 

1 . 66E *02 

1 . 88E.02 

| 2. 1 t £.02 

2. 33E*02 

| 2. 5SE  *02 

2.78E.02 

3. 0BE.02 

Ho x I nun  ■ 

4.30E.02 


mi 


lai 


■ s"  jrsm  tit  n ml 


-n  - **>+***'•*  /f  ; 

> s 4 % t i « //yj 


T*r»  nr  ms 
• * * * <■  +■  — 
''•••*+  #»  ^ ... 
’•••*+  / / ✓ ^ 

• < * > / s S S S 

**  • I * / S S S * 
"U/  / / / ✓ ✓ 
1 //////  / / / 
till  l I l i i i 
J I I l \ \ \ 


-5r^i'K  - 

• ^ Vn  s 

l-V^v'-v 


ll»SI 


■ • - •■  . - 


RSRM  AFT  SLOT/REflOJUSTEO  GRIP  HALF-HEIGHT  I NH I 
V © I o c i tv  Vectors  colored  bv  U~VELOCITV|  20  Dorcoin 

Steodv  Stote 


682 


683 


685 


ASRM  Fwd  Slot  Undeformed  Grain  Port  Pressures  On  the 
Surface  and  at  the  Motor  Centerline 

Upstream  Edge  of  Slot 

890 
880 
870 

Pressure  (psia)  860 
850 
840 
830 

800  850  900  950  1000  1050  1100 

Axial  Distance  (inches) 


Wall 


Centerline 


ASRM  Aft  Slot  Undeformed  Grain  Port  Pressures  On  the 
Surface  and  at  the  Motor  Centerline 


ASRM  Aft  Slot  Undeformed  Grain  Port  Pressures  On  the 
Surface  and  at  the  Motor  Cenerline 


870 

860 

850 

840 

830 

g Pressure  (psia)  820 
810 
800 
790 
780 
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1250 


1300 


1350  1400  1450 

Axial  Distance  (inches) 


1500 


155 


CONCLUSIONS 


1)  CFD  ANALYSES  HAVE  BEEN  COMPLETED  FOR  THE  AFT  AND  FORWARD  SLOTS  AT  0 AND  19 
SECOND  BURN  TIMES. 

2)  THE  PRESSURE  LOADS  AT  0 SECOND  BURN  TIME  ARE  SMALL. 

3)  THE  PRESSURE  LOADS  ON  THE  PROPELLANT  GRAIN  AT  THE  MOTOR  JOINTS  AT  1 9 SECONDS 
BURN  TIME  IS  SIGNIFICANTLY  AFFECTED  BY  THE  INHIBITOR  HEIGHT  AND  ORIENTATION. 

4)  THE  SLOT  REGION  ANALYSES  ARE  BEING  EXTENDED  TO  INCLUDE  ACTUAL  DEFORMED  GRAIN 
AND  ERODED  INHIBITOR  GEOMETRIES. 

5)  INTERACTIVE  CFD/STRUCTURAL  ANALYSES  ARE  REQUIRED  TO  PROVIDE  REALISTIC 
ASSESSMENT  OF  THE  SLOT/PORT  FLOW  INTERACTIONS  AND  RESULTING  PROPELLANT 
LOADS. 


ERC,  Inc.  4/29/92 
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Effect  of  Including  Variable  Gas 
Properties  and  Entrained  Particles  in 
the  Flow  Analysis  of  the  ASRM  Nozzle 


Curtis  D.  Clayton,  Ph.D. 
Aerojet  ASRM  Division 
luka,  Mississippi 


ABSTRACT 

CFD  analyses  of  solid  rocket  motors  typically  use  constant  fluid  properties 
throughout  the  flow  domain.  While  this  may  be  an  acceptable  approximation 
inside  the  motor  chamber,  it  is  probably  not  a good  approach  for  the  expansion 
that  occurs  in  the  nozzle. 

As  the  flow  expands  from  900  psi  chamber  pressure,  the  temperature  de- 
creases by  35%,  viscosity  and  thermal  conductivity  are  reduced  by  similar 
amounts  (25%),  and  the  specific  heat  (Cp)  and  specific  heat  ratio  (y)  change  by 
4%  and  1%,  respectively.  While  the  change  in  y appears  to  be  small,  its  effect 
is  significant  because  of  its  use  as  an  exponent  in  the  isentropic  expansion 
equations. 

The  objective  of  this  study  is  to  determine  the  effect  of  using  constant  gas  prop- 
erties for  the  analysis  of  the  ASRM  nozzle  and  to  gain  more  understanding  con- 
cerning those  types  of  analysis  which  might  require  this  additional  complexity. 

Kinetics  data  for  viscosity,  thermal  conductivity,  specific  heat  (Cp),  and  specific 
heat  ratio  (y)  are  extracted  from  the  Solid  Propellant  Rocket  Motor  Performance 
Prediction  Computer  Program  (SPP)  and  tabulated  as  a function  of  tempera- 
ture. These  tables  are  added  to  the  Aerovisc  CFD  code  in  place  of  the  constant 
gas  property  values. 

The  results  of  a CFD  analysis  of  the  ASRM  48"  motor  with  constant  gas  proper- 
ties will  be  compared  with  an  analysis  which  uses  variable  gas  properties. 
Mach  number,  surface  pressure, and  torque  plots  will  be  presented.  A full  scale 
ASRM  analysis  using  SPP  with  and  without  particle  flow  will  also  be  presented. 
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1992  Workshop  for  CFD  Applications  in  Rocket  Motors 

LOCKHEED  • AEROJET  • RUST 

Overview 

□ CFD  analyses  frequently  assume  that  the  fluid 
behaves  as  an  ideal  gas  with  constant 
properties. 

□ Combustion  gases  in  solid  rocket  motors  are 
not  ideal. 

■ Chemical  and  phase  changes 

■ Entrained  particles  in  the  flow 

■ Extreme  temperature  changes  effect  gas 
properties 

□ This  study  compares  several  analyses  with 
and  without  these  effects  to  determine  their 
impact  on  the  results. 
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LOCKHEED  • AEROJET  • RUST 

Variable  Gas  Properties 

□ Performed  Aerovisc  CFD  analysis  of  48"  motor 
with  a scaled  ASRM  nozzle. 

□ Case  1 was  run  with  constant  gas  properties 
typical  of  chamber  conditions. 

a Case  2 used  local  temperatures  to  determine: 

■ Specific  heats  - Cp  and  Cv 

■ Viscosity,  \i 

■ Thermal  heat  transfer  coefficient,  K 

□ Values  were  obtained  from  the  SPP  kinetics 
module. 


Cp  (ft-lbf/lmb-R) 
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1992  Workshop  for  CFD  Applications  in  Rocket  Motors 

Effect  of  Temperature  Changes  on 
Thermodynamic  Properties 


Property 


Change  for  2000  R 
Temperature  Change 
(4000-6000  R) 


Cp 

Cv 

Conductivity 

Viscosity 


3.5  % 

3.6  % 
33% 
27% 
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Mach  Number 


1992  Workshop  for  CFD  Applications  in  Rocket^  Motors 

48"  Motor  Nozzle  Exit  Plane 
Mach  Number 


Pressure  (psia) 


Asibvy 

LOCKHEED  * AEROJET  • RUST 


1992  Workshop  for  CFD  Applications  in  Rocket^  Motors 

48"  Motor  Nozzle  Pressure  Profile 

Along  Wall 


Variable  Gas  Properties 


Constant  Gas  Properties 


Aerovisc  CFD  Analysis 


10  15 

Axial  Distance  (inches) 
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Effect  of  Pressure  Differences 

□ Fictitious  torque  created  by  applying  different 
pressure  profiles  to  different  halves  of  the 
nozzle. 

□ Net  torque  2.32  K in-lbf 

■ 1 .3%  of  torque  due  to  variable  properties 

■ 342  K in-lbf  if  scaled  to  ASRM 

□ This  would  be  a significant  error  if  it  were 
present  in  the  gimbaled  nozzle  analysis. 


Torque  per  Unit  Length  (Ibf) 


AsBfW 
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48"  Motor  Nozzle  Torque 
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1992  Workshop  for  CFD  Applications  in  Rocket  Motors 


LOCKHEED  • AEROJET  • RUST 

Entrained  Particle  Flow 


□ Compared  SPP  analysis  with  and  without 
particles  in  full  scale  ASRM. 

□ SPP  - Solid  Rocket  Motor  Performance 
Program 

■ Industry  standard  for  motor  performance 

■ Inviscid,  axisymmetric,  real  gas  flow 

■ Method  of  Characteristics 

i 

□ Three  particle  groups  with  diameters  of  6.3|x, 
11.7|i,  and  21.8(l 


Mach  Number 
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48"  Motor  Nozzle  Exit  Plane 
Mach  Number 


CFD  with  Variable  Gas  Properties 
y ^ CFD  with  Constant  Gas  Properties 


SPP  without  Particles 


\\^ 


SPP  with  Particles 


Radial  Distance  (inches) 


Pressure  (psia) 
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Effect  of  Particles  on  Wall  Pressure 
in  Full  Scale  ASRM  Nozzle 


SPP  With  Particles 


SPP  Without  Particles 


Axial  Distance  (inches) 


Torque  per  unit  Length  (K  Ibf) 
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Effect  of  Particles  on  Torque 
in  Full  Scale  ASRM  Nozzle 
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Effect  of  Pressure  Differences 

a Fictitious  torque  created  by  applying  different 
pressure  profiles  to  different  halves  of  the 
nozzle. 

a Net  torque  142  K in-lbf 

s ■ 0.39%  of  torque  due  to  particles 

■ Net  torque  would  have  been  greater  if 
pressure  profiles  had  not  crossed 

□ This  would  be  a significant  error  if  it  were 
present  in  the  gimbaled  nozzle  analysis. 


1992  Workshop  for  CFD  Applications  in  Rocket  Motors 


Torque  Differences 
With  and  Without  Particles 
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Flow  Field  Effects 

□ Particles  slow  the  flow  along  the  center  line. 

■ 33%  mass  fraction 

■ Less  actual  gas  per  unit  volume  to  expand 

■ Velocity  difference  adds  to  drag 

□ Using  an  effective  R does  not  correct  for  these 
problems  because  the  particles  are  not 
uniformly  distributed. 

□ When  the  nozzle  is  gimbaled  the  low-flow, 
particle  entrained  region  comes  closer  to  the 
wall. 


Mach  Number 


I 
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ASRM  Nozzle  Exit  Plane 
Mach  Number 


With  Particles 


Radial  Distance  (inches) 


1992  Workshop  for  CFD  Applications  in  Rocket^  Motors 

ASRM  Nozzle 


LOCKHEED  • AEROJET  • RUST 


Particle  laden  flow  retards  velocity  in  center. 

Flow  along  wall  is  nearly  equal  to  non-particle  values. 
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Particles  influence  flow  near  nozzle  wall. 

Flow  decreases  along  lower  wall;  increases  on  upper 
wall. 
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Summary 

□ Actual  gas  properties  and  characteristics  do 
effect  the  flow  field. 

■ Especially  compressible  flows 

□ Particles  have  a significant  impact  on  the  flow 
field. 

□ Torque  calculations  magnify  these  effects. 
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A TWO-PHASE  RESTRICTED  EQUILIBRIUM  MODEL  FOR  COMBUSTION 
OF  METAUZED  SOUD  PROPELLANTS* 

J.  S.  Sabnis*,  F.  J.  de  Jong  and  H.  J.  Gibeling 
Scientific  Research  Associates,  Inc. 

Glastonbury,  Connecticut 

ABSTRACT 

An  Eulerian-Lagrangian  two-phase  approach  has  been  adopted  to  model  the  multi-phase 
reacting  internal  flow  in  a solid  rocket  with  a metalized  propellant.  An  Eulerian  description  has  been 
used  to  analyze  the  motion  of  the  continuous  phase  which  includes  the  gas  as  well  as  .the  small 
(micron-sized)  particulates,  while  a Laqrangian  description  is  used  for  the  analysis  of  the  discrete 
phase  which  consists  of  the  larger  particulates  in  the  motor  chamber.  The  particulates  consist  of 
Al  and  Al2  O,  such  that  the  particulate  composition  is  100%  Al  at  injection  from  the  propellant 
surface  with  Al3  O,  fraction  increasing  due  to  combustion  along  the  particle  trajectory.  An  empirical 
model  is  used  to  compute  the  combustion  rate  for  agglomerates  while  the  continuous  phase 
chemistry  is  treated  using  chemical  equilibrium.  The  computer  code  was  used  to  simulate  the 
reacting  flow  in  a solid  rocket  motor  with  an  AP/HTPB/AI  propellant.  The  computed  results  show  the 
existence  of  an  extended  combustion  zone  in  the  chamber  rather  than  a thin  reaction  region.  The 
presence  of  the  extended  combustion  zone  results  in  the  chamber  flow  field  and  chemical  being  far 
from  isothermal  (as  would  be  predicted  by  a surface  combustion  assumption).  The  temperature  in 
the  chamber  increases  from  about  2600  K at  the  propellant  surface  to  about  3350  K in  the  core. 
Similarly  the  chemical  composition  and  the  density  of  the  propellant  gas  also  show  spatially 
non-uniform  distribution  in  the  chamber.  The  analysis  developed  under  the  present  effort  provides  a 
more  sophisticated  tool  for  solid  rocket  internal  flow  predictions  than  is  presently  available,  and  can 
be  useful  in  studying  apparent  anomalies  and  improving  the  simple  correlations  currently  in  use.  The 
code  can  be  used  in  the  analysis  of  combustion  efficiency,  thermal  load  in  the  internal  insulation, 
plume  radiation,  etc. 


* This  work  was  supported  by  Phillips  Laboratory,  Edwards  AFB,  under  Contract  F0461 1-86-C-0096 

* Currently  at  United  Technologies  Research  Center,  East  Hartford,  CT 
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PRECEDING  PAGE  BLANK  NOT  FILMED 


OBJECTIVE 


DEVELOPMENT  OF  A PREDICTIVE  TOOL  FOR  REACTING  MULTI-PHASE  FLOW 
IN  SOLID  ROCKETS  WITH  METALIZED  PROPELLANTS 


MOTIVATION 


• MOTOR  PERFORMANCE  AND  DESIGN 

• STABILITY  ANALYSIS 
- EROSIVE  BURNING  MODELING 

• NOZZLE  PERFORMANCE  PREDICTON 

• INSULATOR  DESIGN 


PREVIOUS  ANALYSES 


• POTENTIAL  FLOW  MODELS 

• ROTATIONAL  INVISCID  FLOW  MODEL 

• INADEQUATE  WHEN  VISCOUS  EFFECTS  BECOME  SIGNIFICANT 

• INVISCID  FLOW  MODELS  COUPLED  WITH 
ITERATIVE  BOUNDARY  LAYER  CORRECTIONS 

- REQUIRE  DIVISION  OF  FLOW  INTO  INVISCID  CORE  AND  BOUNDARY  LAYER 

• NAVIER-STOKES  ANALYSES 

• TWO-PHASE  NAVIER-STOKES  ANALYSES 


Scfentfc 

Rtsesrch 

As*ocmt*s 
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CHARACTERISTICS  OF  METALIZED  SOLID 
PROPELLANT  COMBUSTION 


• MULTI-PHASE  REACTING  FLOW  WITH  POLY-DISPERSED  PARTICLES  WHICH 
UNDERGO  CHANGES  IN  SIZE,  TEMPERATURE  AND  COMPOSITION 

• DISPARATE  TIME  SCALES  ASSOCIATED  WITH  THE  DROPLET  COMBUSTION 
AND  THE  GAS  PHASE  CHEMISTRY 

• SPATIAL  VARIATION  OF  GAS  COMPOSITION  DUE  TO 

- SPECIES  TRANSPORT 

• COMBUSTION  OF  ALUMINUM  DROPLETS 

- CHANGE  OF  THERMODYNAMIC  STATE 

• GAS  COMPOSITION  DETERMINES  MIXTURE  MOLECULAR  WEIGHT, 

SPECIFIC  HEAT  ETC. 


<2 
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TWO  PHASE  REACTING  FLOW  MODEL  FOR 
SOLID  ROCKET  INTERNAL  FLOWS 


CONTINUOUS  PHASE  CONSISTING  OF  PRODUCTS  OF 
COMBUSTION  FROM  AP,  BINDER  AND  REACTED  ALUMINUM. 

> DISCRETE  PHASE  CONSISTING  OF  DROPLETS  CONTAINING 
UNREACTED  ALUMINUM  AND  Afe03  CAPS. 

• KINETIC  TIME  SCALES  FOR  CONTINUOUS  PHASE  REACTIONS 
SIGNIFICANTLY  SMALLER  THAN  FLUID  DYNAMIC  TIME  SCALES. 
HENCE  CHEMICAL  EQUILIBRIUM  ASSUMED  FOR  CONTINUOUS 
PHASE  ANALYSIS. 

> COMBUSTION  OF  ALUMINUM  DROPLET  LIMITED  BY  AVAILABILITY 
OF  OXIDIZING  SPECIES  AT  THE  DROPLET.  HENCE  MASS  TRANSPORT 
CONTROLLED  COMBUSTION  OF  ALUMINUM  DROPLET  ASSUMED. 

1 AN  EULERIAN-LAGRANGIAN  ANALYSIS  ADOPTED  TO  SIMULATE  THE 
MULTI-PHASE  REACTING  FLOW. 


CONTINUOUS  PHASE  ANALYSIS 


CONTINUITY  EQUATION 


+ V-(apt/)  = mv 


MOMENTUM  EQUATION 


+ V-(apUU)  = V{ap)  + Vcn  + mvUp+FD 


ENERGY  EQUATION 


Scientific 

Research 


+ V-(apUh)  = + af  + Vq  + qv 


-UR.FD+mv[hv+l(JFrUR>j 
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CONTINUOUS  PHASE  ENERGY  EQUATION 


MIXTURE  ENTHALPY 


I Yfh, 
/-I 


0 5 / 

ht  = hi  + £ 

/=i  J 


HEAT  FLUX  VECTOR 

q=  kVT  - 2h,pDVY? 

1*1 

IF  TURBULENT  AND  LAMINAR  LEWIS  NUMBERS  ARE 
ASSUMED  TO  BE  UNITY,  THIS  CAN  BE  SIMPLIFIED  TO 


CONTINUOUS  PHASE  MASS  TRANSPORT  ANALYSIS 

<?  (p  Vs  ) 

-L-JJ  * -v  (pD,VVf  ) * »,♦  mVi 


NOTES: 


DEFINE 


1.  2>,  = 0 

/«1 


2.  my  m 0 FOR  l>2  (l-1=>AI;l-2=*A^03) 


ow  = MASS  FRACTION  OF  ELEMENT  k IN  SPECIES 
Yk  = IawY/S  = MASS  FRACTION  OF  ELEMENT  k IN  GAS 
aM  - ak2  - 0 FOR  k>2  (k-1=>AI;  k-2=»Q 
ASSUME 

D,  = D\  / = 1,2,...n 


CONTINUOUS  PHASE  MASS  TRANSPORT  ANALYSIS  (cont'd) 


d(pYk)  n 

-^YT  + V (pUY^  = V(pDVy*>  + 


BOUNDARY  CONDITIONS 

• PROPELLANT  SURFACE:  Yk  SPECIFIED  BY  PROPELLANT 

COMPOSITION 

3Yk 

• INERT  SURFACE:  — — = 0 

dn 

• IF  PROPELLANT  COMPOSITION  IS  UNIFORM  THEN  ELEMENTAL 

TRANSPORT  EQUATIONS  FOR  k = 3 t ARE  PROPORTIONAL 

TO  EACH  OTHER  AND  NEED  NOT  BE  SOLVEDI 

• IF  I1V|  ocmy2  THEN  ELEMENTAL  TRANSPORT  EQUATIONS  FOR  Al 
AND  0 ARE  PROPORTIONAL  AND  ONLY  ONE  NEEDS  TO  BE  SOLVED. 


CONTINUOUS  PHASE  CHEMISTRY  ANALYSIS 


ELEMENTAL  MASS  FRACTIONS  AND  TWO  STATE  VARIABLES 
DEFINE  THE  THERMODYNAMIC  STATE. 

COMPUTE  SPECIES  MASS  FRACTIONS  USING  SUITABLE 
CURVE  FITS  OBTAINED  FROM  EQUILBRIUM  CODE. 

COMPUTE  TEMPERATURE  FROM  MIXTURE  COMPOSITION 
AND  ENTHALPY  USING  POLYNOMIAL  COEFFICIENTS 
(ITERATION  REQUIRED). 

COMPUTE  MIXTURE  MOLECULAR  WEIGHT  AND  SPECIFIC 
HEAT. 


SOntHc 

R0S*MTCh 

Associates 
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DISCRETE  PHASE  ANALYSIS 


• COMPUTATIONAL  PARTICLES  USED  TO  REPRESENT 
COLLECTION  OF  DROPLETS  CONTAINING  ALUMINUM  AND 
ALUMINUM  OXIDE 

• DISTRIBUTION  OF  PARTICLES  IN  THE  DOMAIN  COMPUTED 
USING  LAGRANGIAN  ANALYSIS  IN  COORDINATE  SPACE 

• SOURCE  TERMS  FOR  EULERIAN  ANALYSIS  COMPUTED 
FROM  THE  COMBUSTION  RATE,  DRAG  FORCE  AND 
HEAT  TRANSFER  FOR  THE  PARTICLES 


DISCRETE  PHASE  ANALYSIS 


• EQUATION  OF  MOTION  FOR  PARTICLE 


Fp  = ftlXp 


Scientific 
Research 
Associates 


COORDINATE  TRANSFORMATION 


= hL 


y1  = y,(x1,x2,x3) 

Yp  = JXp  where  J 


INTEGRATION  YIELDS 

AYP  = + AtXp\ 

p 2 m P|f0 
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DISCRETE  PHASE  DROPLET  COMBUSTION  MODEL 


IDEALIZED  ANALYSIS  DIFFUSION  CONTROLLED  COMBUSTION 
OF  DROPLET  WITH  SURFACE  REACTION  YIELDS 


WHERE 


mb  =2wDppDln(1+*  Y0-) 

X - STOICHIOMETRIC  FUEL  TO  OXIDIZER  MASS  RATH 
Y0>-  = MASS  FRACTION  OF  OXIDIZER  IN  FAR  FIEL 


FOR  MASS  TRANSPORT  CONTROLLED  COMBUSTION  OF 
ALUMINUM  DROPLET  CALCULATE  BURNING  RATE  FROM 


• _ * k -?~n 

mAI~  2PAI  nDP 


- REDUCES  TO  RESULTS  OF  THE  IDEALIZED  ANALYSIS 
FOR  n=2  AND  APPROPRIATE  EXPRESSION  FOR  k. 


DISCRETE  PHASE  DROPLET  COMBUSTION  MODEL  (Cont'd) 


Scientific 
Ammic/i 
Associate* 


-REDUCES  TO  HERMSEN  MODEL  WITH 

n =18 

k =8.3314x10"*  •A°sp^27Rk 
Rk=  2.7 

Ak=  100ZX,  ; i = C02,  H20,02,0H,0 

-CAN  ACCOUNT  FOR  EFFECTS  OF  FORCED  CONVECTION, 
DROPLET-DROPLET  INTERACTION,  VARIATION  OF  OXIDIZER 
CONCENTRATION  ALONG  TRAJECTORY  ETC.  ON  BURNING 
RATE  WITH  REDUCED  EMPIRICISM. 

FOR  EXAMPLE 

k«1  + 0.24  Re^2Sc^3 


n 

1 E*,X,S 
W-1 


IN  GENERAL 


k=k(.D,  Sh,  a,  */,  Y,s ) 
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APPLICATION  CASE  STUDY 


• GEOMETRY 


PROPELLANT  DATA 
COMPOSITON 
AP 

BINDER  (HTPB) 

A I 

m 

DENSITY 
BURN  RATE 


71.0% 

14.0% 

lUfMt 


1794.6  kg/m3  (112.0  lbm/ft3) 
9.0678  x 10-3  m/s  (0.357  in/s) 


APPLICATION  CASE  STUDY(Continued ) 


CASE  I - SURFACE  COMBUSTION  SIMULATION 

• ALL  ALUMINUM  ASSUMED  TO  BURN  AT  PROPELLANT  SURFACE 

• FLAME  TEMPERATURE  - 3435  K 

• 20%  AI4O3  ASSUMED  TO  BE  IN  CAPS 

- GAS  PHASE  ELEMENTAL  ALUMINUM  MASS  FRACTION 
GIVEN  BY  (0.8  x 0.15)  / (1  - 02  x 0.15  X 102  / 54)  » 0.1272 

CASE  II  - DISTRIBUTED  COMBUSTION  SIMULATION 

• ALL  ALUMUNUM  INJECTED  AS  DROPLETS  IN  DISCRETE  PHASE 

• FLAME  TEMPERATURE  a 25924  K 

• LOG-NORMAL  SIZE  DISTRIBUTION  WITH  MEAN  DIAMETER  = 150|im 
AND  LOG,0  STANDARD  DEVIATION  = 0.2 
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PARTICLES  ASSUMED  TO  RETAIN  20%  OF  REACTED  ALUMINUM 
IN  FORM  OF  Al/X,  CAPS 
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ALUMINUM  MASS  FRACTION  CONTOURS 
(Y^s  17.0%,  Ymi„  = 0.0%,  AY  = 1.0%) 


DISTRIBUTED  COMBUSTION 


SURFACE  COMBUSTION 


SUMMARY 


• A TWO-PHASE  DISTRIBUTED  COMBUSTION  MODEL  DEVELOPED 
TO  SIMULATE  COMBUSTION  OF  METALIZED  SOLID  PROPELLANTS 

• CALCULATED  RESULTS  SHOW  EXISTENCE  OF  AN  EXTENDED 
COMBUSTION  REGION  IN  THE  MOTOR  CHAMBER 

- SIGNIFICANT  SPATIAL  VARIATION  IN  TEMPERATURE,  COMPOSITION 
AND  DENSITY  IN  THE  CHAMBER 

• EXPERIMENTAL  DATA  NEEDED  FOR  INmAL  PARTICLE  SIZE 
DISTRIBUTION  AND  FRACTION  OF  METAL  THAT  BURNS  AT 
SURFACE  FOR  FURTHER  CODE  VALIDATION 

• CODE  CAN  BE  EFFECTIVELY  USED  IN  PARAMETRIC 
STUDIES  AT  PRESENT 

• PRESENT  APPROACH  CAN  BE  READILY  MODIFIED  TO  STUDY 

EFFECTS  SUCH  AS  RADIATION  AND  PARTICLE  SIZE  CHANGES 
DUE  TO  BREAKUP  AND  COALESCENCE 
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